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INTRODUCTION 


The  quantum  theory  of  atoms,  molecules  and  solids,  in  its 


non-relativistic  formulation,  stems  out  of  the  time  dependent 


Schrodinger  equation: 


and  from  its  time  independent  form  for  the  stationary  states: 


HT  = E f 


n n n 


The  Hamiltonian  H consists  of  the  kinetic  energy  operator  plus  the  po- 
tential energy  operator,  constructed  with  the  well-known  prescriptions 


charges  g^  and  interacting  through  coulombic  potentials,  the  Hamil- 
tonian may  be  written: 


It  is  a truism  to  say  that  the  SchrBdinger  equation  very  seldom 


may  be  solved  exactly.  In  fact,  the  vast  literature  in  this  field 
during  the  last  forty  years  contains  a large  variety  of  approximate 
methods  that  have  and  are  being  used  to  solve  this  complicated  eigen- 
value problem.  As  exact  solutions  are  not  easily  available,  it  is 
desirable  that  the  approximate  methods  provide  estimates  of  the  error, 


Thus,  for  a system  of  N particles,  with  masses  nu, 


2 


namely,  upper  and  lower  bounds  to  the  eigenvalues  E , and  criteria  of 


accuracy  for  the  approximation  to  the  eigenfunctions  Y . The  latter 
will  be  of  two  types:  local,  and  in  some  average  sense  to  be  defined 

in  each  case. 

The  approximation  methods  used  can  be  classified  as  falling 
into  two  types:  variation  and  perturbation.  Variation  methods  provide 

upper  bounds  to  eigenvalues  through  the  solution  of  secular  equations. 
They  do  provide,  too,  some  estimates  of  the  overlap  of  the  solutions 
of  the  secular  equation  with  the  exact  eigenfunctions.  By  overlap  we 
an  the  integral  fw  $ (dv).  It  is  possible  to  obtain  lower  bounds 


me 


from  variational  solutions  by  additional  computation  and  additional 

information:  in  Temple's  method^ ^ the  computation  of  the  expectation 

value  of  and  the  knowledge  of  E^,  the  first  excited  state  of  H is 

needed  to  compute  a lower  bound  to  the  ground  state,  E . Similarly, 

the  knowledge  of  E^  provides  a bound  for  the  overlap  of  with  an 

[21 

approximate  in  Eckart's  method.  . , Lower  bounds  from  variational 

methods  can  be  constructed  by  the  technique  of  intermediate  problems, 

[3] 

of  which  Weinstein,  Aronszajn,  Bazley  and  Fox:  , are  mainly  responsible. 

Perturbation  methods  can  provide  upper  and  lower  bounds  to  the 
quantities  being  approximated,  yet  this  is  not  true  for  all  pertur- 
bation schemes. 

[4] 

The  partitioning  technique,  developed  by  LBwdin  f provides  a 
synthesis  of  different  perturbation  approaches  and  relates  them  to 
variation  methods. 


Through  use  of  a bracketing  theorem,  upper  and  lower  bounds  to 
the  eigenvalues  and  local  and  norm  errors  for  the  wave  function  can  be 
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estimated.  This  has  been  extensively  discussed  by  LBwdin  in  his  series 

[4] 

on  perturbation  theory  . 

In  this  dissertation  we  try  to  explore  some  of  the  aspects  of 
partitioning  technique  combined  with  operator  inequalities  and  the 
techniques  of  inner  and  outer  projections.  A review  of  the  mathematical 
background  is  presented  in  Chapter  I,  which  indicates  tools  to  be  devel- 
oped further.  In  Chapter  II  we  explicitly  relate  the  widely  used 
Temple  formula  with  a specific  approximation  to  the  partitioning  expres- 
sions.' By  use  of  successive  partitioning  technique,  in  Chapter  III,  we 

1-rI 

establish  the  connection  with  Feenberg's  ^ perturbation  schemes.  This 
derivation  emphasizes  the  conceptual  relationships  between  the  pertur- 
bation schemes.  The  combination  of  the  inhomogeneous  equation  approach 
of  Dalgarno,  Hirschfelder  et  al[^  with  partitioning  schemes  proves  to 
be  very  fruitful  in  formulating  explicit  solutions  to  the  intermediate 
problems  in  Chapter  V.  Chapter  IV  contains  a new  derivation  of  upper 

and  lower  bounds  to  even  orders  of  perturbation  series,  obtained  origi- 

[7] 

nally  by  Frager  and  Hirschfelder  through  analogy  with  electromagnetic 

theory.  Variation  principles  for  excited  states,  via  perturbation 

theory,  are  obtained  in  this  fashion  too.  The  approach  is  geometrical 

r 0 1 

in  nature  and  follows  the  work  by  Synge  . Chapter  VI  discusses  a 
transformation  of  the  SchrBdinger  eigenvalue  problem  which,  for  one- 
electron  problem  yields  a complete,  discrete  spectrum  as  shown  by 
Hylleraas^^  and  SchrBdinger^  For  many-electron  problems,  the  trans- 
formation allows  the  use  of  the  operator  inequalities  reviewed  in 
Chapter  I thus  yielding  new  lower  bound  formulae.  In  a simple-minded 
approach  for  this  transformation  proved  to  be  very  effective^  ^ , 
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as  compared  with  the  Bazley  and  Fox  method.  Several  ways  of  computing 
lower  bounds  appear  to  be  feasible,  including  the  application  to 
molecular  systems  by  the  construction  of  an  Hermitian  positive  oper- 
ator . 

Some  numerical  aspects  of  partitioning  derivatives  of  the 
bracketing  function,  etc.,  in  connection  with  the  He  atom  are  discussed 
in  Appendix  I. 

The  discussion  will  be  based  on  the  areas  of  future  research 
and  applications  that  are  suggested  by  the  results  presented  in  this 


dissertation. 


CHAPTER  I 


MATHEMATICAL  BACKGROUND 


1.1.  Partitioning  Technique 


In  what  follows  we  present  a review  of  the  aspects  of  parti- 
tioning technique  pertinent  to  our  work.  For  a more  detailed  account 
see  the  papers  by  LBwdin  on  perturbation  theory,  which  will  be  referred 
to  as  PTI  - PTXI,^. 

Consistent  with  the  aims  of  this  study,  which  is  to  exploit 
the  operator  formalism  in  solving  the  SchrBdinger  equation,  we  will 
discuss  partitioning  in  operator  form.  The  object  of  study  is  the 
time  independent  eigenvalue  problem: 


HT  = E f 
n n n 


(1.1) 


We  seek  upper  and  lower  bounds  to  the  E^  and  well  defined  approximations 
to  the  'S  . We  shall  restrict  our  attention  to  bound  states,  although 
the  formalism  allows  the  treatment  of  scattering  states  as  well  as  time 
dependent  problems. 

The  building  blocks  of  the  theory  are  an  energy  variable  , 
and  a reference  space  of  linearly  independent  functions.  For  our  pre- 
sent purpose  it  is  enough  to  make  the  reference  space  one-dimensional, 

i.e.,  we  consider  a reference  function  Y,  assumed  to  be  normalized. 

[12] 

LBwdin  J discussed  the  multidimensional  reference  space,  and  appli- 

[13] 

cations  have  been  made  and  are  under  progress. 
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The  projection  operator  associated  to  \|r:  0 = |\|/)(\|r|  satisfies 
the  usual  relations 

02  =0;  0+  = 0;  Tr(0)  = 1 . (1.2) 

The  orthogonal  complement  to  \Jr  is  characterized  by  a projection 
operator  P with 

p2  = P.  p+  = p;  p = 1 - 0;  PO  = OP  = 0 . (1.3) 

Given  the  reference  function  and  the  energy  variable  £ , we 
introduce  the  reduced  resolvent  associated  with  (l.l): 

T = P[a0  + P(£-H)P]_1P  (1.4) 

where  a is  any  constant  different  from  zero.  The  properties  of  T 
give  a clue  to  its  nature:  it  is  the  inverse  of  £ - H in  the  subspace 
projected  by  P.  P(<£  -H)P  itself  has  no  inverse,  because  it  is  zero  in 
the  subspace  0.  By  adding  QO  to  P(£  - H)P  we  may  safely  construct  the 
inverse,  and  then  project  out  the  desired  component  in  P.  If  this  is 
to  be  consistent  at  all,  T should  be  independent  of  d: 

H = p[ao  + p(£- h)p]-1o[qo  + P(£- h)p]_1p  = 0 . (1.5) 

Notice  that  T "blows  up"  whenever  £ is  equal  to  one  of  the  eigenvalues 
of  PHP . The  properties  that  make  T interesting  are: 

0T  = TO  = 0,  pt  = TP  = T,  T = T+  . (1.6) 

If  one  multiplies  by  P the  following  expression: 

[OD  + P(£-  H)P]  [qO  + P(£-  H)P] _1  = 1 , 
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the  important  result 


p(£-  h)t  = p 


(1-7) 


follows.  It  is  customary  to  use  the  symbolic  notation 


T 


P 

e - h 


(1.8) 


although  (1.4)  is  meant  by  it. 

A trial  wave  function  associated  with  \|r  and  £ is  defined 

by: 


= \|r  + TIty  = (1  + TH)\|r  = f + 0£  . (1-9) 

This  definition  implies  the  adoption  of  intermediate  normal- 
ization: 

<1r|fg>  = (*!♦>  + = 1 + (i|f|PTH«|f)  = 1 (1.10) 

(Ye|yg)  = 1 + (*|HT2H|t)  = 1 + (<tL  |0£  ) 

If  <0£|0£)  exists,  Yg  belongs  to  the  discrete  spectrum,  otherwise,  it 
is  a scattering  state. 

We  now  consider  the  expression  (H  - S-J'Fc  • First  notice  that 
¥ is  a solution  of  the  SchrBdinger  equation  in  the  subspace  P: 

P(H  - £)l£  = P(H  - £.)t  + P(H  - <f)Tlty 


= Pity  - Pity  = 0. 


Therefore : 
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(H  - S)WL  = (0  + P)(H  - £)y£  = 0(H  - £)y 

= i|r<i|rlH  - E\y£)  = t(<t|H  | **>-  £)  . 

Denote  (\);|h|t^)  by  £^.  It  is  apparent  that  Wr  satisfies  an 
inhomogeneous  equation: 


(H  -£)w£  = (£1  -Eh  (1.11) 

If  = £ = E it  turns  out  that  is  an  eigenfunction  corresponding 
to  the  eigenvalue  E.  Obviously,  for  an  arbitrary  value  of  <? , f E, 

yet  the  relation  between  and  E is  of  particular  significance: 

E]_  = <i|f|H|y£)  = |H  + HTH  | \|r ) = f(£).  (1.12) 


Consider  the  slope  of  the  function  £ ^ = f (£) 

^ = (*|h||h|^)  = - 4|ht2h|i|/) 


a& 


(1.13) 


= - (0£|0£><O. 


Therefore,  for  values  of  £ in  the  discrete  spectrum  region. 


d £ 1S 
for  all 

See  PTI. 


negative.  This  leads  to  the  bracketing  theorem,  which  holds 
iteration  schemes  of  the  type  £^  = f(£)  with  f'(£)  <0. 


Let  E be  an  eigenvalue  of  H,  i.e.,  - f(d)  - E - 

Consider  a pair  of  values  £,  £^.  We  can  write: 


E.  = E + 6 ; £]_  = E + 6L 


(1.14) 
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6^  and  £ not  necessarily  small.  Then  it  follows: 

£l  = E + £1  = f(£)  = f(E  + €)  = f(E)  + £f'(E  + 9€) 

which  reduces  to: 

= f ' (e  + 06)  . (1.15) 

We  have  used  Lagrange's  mean  value  theorem,  and  0^  0 £ 1.  As 
f'(6)  is  negative,  £ and  £^  have  different  sign,  and  the  bracketing 
theorem  is  proved:  £ ^ and  £ bracket  one  eigenvalue  E to  H.  (See  Fig.  1.) 

Clearly  one  may  bracket  several  eigenvalues  to  H. 

An  immediate  application  of  the  bracketing  theorem  follows  from 
(1.12)  by  setting  £=  “00.  Then  = (\|/-|h|i|/)  is  an  upper  bound  to  the 
ground  state  eigenvalue  E^,  a result  that  follows  from  the  variational 
principle  otherwise. 

The  curve  = £(£)  is  interesting.  The  eigenvalues  are  the 
crossings  of  £ = f(£)  with  the  line  = £.  The  curve  has  negative 

slope  everywhere  in  the  bound  state  region  and  it  has  vertical  asymp- 
totes for  B equal  to  the  eigenvalues  of  PHP  s H,  which  we  shall  label 
by  E^.  (See  Fig.  2.) 

The  conventional  perturbation  schemes  stem  out  of  power  series 
expansions  of  (l.ll),  and  provide  approximations  to  E and  Y.  The  aim 
of  the  present  method  is  to  study  the  interval  (£,£^)  instead  of  the 
special  points  E.  This  is  significant  for  two  reasons:  Upper  and 

lower  bounds  to  the  eigenvalues  are  obtained.  Further,  the  trial  wave 
function  satisfies  the  SchrBdinger  equation  with  a local  accuracy 
(£x  - fc)  , as  indicated  by  (l.ll)  and  with  an  average  accuracy 

||(h  -£)i£||-  tx  - £ (1.16) 

which  should  be  useful  in  the  study  of  properties  other  than  energy. 


The  following  digression  is  not  out  of  place  here.  It  was 
r q i 

pointed  out  by  SyngeL  J that  many  boundary  problems  of  mathematical 
physics  can  be  stated  in  the  form  "It  is  required  to  find  the  inter- 
section of  two  orthogonal  subspaces  of  a function  space."  The  usual 
method  of  setting  up  an  infinite  process  (power  series  solutions)  has 
the  required  intersection  as  a limit.  The  individual  terms  do  not 
satisfy  in  general  the  boundary  conditions,  yet  they  are  satisfied  by 
the  limit.  An  example  of  this  would  be  Fourier-type  description  of  the 
solution  of  a boundary  value  problem.  It  is  quite  clear  that  the  par- 
titioning scheme  implies  two  orthogonal  subspaces,  in  a somehow  different 
sense  from  Synge's,  yet  with  a direct  implication  that  in  order  to 
bracket  solutions  one  must  keep  track  of  the  SchrBdinger  equation  in 
both  subspaces  0 and  P.  The  variational  principle  in  its  conventional 
form  yields  upper  bounds.  This  is  a consequence  of  working  in  a certain 
subspace  only.  In  perturbation  theory,  it  is  possible  to  give  upper 
and  lower  bounds  provided,  as  we  shall  see  later,  that  we  "keep  track" 
of  the  orthogonal  complement;  if  we  neglect  it  we  are  left  with  approxi- 
mations that  do  not  bracket  the  energy.  The  truncations  in  Chapter  III 
are  an  example  of  this  situation. 

In  PTI  LBwdin  pointed  out  the  connection  of  partitioning  tech- 
nique with  iteration  variation  procedures  in  general,  and  with  the 
Newton-Raphson  method  in  particular.  The  following  result  is  important. 

Let  £ be  an  upper  bound  to  the  ground  state  energy,  then 

£ £ E ^ £.  Clearly,  from  the  variational  principle: 

1 o 
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As  in  PTI,  it  is  eqsy  to  see  that 


£,  6 E ±e 
1 o 


Proof : 


<vIh1y£>  + (%|h|%) 
1 + <0£|0£> 


£,  +£<^|fl»e.) 
i + 


£ - 
1 + 


^ 6 . 


Q.E.D. 


The  practical 'application  of  these  results  depends  on  the  possi- 
bility of  evaluating  and  lg) • We  can't  do  it  exactly,  but  it  is 

possible  to  obtain  approximations  that  preserve  the  inequalities. 
Instead  of  computing  it  is  possible  to  evaluate  6^,  which 

still  brackets  E if  E f £. 

We  should  like  to  point  out  the  following  theorem: 


Theorem  1.1: 

If  £ ' ^ £•,  — E ft  — £.  — E.. , where  E is  the  first  eigenvalue 
1 1 o 1 1 

of  PHP,  g (£)  0£>  satisfies  the  relations: 

g(£')^  g(&)  ^g(E  ) & g(£*)  ^ g(£).  (1-18) 

d 

The  proof  is  very  simple;  it  relies  on  the  fact  that  g(£)  = - f'(£)  is 

a monotonous  function  for  -oo  £ -t-  E..  ; $im  g(£)  = 0;  J5im  g(£)  = °°. 

1 t--°°  6.-»E]_ 

t-* 

.Further  g(£)  is  non  decreasing;  hence;  g(£A)  £ g(£g)  f°r  ^ and 

this  completes  the  proof. 

Neither  g(£1)  nor  g(EQ)  are  easily  available.  Equation  (1.18) 
provides  bounds  in  terms  of  £^ ' and  £ which  ate  known.  The  problem 


of  evaluating  g(£)  itself  will  be  discussed  later.  A generalization 

of  Theorem  1.1  will  be  given  in  Chapter  V. 

The  precision  of  the  bracketing  is  determined  by  Eq.  (I.I5). 

If  |f'(E+©£,)j  <.  1 the  lower  bound  is  going  to  be  better  than  the  upper 

bound,  or  in  virtue  of  Eq.  (1.18),  if  g(£)  <■  1.  In  terms  of  Eq.  (I.9) 

and  Eq.  (1.10)  this  means  that  if  the  reference  function  \|r  contributes 

more  than  50  per  cent  towards  the  trial  wave  function  ¥ ^ > the  lower 

bound  is  better  than  the  upper  bound. 

To  conclude  the  general  remarks  on  partitioning  it  should  be 

pointed  out  that  if  the  reference  function  i|/  is  orthogonal  to  a partic 

ular  'H  , the  associated  eigenvalue,  E , is  not  bracketed  at  all,  thus 
n n 

yielding  the  only  restriction  on  i|f.  The  sharpness  of  the  bracketing 
depends  only  on  the  propinquity  of ^ to  Y . 

1.2.  Full  and  Reduced  Resolvents  Weinstein  Functions 

Given  the  Hamiltonian  H,  and  the  eigenvalue  problem  Eq.  (l.l), 
there  is  a "resolution  of  the  identity"  and  a spectral  resolution  of  H 
in  terms  of  the  eigenoperators  0^: 

HO,  = E.  0, 
k k k 


* =1°, 
k 


k 


The  operators  0^  are  self-adjoint  and  idempotent.  Furthermore  the 
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resolutions  in  Eq.  (1.19J  imply  not  only  a sum  over  the  discrete 
states,  but  the  Stieltjes  integral  over  the  continuum  spectrum.  The 
conditions  under  which  Eq.  (I.I9)  hold,  and  their  proofs  can  be  found 
in  von  Neumann^-  l^\  Riesz,  Nagy^l5\ and  Lengyel  For  our  present 

purposes,  it  is  sufficient  to  assume  that  Eq.  (1.19)  holds  for  H and 
Hq,  both  self-adjoint  operators. 

Functions  of  operators  may  be  defined  in  terms  of  their  spec- 
tral resolutions  and  we  have  then: 

f(H)  = ^ f(Ek)Qk  (1.20) 

k 

(e-H)-1  - £ (£-V'\  • 

k 

The  coefficient  of  0^  blows  up  for  £ ■»  E^.  The  significance  of  this 
follows  from  considering 

X=(£.-H)_1t  (1.21) 

for  a given  \|r. 

Jim  * _ Yn  . v 

t-En  (\|r|X)  (i|f|in)  (!-22) 

The  connection  between  the  full  resolvent  and  the  reduced 
resolvent  was  pointed  out  by  LBwdin  in  PTX.  It  arises  from  the  inhomo- 
geneous equation  (l.ll): 

(h-£)*£  = (£1  -&h 

Therefore : 

= -S.)(E-£)_1t  - (£-£i)x  . (1.23) 
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The  intermediate  normalization  implies: 


1 


(£  - £l)  <v|x> 


and  then 


The  full  resolvent  equations  (1.24)  and  (I.25)  are  associated 
with  the  inhomogeneous  equation  (1.2l).  The  inhomogeneous  equation 


approach  has  been  very  successful.  For  a review  of  its  applications 
see  [6].  It  was  used  in  connection  with  perturbation  theory.  Within 
the  same  context  we  shall  make  appeal  to  it,  but  to  solve  equations  to 
all  orders,  instead  of  a specific  order  of  the  perturbation  series. 
This  will  be  done  in  Chapters  III  and  V. 


Weinstein  functions,  and  in  general  Weinstein  determinants,  have  been 

[3] 

introduced  in  the  theory  of  "intermediate  problems"  by  Weinstein 


B of  PTX,  LBwdin  discusses  its  relation  to  variations  subject  to  a 
constraint.  We  should  expect  a connection  with  the  eigenvalues  of  PHP 
as  being  determined  by  constrained  variations.  We  shall  require  some 
use  of  the  properties  of  PHP  and  T,  therefore  it  is  convenient  to 


The  expression  (\|r|X)  is  called  the  Weinstein  function: 


W(£ ) = <i|r|(£  -H)_1|t)  • 


(1.26) 


study  W in  some  detail. 

W(\)  is  monotonously  decreasing: 
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w'(£)  = - (t  I(£-h)‘2U)  4 o 


(1.27) 


It  has  vertical  asymptotes  for  £•  = E^,  the  eigenvalues  of  H.  This  can 
be  seen  by  expanding  the  reference  function  \Jr  in  terms  of  the  eigen- 
functions to  H,  thus  obtaining: 


W 


u)  - Y ] 


00  ,2 

a. 


k=0 


(1.28) 


with  ak  = (¥k|i|f)  • 

What  are  the  zeroes  of  W(£)?  Consider  the  eigenvalue  problem  of  H = PHP. 

Hf  = 

As  py  = $ 

p(h-e)?  = 0 

(H-E)i  = \|r(i|r|H -E|5)  = i|r<i|r|  Hj  ?> 

i = (H -Ej'^^lHjf) 

Finally,  from  the  orthogonality  between  \jr  and  Y: 

0 = W(E)  (t|H|  w) 


Now,  if 

(t|H|7)  + 0, 

W(E)  = 0 . 

Therefore,  the  Weinstein  function  vanishes  for  those  eigenvalues  of  PHP 
which  are  not  simultaneously  eigenvalues  of  H.  The  zeroes  of  W(£) 
the  EH  . (See  Fig.  3-  ) 


are 
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For  £ < E 


P+1’ 


the  following  simple  estimates  of  W (£)  are  valid: 


(1.29) 


We  have  come  this  far  in  discussing  the  Weinstein  function  and  found 
that  its  zeroes  are  the  eigenvalues  of  PHP.  This  effort  becomes 
meaningful  when  we  consider  the  spectral  resolution  of  T: 

A starting  point  is  Eq.  (l.j): 

P(£  -H)T  = P (1.7) 


or  simply 

(£-  h)t  = p 

T = (8  -H)_1P  = £ (8  -Ek)_1OkP  . 
k 

Now  in  the  subspace  P,  the  resolution  of  the  identity  is  replaced  by 

p ■ £ v v ■ p5k  ■ 5k 

k=l 


where  0k  = |fk)(fk|  is  the  projection  operator  associated  with  the 
kt*1  eigenfunction  to  H. 

We  obtain  then  the  spectral  resolution: 
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Finally,  if  £<  E the  bounds 
P+1 


F A F 

E A— E 

k=l  k=l 


<5-Ek 


+ 


p - 7 \ 

k=l 


£ - E 


P+1 


= T(p) 


(1.30) 


hold,  and  are  similar  to  the  bounds  ' in.  Eq . (1.29)  for  W(£).  In  fact, 


spectral  resolutions  and  bounds  for  the  operator  T = 


o & - H 


-,  defined 


in  a similar  way,  are  of  analogous  nature.  One  has  to  replace  H by  H^, 

Y by  Y°,  E by  E°,  and  of  course  H Y°  = E°T°  must  be  satisfied, 
n nn  n onnn 

It  should  be  pointed  out  that  there  is  another  expression  for 

the  reduced  resolvent  T,  in  terms  of  the  Weinstein  function  and  the 

eigenvalues  of  H,  given  by  Lowdin  in  Appendix  A of  PTXI.  We  give  it 

here  to  illustrate  a further  relation  between  the  E.’s  and  the  E,  ’ s 

1 k 

through  W(£).  The  starting  point  is  again  Eq.  (1.7): 


P(£-  H)T  = P 

According  to  Eq.  (1.3): 

(£-  H)T  + | x(X\|r  1 HT  = 1 - |\|r)(i|r| 


(1-7) 


and 

T = (£-H)_1  - (<£  - - (f>-  I!)-1 ! \!r)  (\|r  i HT  . (1-31) 

If  we  multiply  by  (X  : 

0 = (i|/|(£  - H)  -1  - W(£)(i|r|  “ W(£)  (\Jc  | HT 

which  means 

(v|HT  = W(£)"1<1J/|(£- H)"1  - (t| 
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Replacing  in  Eq.  (I.3I): 

T = (£  - H) _1-  (£  - H)_1j\|/)(\|r|  - [W(£)]"1(5-H)“1|1(f}(tl(£ -H)"1 


Therefore, 


T = (£  -H) 


-1  (k-H)"1 

W(E) 


The  spectral  resolution  of  T in  terms  of  the  eigenfunctions  of  H is 
particularly  interesting.  For  \|r  = E a^  f we  have: 

W(£)(£  - H)"1  - (£-H)_1  0(£-H)_1 


V 3k"  V °e  V 

L>  £ -E.  L -Eo  " 
k e k,  e 


0,  00 

k e 


(£-Ek)(£-Ee) 


a,  c0  -0,  00 


Zr 1 dL.  V 1 V 1 ^ 0 0,  ^0  “0  00, 

\ k ,e  k e _ _1\  \ keekkeek 

L (£.  -Ek)(£-Ee)“  2 L L (£-Ek)(£-Ee) 
k e k e 


ke 


2 k e (£-  “ E,  )(£  - E ) 


where  the  operators  are  defined  as: 


/!,.  = a (t,  i - a.  | = -fl  . 

ke  e'k*  k'ie‘  a ek 


<A+kAe)2  - 


ke"^  keA  ke  ^ke 


(1.32) 


(1-33) 


^le^ke^^kJ^^e*2) 


(1.3^) 
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The  SL  ^ /^e  are,  except  for  a constant  factor,  projectors  and  satisfy 
the  important  relation: 


p ■ 5 E I 


ke 


(1-35) 


The  reduced  resolvent  T can  be  written  as: 


t ■ v^'1  iEE 


-a 


ke  4ke 


(£  -E,  )(£  -E) 


(1-36) 


which  should  be  compared  with  Eq.  (1.29)*  A similar  expression  should 

hold  for  T , in  terms  of  E°,  T°  and  the  analogous  JX?  . 

o e e & ke 


1.3.  Wave  and  Reaction  Operators 

These  concepts  arise  from  the  splitting  of  the  Hamiltonian 
H into  an  unperturbed  part  Hq  and  a perturbation  V:  H = Hq  + V,  where 
V is  by  no  means  necessarily  small.  Let  \|/q  be  an  eigenfunction  to 
corresponding  to  the  eigenvalue  Eq.  Then  the  Schroedinger  equation  (l.l) 


may  be  transformed,  using  intermediate  normalization  ((\|r  jy)  = l): 

(H  +V)y  = ET  (1.37) 

Eq  + <1rQ|v  y)  = E . (1.38) 

If  we  define  an  operator  W,  the  wave  operator,  by 

w^o  = y (1.39) 

and  another  operator  t,  the  reaction  operator,  by  VW  = t we  have, 

E = E + (y  |vw|\|r  ) = E + <t|t  it  k ) . (1-40) 

o Yo 1 ‘ Yo  o Yo * ‘ o ' 
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These  operators  have  been  extensively  studied  by  many  authors  in  con- 
nection with  bound  and  scattering  states,  stationary  and  time  dependent 

m 

problems.  See  reference  [18]  and  the  book  by  Kumar  [I9]  which  contains 
abundant  references  to  the  original  papers.  In  PTIV  and  PTV  they 
were  discussed  in  connection  with  infinite  order  perturbation  theory 
and  exact  selfconsistent  fields,  as  well  as  in  the  derivation  of  the 
Rayleigh-Schrbdinger  and  Brillouin-Wigner  perturbation  series.  The 
freedom  obtained  by  letting  £ become  a variable  at  our  disposal  allowed, 
in  PTIX,  the  evaluation  of  upper  bounds  of  energy  values  in  Rayleigh- 
SchrBdinger  perturbation  theory,  and  in  PTX  and  PTXI,  lower  bounds  to 
the  energy  of  ground  and  excited  states.  There  are  several  possible 
extensions  of  Eqs.  (1.39)  and  ( 1.40) ; one  may  have  an  arbitrary  refer- 
ence function;  it  may  be  a finite  linear  combination  of  eigenfunctions 
to  H , or  not  expressible  as  a finite  linear  combination  of  eigenfunc- 
tions to  Hq.  Finally,  one  may  have  a reference  space,  spanned  by  a 
set  of  such  functions. 

Consider,  for  the  time  being,  a single  reference  function. 

Define  T and  Tq  as  in  Eq.  (1.4),  and  the  wave  and  reaction  operator  as 
follows : 

W = 1 + TV  ( 1.41) 

t = VW  = V + VTV  . (1.42) 

If  we  use  the  fundamental  identity,  key  to  many  relations  in  the 
(A-B)"1  = A_1  + A_1B(A-B)_1, 


theory 


\)e  obtain 
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T = T +TVT=T( 1+VT ) 
o o o'  ' 


TV  = T t 
o 


T = T + T tT 
o o o 


and  for  the  reaction  operator: 


t = V + VT  t 
o 


which  may  be  written: 


(1.43a) 

(1.43b) 

(1.43c) 


( 1.44a) 


t = (l-VTo)_1V  = V(l-ToV)_1  . (1.44b) 

Finally,  for  the  wave  operator  we  have: 

W=(l-ToV)_1  . (I.45) 

Equation  (1.43)  Is  the  Lippmann~Schwinger  equation^^  derived  in 

[211 

this  way  by  K.  OhnoL  . If  V has  an  inverse,  we  obtain 

t = (v"1-^)"1  (1.46a) 


or  simply: 

t"1  = v_1-To  . (1.46b) 

Equations  (I.45)  and  (1.46)  will  be  the  starting  point  of  our 
discussion  in  Chapters  III  and  V.  Quite  clearly  the  evaluation  of 
Eqs.  (I.45)  and  (1.46)  implies  solving  the  Schrbdinger  equation  through 
Eq.  (I.39)  and  Eq.  (1.40).  Actually,  we  need  expressions  that  replace 
Eqs.  (I.39)  and  ( 1.40)  when  the  reference  function  is  arbitrary  and  not 
a single  eigenfunction  to  H . Multidimension  partitioning  techniques 

and  the  need  to  improve  convergence  through  Eq.  (I.I5)  establish  the  need 
of  such  an  extension. 
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The  starting  point  is  the  trial  wave  function  ^ , as  defined 
in  Eq.  ( 1.9) : 

= ♦ + THt  = (l+TH)^  (1.9) 

= [ 1 + ( l-T  V)_1T  H]\J; 

o o 

- H-nl'V-yn  ilt 

*e  =W(1«0H0)*  . (1.1*7) 

LBwdin  defined  a modified  reference  function  ijr  by: 

i = (1+T  H )>|r  , 

00 

It  is  clear  that  if  the  reference  function  ^ is  a linear 
combination  of  eigenfunctions  to  Hq,  \j?  and  a|/  ard  spanned  by  the  same 
space.  In  general,  we  have 

*£  = wt  (1.48) 

which  contains  Eq.  (I.39)  as  a very  special  case,  namely,  \|r  = \|r  and 
= £ = E.  We  need  an  expression  for  £ ^ . Hence,  we  consider 

H + HTH  = H + V + HTH  +HTV  + VTH  + VTV 
o 000  o 

= H + t + H- (T  +T  tT  )H  + H T t + tT  H 
0 o'  00  o'  o 00  00 

= H (1+T  H ) + (1+H  T ) t ( 1+T  H ) . (1.1+9) 

ov  00'  v 00'  ' 00'  ' ~ ' 

The  bracketing  function  becomes  now 


<^1  - (tlH0lt>  + <tMt> 

which,  again,  contains  Eq.  (I.38)  as  a very  special  case. 


(1.50) 
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We  conclude  by  stating  that  \|r  is  not  normalized  to  unity,  as  \J r: 

(^It)  = 1 + (t|HoTo2Ho|t)  (1.51) 

and  that  \j } depends  on  the  energy  variable,  £. 

Equation  (I.50)  is  very  general  and  has  been  the  starting  point 
fl3l 

of  several  applications  . The  applicability  depends  on  evaluating 

I 

it  exactly,  or  on  obtaining  a lower  bound  to  it,  for  the  case 

when  £ is  an  upper  bound  to  an  eigenvalue  E.  The  basis  for  this  is 
discussed  in  the  next  section. 

1.4.  Operator  Inequalities 

An  operator  A is  said  to  be  greater  or  equal  than  B,  if  for  all 
functions  belonging  to  the  intersection  of  the  domains  of  A and  B,  the 
expectation  values  of  A are  larger  or  equal  than  those  of  B,  that  is; 

Definition  1.1:  A^B  iff  for  all  f e D D D , (f  J A | f ) 5-  (f  |B  J f ) 

where  denotes  the  domain  of  A. 

An  operator  is  positive  definite,  A > 0,  iff  (f  | A | f ) > 0. 

An  operator  is  non-negative  definite,  A 0,  iff  (f  j A j f ) ^0. 

For  a finite  matrix,  a necessary  and  sufficient  condition  for 
positive  definiteness  is  that  all  the  determinants 

|a^1  ^ij  J 5 J — 1,2,  ...,1c,  1c  1,2,  ...,n 

must  be  positive.  A necessary  condition  is  that  the  determinant  | A j >0, 
and  in  general,  an  hermitian  positive  definite  operator  must  necessarily 
have  positive  eigenvalues. 
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Theorem  1.2:  (Ordering  theorem) 

If  A B,  the  eigenvalues  to  A,  a^>  and  the  eigenvalues  to  B, 
b^_,  satisfy  a^  £ b^>  in  order.  The  proof  is  simple  and  resorts  to 
the  variation  principle: 


Auk  ■ W Bvk  = Vk>  k ■ 1-2> 


a1  6 {v1|a|v1)-  (v1|b|v1)  = b:  . 


For  the  first  excited  state,  one  repeats  the  reasoning  but  using  a 
trial  function  ^ = cc^v^  + Q^v^,  orthogonal  to  u^,  that  is 
(tU)  = 1;  (tlu1)  = 0-  Then: 

a2  — (f I Alt)^  <tlBk)  = laj  2b1  + \a2\ %2  £ b£  . 

3 


The  inequality  a £ b follows  from  taking  \|r  as  i|r  = ^ CC.  u.,  normalized 

J J i=l  1 L 

and  orthogonal  to  u^  and  u^.  This  process  can  be  continued  for  higher 

eigenvalues.  The  converse  of  the  theorem  is  not  true.  If  a,  £ b , it 

iC  kC 

does  not  follow  that  A ± B. 


Theorem  1.3: 


.+ 


If  A > B,  it  follows  that  XL  A -O.  > £k  BA.  for  flf  e D f)  D . 

A D 


< f |a  1 f ) > (f  |b  |f ) ; take  f = i\g,  then 


(g|-0-+AA|g)  > (gln^AIg)  . Q.E.D 


+ 
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Inequalities  of  this  type  were  discussed  by  Rellich  and  Heinz  . 
[24] 

Kato  discusses  slightly  different  inequalities  for  linear 

operators,  which  are  useful  in  the  present  context.  These  inequalities, 
in  the  sense  of  Kato,  will  be  denoted  by  instead  of  Let  S and  T 

be  linear  operators. 
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Definition  1.2:  S T if  D 2 D and  ||  Sx//£  ||Tx|/  for  all  x e D . 

b I T 

If  S and  T are  self-adjoint,  S <£.  T is  equivalent  to  S2  < T2  in  • 

the  sense  of  Rellich.  Kato  proves  the  following  two  theorems  for  A and 

B self-adjoint: 

Theorem  1.4. : 

If  A « B and  A ^ exists,  it  follows  that  B ^ exists  and 
B-1  « A"1. 

Theorem  l.b.  Let  A > 0,  B ? 0,  A » B.  Then  AV » BV  for  0 < vi  1. 

The  concept  of  positive  definite  operators  will  prove  to  be 

very  important  in  what  follows.  We  may  qualify  somehow  definition  (l.l) 

by  introducing  the  concept  of  domain  of  positivity  R,  associated  with  a 

given  operator  A.  We  say  that  f,g  e R if  this  implies  that  (f,Ag)  0. 

T25 1 

This  was  discussed  by  Koecher  The  significance  of  this  will 

become  apparent  when  considering  inner  projections  and  expectation  values 
of  t for  V positive  definite. 

Inequalities  for  inverses  of  operators  appear  to  be  useful. 

-1  -1 

If  A >.  B > 0,  it  follows  that  A and  B exist,  moreover  we  have: 

Theorem  1.6. 

If  A > B > 0,  it  follows  that  0 < A_1  ^ B-1.  In  PTX  there  is 
a simple  proof  of  this  theorem  based  on  the  identity 

(A-B)"1  3 A-1  + A-1BA_1  + A_1B(A-B)_1BA-1  (1-52) 

which  can  be  used  in  the  following  way: 

b"1  s [a  - (a-b)]-1  = a"1  + a'1(a-b)a“1  + a_1(a-b)b_1(a-b)a"1  . 
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If  A-B  > 0,  A (A-B)A  > 0 by  Theorem  I.3.  Similarly  as  B ^ > 0, 

A (A-B)B  (A-B)A  ^ > 0.  But  then  the  theorem  follows:  B ^ > A 

Theorem  1.6  states  an  order  relation  between  A ^ and  B ^ 
starting  from  an  order  relation  for  A and  B.  The  question  arises  on 
order  relations  for  f(A)  and  f(B),  where  f(A)  is  an  arbitrary  function 
of  the  operator  A. 

Loewner^0^  studied  the  problem  of  what  are  the  functions  f(z) 
such  that  A B implies  f(A)  f(B).  It  turns  out  that  these  functions, 
called  positive  real  functions,  are  generated  by  two  simple  functions; 
cz  + Cg,  where  c^  > 0 and  -z  The  second  one  corresponds  to 
Theorem  1.7. 

[27] 

Lengyel  showed  that  if  [A,B]  = 0 and  A,B  are  positive,  it 
follows  that  AB  itself  is  positive.  This  proof  is  very  general,  and 
does  not  invoke  the  spectral  resolutions  of  the  operators.  Let's  assume 
that  A and  B are  positive,  self-adjoint  permutable  operators.  There- 
fore we  can  choose  a complete  set  of  simultaneous  eigenfunctions  to  A 
and  B.  Then  expanding  an  arbitrary  function  in  terms  of  this  set  we 
have : 


<f  | AB  | f ) = ^ ckCe^UJABlUe^ 
k,  e 

■ L CkVkb  Ae  * Z AVk a 0 • 

k,  e k 

Thus  we  proved: 

Theorem  1.7: 

If  A > 0,  BkO,  A = A+,  B = B+,  [A, B]  = 0,  if  follows  that  AB  > 0. 
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We  used  the  fact  that  a positive  operator  has  non-negative 

eigenvalues.  Lengyel's  proof  is  more  general  than  this  one. 

We  are  now  in  & position  to  examine  the  question:  if  A B, 

2 2 2 

under  what  conditions  A >.  B ? Clearly  z is  not  a monotone  matrix 

2 2 

function,  yet  we  can  state  that  a sufficient  condition  for  A 2-  B , 
if  A 2 B is  that  [A,B]  = 0. 

If  [A,B]  = 0,  (A-B)(A+B)  > 0 by  Theorem  1.8.  Then  (A-B)(A+B)  = 


Corollary : If  A >.  B,  [A,B]  = 0,  then  A^  2-  B^. 

2 4- 

Definition  1.3:  Let  Q = Q = Q , A = QAQ  is  called  the  outer  projection 

of  A,  where  A is  self-adjoint,  bounded  from  below.  Clearly,  all  functions 
in  the  orthogonal  complement  to  Q,  are  eigenfunctions  to  A associated 
with  the  eigenvalue  zero. 

Theorem  1.8: 

There  is  an  order  relation  between  the  eigenvalues  of  A and 
those  of  A,  namely  that,  a,  < a in  order,  where  Au  = a u , and 

K K tC  k K 

Au^  = a^u^.  This  is  proven  by  LBwdin  in  PTX.  The  proof  is  simple 
and  uses  the  variational  principle: 

a1  (u1|a|u1>  = <u1JQAQ|u1>  = a1 

For  the  first  excited  state  consider  a normalized  variational  function, 

orthogonal  to  u^,  which  is  simply  a linear  combination  of  u^  and  u^. 

As  in  the  proof  of  Theorem  1.2,  it  follows  that  a^  - a^  and  we  obtain 

a^»  etc*  Q is  entirely  arbitrary.  In  particular  it  may  be  the 

projection  operator  associated  with  a linear  manifold  = (f.,f_, ...,f  ) 

'12  n 
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where  the  (f^)  are  linearly  independent  functions.  Then 

Q = |I>^-1«E|  = ^ lfkH^’1)ke(f^l  (1-53) 

kye 

where 

In  Appendix  A of  PTX,  LBwdin  discusses  these  projection  operators. 
Furthermore,  he  pointed  out  that  for  an  orthonormal  basis,  I = (<tv<t2, 
...,0n),  the  operator  QAQ  is  represented  by  a matrix  having  the 
elements  A^g  = (0k|A|0  ).  Theorem  1.8  then  asserts  that  the  eigen- 
values of  this  matrix  are  upper  bounds  to  the  eigenvalues  of  A,  and 
further,  that  these  bounds  improve  when  the  set  is  enlarged.  This 
is  a new  proof  of  the  Hylleraas  and  Undheim^^  and  McDonald^-^  result 
which  actually  goes  beyond  Theorem  1.8. because  it  states  that  the 
eigenvalues  of  approximation  k ( k-dimens ional  matrix)  are  separated  by 
those  of  approximation  k-  1.  Perhaps  it  should  be  pointed  out  that 
this  theorem  is  originally  due  to  Cauchy^^  and  Poincare^^ . A dis- 
cussion of  the  accuracy  of  the  variational  approximations  to  the 
excited  states  is  left  for  Chapter  II. 

Definition  1.4!  If  A is  a positive  definite  self-adjoint  operator,  we 

• i ! 

define  its  inner  projection  A'  by  A'  = A2QA2,  where  Q is  an  arbitrary 
projection  operator. 

Note  that  if  A is  a multiplicative  operator,  its  square  root  is 
simply  the  square  root  of  the  function  characterizing  the  operator.  If 
A is  a general  differential  operator,  the  square  root  may  be  defined  in 
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terms  of  the  spectral  resolution  A = £ ak^Uk^UR^’  where  ak  >0, 

JL  _ i, 

then  A2  = ak2|uk)(uk|  is  the  "positive  square  root." 

The  inner  projection  may  be  defined  for  a non-positive  definite 
operator,  but  only  in  its  domain  of  positivity  R.  For  instance,  let  0 
be  the  projection  operator  associated  with  the  subspace  R.  Then  0A0  > 0, 
and  A has  a positive  square  root.  The  important  property  of  inner 
projections  is: 


Theorem  1.9: 

If  A'  is  an  inner  projection  of  A,  then  0 < a'  < a.  As  both  Q 
and  1 - Q are  positive  and  bounded  by  1,  i.e., 


i s.  (tlQlt)  = <*|Q2k)  = (Q'Hq'iO^  o 
4 1 1 -Q 1 1 ) = l - (t|Q|t)^  0 

it  follows,  as  a consequence  of  Theorem  1.3,  that: 

0 < (t|A'|t)  = (t|A2QA2|\|r)k  (\|/|A|t)  . (1.54) 

Again  Q is  entirely  at  our  disposal.  One  may  define  inner  pro- 
jections on  linear  manifolds.  Lbwdin  introduced  the  inner  projection 
concept  in  PTX.  Specific  choices  of  the  projection  yield  formulas 
that  have  been  used  by  some  authors.  Consider  a linear  manifold  •£, 
with  metric  A = (C|  (£).  Then,  a generic  inner  projection  is: 

A'  = A^l  £>  A'1  <£|a2,  a = < ff|  ff>  . (1.55) 

By  the  transformation  A 2 ) = |g),  we  obtain 

A'  = A |g  > A ^’(g  | A , A = <g  | A | g ) , (1-56) 


[32] 


[33] 


the  Aronszajn  projection 


, of  extended  use 
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The  transformation  A2j-f)  " |h)  yields  the  Bazley  projection: 

A'  = |h)A>"1(h|,  A=  <h|  A"1  J h>  (I.57) 

This  form  was  successfully  used  by  Bazley  and  Fox  in  pioneering  work 

rui 

on  lower  bounds  to  the  eigenvalues  of  the  SchrBdinger  equationL  . 

The  application  of  the  operator  inequalities  introduced  in  this 
section  is  of  two  types: 

One  may  use  Theorem  1.2,  and  have  an  order  relation  between  the 
eigenvalues  of  A and  B,  A B.  This  is  the  "leitmotiv"  of  the  Weinstein- 
Aronsza jn-Bazley-Fox  methods.  This  implies  counting  the  lower  eigen- 
values and  assuring  that  no  eigenvalue  is  lost  in  the  construction  of 
the  "intermediate  Hamiltonians"  h,  with  E , in 

n — n 

order . 

The  alternative  is  to  use  the  partitioning  formula  in  Eq.  (l.jjl), 
and  compute  lower  bounds  to  (\|r|t|\jr).  If  t > 0 it  is  possible  to  con- 
struct an  inner  projection  t'  < t.  But  we  don't  actually  need  that 
t > 0;  in  fact,  as  we  shall  see,  it  is  possible  to  construct  lower 
bounds  to  through  inner  projections  of  V,  or  truncations  of  T or 
Tq,  of  combinations  of  these,  as  it  was  discussed  by  Lowdin  in  PTX, 
and  PTXI.  We  shall  present  some  further  comments  on  this  and  on  the 
connection  with  the  intermediate  problems.  Again,  an  inner  projection 
of  t can  be  constructed  within  its  domain  of  positivity,  the  condition 
t > 0 being  too  strong. 


CHAPTER  II 


ACCURACY  OF  EIGENVALUES  AND  EIGENFUNCTIONS 

There  has  been  considerable  work  on  the  accuracy  of  eigen- 
values and  eigenfunctions.  For  the  former,  upper  and  lower  bounds  are 
the  natural  estimates.  Local  deviations  .and  norm  deviations  of  the 
type  in  Eq.  (I.l6)  are  used  for  the  wave  functions.  In  this  section, 
we  shall  discuss  Temple-like  formulae  for  lower  bounds  and  Eckart-type 
criteria  for  accuracy  of  variational  wave  functions.  The  motivation 
is  that  Temple-like  criteria,  even  though  not  fully  ab ' initio,  have 
been  widely  used  and  are  extremely  accurate.  Approximate  variational 
functions  are  used  to  compute  all  kinds  of  properties,  for  ground  and 
excited  states,  and  these  depend  heavily  on  the  accuracy  of  the  func- 
tions. The  connection  between  Temple-like  formulae  and  partitioning 
formulae  will  be  established  and  the  use  of  lower  bounds  for  excited 
states  in  estimating  the  accuracy  of  excited  variational  wave  functions 
will  be  emphasized. 

2.1.  Temple  Formula 

It  was  obtained  by  Temple  ^ ^ , and  discussed  by  Kato^^ , 

Hylleraas^ ^ , and  LBwdin^^  . It  requires  the  calculation  of  expec- 

2 

tation  values  of  H , which  is  a computational  complication,  and  esti- 
mates of  neighbouring  eigenvalues,  not  necessarily  available.  This 

F "1  F ^ A 1 

last  requirement  is  not  needed  in  the  Weinstein1-3  J and  Crawford1,3 
formulae,  which  we  discuss  briefly:  let  E^  be  the  ground  state  of 
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Eq.  (l.l).  Then  the  variational  principle  states,  for  a normalized 
trial  function  \|r. 

Eo£  <*|H|t)  = (H)  . (2.1) 

If  (h)<’^(Eo+E1) , the  follow  ing  expression  is  positive: 

J = 4|[(H-(h>)2  - (Eo-(H))2j|t)  > 0 . (2.2) 

This  can  be  seen  by  expanding  ijr  in  terms  of  the  eigenfunctions  to  H, 
and  using  the  property  that  H is  self-adjoint: 


J - £ [(Ek-<H>)2  - (Eo-<H))2] 


\ >° 


k=0 


where  the  c,  ’ s are  the  coefficients  in  the  expansion  of  \Jr  in  terms  of 


the  * s . 
k 


or 


which  yields 


(h-(h) ) > (eq-(h))£ 


(AH)2  = (H2)  - (H)2  > (Eq-(h))2 


(2.3) 


E < <H)  - AH 


(2A) 


which  is  the  Weinstein- formula.  In  order  to  apply  it,  we  must 

know  that  (h)  is  "good  enough,"  in  the  sense  that(H)  ^(E^+E^). 

Temple's  formula  is  easily  derived  from  the  consideration  of 
the  operator  JX  = (H-E^)(H-E  ),  where  E^  and  Ej  are  neighbouring 
eigenvalues,  i.e.,  i = j+1.  and  H have  the  same  eigenfunctions,  and 
the  eigenvalues  are  c»^  = (E^-E^)(E^-Ej ) ^.0.  Then  Cl>~  0,  and  for  a 
normalized  trial  function  ljr: 
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and  let  j = i + 1.  This  yields  the  inequalities: 


(2-5) 


In  particular,  for  the  ground  state,  this  reduces  to: 


(2.6) 


the  well-known  Temple  result.  This  formula  has  been  used  by  several 


symmetry  is  needed,  or  at  least  a lower  bound  to  it.  This  limits  its 
applicability,  as  well  as  a fully  a priori  calculation  of  the  errors 


gave  a proof  of  Eq.  (2.6)  based  on  the  intermediate  problems.  In  PTIV 


in  Eq.  (1.12).  The  connection  has  not  been  established  explicitly 
before,  and  it  is  interesting  in  the  sense  that  it  allows  comparison 
of  the  partitioning  method  for  computing  lower  bounds  which  is  an 
a priori  method,  with  Eq.  (2.6)  which  is  an  accurate,  though  a non 
a priori  formula. 

In  order  to  derive  Eq.  (2.6)  from  Eq.  (1.12),  we  need  a bound 
for  T:  If  £ < E^,  Eq.  (1.3))  reduces  to: 


Lbwdin  suggested  that  (2.6)  should  follow  from  the  partitioning  formula 
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T > 


<£  - E, 


(2.7) 


and  theorem  1.3  implies  that 


HTH  > 


HPH 1 
f-  Ei 


Therefore  for  E < £ 4 E.,  we  have  a lower  bound  to  : 
o I 1 

£ = <tlH+HTH|^)>  <i|f|H+~^-U)  = 

1 e-El  1 


(2.8) 


which  means 


(ah)* 


E > <H)  - ±r 
° E,  - £ 


(2.9) 


E £ £ — E''  — E. 
o 11 


This  is  already  very  similar  to  Temple's  formula,  and  the  best 
possible  bound  it  can  give  is  for  £.  = Eq . We  do  not  know  E^,  but  any 
lower  bound  to  it  will  do,  because  the  inequalityv  is  preserved  if 

We  shall  present  some  lower  bounds  to  E^  derived  by 

LBwdin  in  PTX. 

Consider  Eq.  (I.29)  for  ^ E^,  which  reduces  to: 

| a |2  |a  |2  1 - | a |2 

o o 1 

where  a^ 

As  W(£.)  has  negative  slope  and  its  zeroes  are  the  eigenvalues 

of  PHP,  the  E.’s,  we  conclude  that  if  we  set  the  left  hand  side  of 
1 

Eq.  (2.10)  equal  to  zero  we  obtain  an  upper  bound  to  E^,  and  if  we  set 
the  right  hand  side  equal  to  zero  we  obtain  a lower  bound  to  E^.  Then 


(2.10) 


(d-EL)|ao|2  + (£.-Eo)(l  - |a f)  = 


0 


the  solution  of 
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is  a lower  bound  to  E^,  that  is 


E„(1'la0l2)  +Ell‘ol2<  E1 


(2.11) 


r 2 1 

At  this  state,  it  is  convenient  to  introduce  Eckart's  criterion1  J 
for  the  accuracy  of  an  approximate  wave  function.  In  its  operator  form 
it  states: 


(H-E  ) 

srrr  - 1 - Mv 

1 o 


The  proof  follows  from  considering  the  operator  H-Eq: 


(2.12) 


H 


•Eo  ■ z <Ek-Eo>°k  D <ErEo)0k  - (Ei-EoKi-0o)  a »• 

ktO  k^O 


Applying  Eq.  (2.12)  to  \|r  we  obtain  the  familiar  expression: 

(H ) — E ? 

nr*1-  KI2  • 

1 o 

Substitution  into  Eq.  (2.1l)  yields: 

E.  > E1  + (E  -E  ) ( 1 — | a |2)  > E.  + E - <H>  . 

1 1 ' o 1 ' ‘o'  1 o 


(2.13) 


(2.14) 


The  question  now  is  how  to  transform  Eq.  (2.9)  into  Eq.  (2.6).  We 
will  show  that  there  is  a value  of  £ , compatible  with  the  constraints 


E <.  £ 1 E,,  that  verifies 
o 1 

Ef  - £ = Ei  -(h) 

According  to  Eq.  (2.I5): 

£ - Ex  + (H)  - Ex  <.  Ex 

if  (h)<e1.  F urthermore,  Eq.  (2.14)  implies 

£ > E.  + E - (h)  + (h)  - E * E 
1 o I o 


(2.15) 
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Therefore  Eq.  (2.I5)  is  compatible  with  the  constraints  on  £, 

and  Eq.  (2.6)  follows  from  Eq.  (2.9).  If  we  could  have  better  lower 

bounds  to  E^  than  Eq.  (2.14),  Eq.  (2.9)  could  be  used  to  bypass  Temple's 

formula.  In  particular,  bounds  to  E^  independent  of  Eq  and  E^  will 

give  a modification  of  Temple's  formula  that  is  a priori  in  the  sense 

that  it  will  not  depend  on  E^,  and  we  will  present  those  in  Chapter  V. 

We  have  used  Eq.  (2.12),  which  is  a simple  form  of  Eckart's 

criterion.  This  shows  that  the  problem  of  obtaining  lower  bounds  to 

E.  and  E.  is  coupled  to  the  determination  of  the  coefficients  a.  in 
11  1 

y = Ea^T^.  But  this  is  an  important  problem  per  se:  given  a trial 
function  for  the  ich  excited  state,  how  close  is  it  to  y^?  Before  we 
may  proceed,  we  should  consider  this  problem: 

2.2.  Eckart-type  Formulae 

Our  goal  is  to  obtain  bounds  for  ||o^-y^||  = J"  | |^dv  where 

Hy,  = E,  y,  , and  is  the  k*"^  normalized  eigenvector  of  a secular 
k k k k 

problem  of  order  k.  Furthermore,  it.  will  be  interesting  to  have 

pointwise  bounds,  that  is,  bounds  for  l^(T)  ~ y^(^)l  . 

There  has  been  a considerable  amount  of  work  in  this  area.  To 

[41] 

consider  but  a few  of  the  authors,  we  should  mention  Shull  and  LBwdin 
Weinberger^-  , and  Lengyel^^  for  norm  bounds,  and  Diaz  and  Greenberg 
Kinoshita^^,  Synge and  Ridei^^  for  point  wise  bounds. 

The  bounds  presented  by  LBwdin  and  Shull  are  extensions  of 
Eq.  (2.13),  which  is  a bound  for  1 - |aj2  in  terms  of  (h)  - Eq.  Let 
us  denote  by  J the  kth  eigenvalue  of  the  secular  equation  of  order  M. 

K 
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Then 


<4lJHIV  ‘ We;  <W  ’ «k. 


(2.16) 


What  follows  depends  on  the  assumption  that  H is  a self-adjoint 
operator  bounded  from  below,  with  a partially  discrete  spectrum. 

Further,  the  results  are  valid  for  the  discrete  part  of  the  spectrum 
only,  below  the  continuum.  The  multiplicities  of  the  eigenvalues  must 
be  taken  into  account. 

Consider  J^-E^,  and  expand  in  terms  of  the  complete  set  ^ej» 


$ 


■ E 


. , a.  Y : 

k t—>  ke  e 
e=0 


jk  - \ - <\i“-Eki\>  - E iwi' 


ek 1 


e=0 


where 


ek 


- fr.iV5  E 


aek'  ' 


e=0 


(2.17) 


We  may  transform  Eq.  (2.17)  into: 

k-1 


Jk  ' Ek  ' E <VEk>laeJ2  + VI  <VEk>laek< 


e=0 


e=k+l 


k-1 


k-1 


E (Ee-Vlaekl2  + (Ek+rEk)(i  - E 


i%kl2  - lakkl2> 


e=0 


e=0 


k-1 

- - E (’WVkJ2  + <'VrV<1-l-kkl2>  • 

e=0 
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Therefore : 


k-1 


VEk  + 
1 - |akJ2-  


5-l  ^Ek+1-Ee^ ‘“ek1 


(2.18) 


Ek+1  - Ek 


LBwdin  and  Shull  point  out  that  the  problem  is  to  estimate  the 

la  , |2  in  terms  of  the  J.-E.,  We  should  note  that  the  first  term  in 
1 ek 1 ii 

the  sum  in  the  numerator  is  ( E,  , -E  ) |a  , |2.  The  nature  of  the  spectra 

' k+1  o ‘ ok ' 

of  quantum  mechanical  problems  is  such  that  the  differences  E^+^-Eq 
are  rather  large;  on  the  other  hand,  if  the  approximation  is  reasonable, 
|a  should  be  small  and  the  two  effects  cancel  out.  The  approxi- 
mations  involved  in  estimating  |ae^|  may  destroy  this  feature,  as  we 
will  see  now.  Following  LBwdin  and  Shull,  we  introduce  an  auxiliary 
function 

k oo 

®-  £ 


Y d 
e e 


(2.19) 


0=0 


e=0 


with  the  restriction: 


Clearly 


tv. 

^ Imo|2  = 1. 


a 


de  zLj  ^o^a* 

o=o 


d r = i 
e 1 


e=0 


(2.20) 


Let  us  choose  the  multipliers  such  that 


d = d,  = ...  = d,  , 
o 1 k-1 


= 0. 


(2.21) 
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Consider  now: 


<®|h|.®>  = £ 1 4J2ja  = X *d 


0=0 


e=0 


00  00 

aZ  |de|2Eea.Ek^|de|a>  Ek 
e=k  e=k 


which  may  be  rewritten: 
k-1 

Z <Je‘VI‘ie|2^  <Ek-Jk)l\|2  • . <2'22> 

e=0 


Assume  that  the  approximation  is  so  good  that  £ . . . £ ^ E^> 


hence  E,  -J  21  0 for  all  e < k and  Eq.  (2.22)  implies: 

1C  6 


2 J.  -E. 
x.  k k 

E.  -J 
k e 


(2.23) 


But  Eq.  (2.20)  and  Eq.  (2.2l)  yield  the  following  expression 


for  A . : 
ek 


hence, 


k-1 


, a 


ek 


7 a IS 
Z_>  ea  (i. 


a= 0 


(2.24) 


k-1 


k-1 


3ekl  ~ 


O 


hr 


0=0 


J,  -E  )2  i 

1 = (Jk'Ek)2  A (2-25) 


e=0 


E.  -J 
k e 


where 


k-1 

-z 

e=0 


<VJe>'Z 
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Replacing  this  in  Eq.  (2.18)  we  finally  obtain  the  LBwdin-Shull  formula: 


la.  . I2-  VEk 


kk  - Ek+rEk 


k-1 

M[> 

^ e=0 


k+l-Ee^ 


(2.26) 


,2  . 


If  A is  small,  which  is  the  case  when  J.  is  much  closer  to  E. 

l 1 

than  to  E^+^,  the  sum  of  the  terms  (E^+1-Eo)  is  multiplied  by  a small 


k+1  eJ 


number,  thus  compensating  the  fact  that  E^+^-Eq  is,  in  general,  rather 
large  in  atomic  and  molecular  spectra.  We  can  say  that  Eq.  (2.26) 
will  give  worse  bounds  with  increasing  k.  By  making  use  of  spectral 

[431 

representations,  Lengyel  J derived  a somewhat  simpler  formula,  which 
is  not  a consequence  of  Eq.  (2.2 6) 


1 - a 


2 . VEk  ( 


kk'  - Ek+rEk 


l1 


+ 


Ek+rE0 

Jk_Jk-i 


(2.27) 


The  same  considerations  apply  to  Eq.  (2.27):  The  fraction 

(E,  , , -E  ) (J,  -J,  ,)  increases  with  k and  can  become  very  large  even  with 

k k+1  o ' k k-1 

small  k,  limiting  the  usefulness  of  Eq.  (2.27)>  which  is  nevertheless 
slightly  better  than  Eq.  (2.26). 

I 

If  we  denote  by  E^  a lower  bound  to  E^,  Eq.  (2.27)  may  be 
replaced  by: 


0 J.  -E.  C J,  . . -E ' 
1 - la,  J2^  i 1 + k+1  ° 


*kk'  = Jk+l'Ek 


Jk-Jk-1 


(2.28) 


which  removes  the  dependence  on  the  exact  eigenvalues,  by  introducing 
the  lower  bounds  to  the  ground  state  and  the  state  under  consideration. 

The  last  equation  shows  that  the  problem  of  estimating  the 
accuracy  of  excited  state  trial  wave  functions  is  coupled  to  the  evalu- 
ation of  lower  bounds. 


4l 


Lengyel  also  presented  a bound  for  1 - |a.,  based  on  a 

K.K. 

successive  determination  of  bounds  for  1 - Ueg|2»  e=0,l,...,k,  but 
its  accuracy  diminishes  with  increasing  k.  On  the  other  hand,  Eq.  (2.27) 
and  Eq.  (2.28)  are  interesting  applications  of  the  spectral  theorem 
to  this  problem.  For  completeness,  we  present  a formula  due  to 

[42] 

Weinberger  who  made  a careful  estimate  of  the  terms  appearing  in 
Eq.  (2.18).  In  simplified  form  it  reads: 


1 


2 ^ Jk"Ek  f. 

‘ K+i-K  [ 


<Jk-Ek»Jk-rE;>) 

<vjk.i><Ek-E:>; 


(Jk-Ek)(Jk-rE4 

(Jk-Jk-i)(Ek-Ei) 


(2.29) 


Equation  (2.29)  has  the  advantage  of  having  the  factors  (j.  1 -E ' ) 
and  (E'-E')  explicitly  as  a ratio,  thus  there  are  no  unduly  large  terms 
spoiling  the  bounds,  and  it  reduces  to  Eckart's  formula  in  the  case  of 
the  ground  state. 

Formulae  of  the  type  in  Eqs.  (2.26)  - (2.29)  have  not  been  used 
very  extensively  because  we  lack  knowledge  of  accurate  lower  bounds  E|, 
but  should  be  kept  in  mind  for  future  reference,  because,  through 
Eq.  (1.27)  they  yield  bounds  for  the  Weinstein  function.  Furthermore, 

Eq.  (2.13)  as  well  as  Eqs.  (2.2 6)  - (2.29)  are  °f  the  type  |a^|^>  b _>  0. 
It  certainly  is  important  to  find,  if  at  all  possible  bounds  of  the  type 

1 I 2 [44] 

1 >c  2.  la^l  • Work  in  this  direction  was  done  by  Rayner  , which 
we  will  discuss  in  the  future,  in  connection  with  properties  other  than 
energy.  The  estimates  of  the  la^  I will  prove  useful  in  applying 
formulae  involving  the  Weinstein  function,  which  connects  the  eigenvalues 
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of  H with  the  eigenvalues  of  PHP,  hence  whenever  reduced  resolvents  for 
arbitrary  reference  functions  are  employed. 

2.3.  A Numerical  Appraisal  of  Eckart-type  Formulae 


Before  going  into  quantum  mechanical  applications  of  Eq.  (2.13) 

and  Eq.  (2.26)  to  Eq.  (2.29)>  it  is  necessary  to  realize  that  they  may 

give  worse  bounds  than  expected.  Both  the  Eckart  formula  in  Eq.  (2.13) 

and  the  Temple  formula  give  sharp  bounds.  By  this  we  mean  that  they 

become  equalities  instead  of  inequalities  under  special  circumstances, 

namely,  that  the  trial  function  considered  is  a linear  combination  of 

the  ground  and  the  first  excited  states  of  H . But  this  event  is  very 

o 

unlikely  in  general.  It  is  important  then,  to  compare  the  different 
formulae  in  a case  where  we  know  the  exact  solutions.  For  this  purpose 
it  is  enough  to  consider  a simple  3*3  matrix: 


(2.30) 


The  eigenvalues  can  be  obtained  by  solving  the  characteristic 
equation  or,  most  easily,  by  partitioning.  They  are: 

= 0.579166 

a2  = 3-370846  (2.31) 

a = 23.  049988 


The  first  two  eigenvectors  are: 
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ui  = 


I o. 921570 \ 
-0.387829 
^ 0.01727 6j 

' 0 . 338207 

0.520378 
^-0.046889 

We  shall  consider  a projection  PAP  with 


(2.32) 


,>.is  gives  for  A: 


I 1 0 0 


P = 0 1 0 


A = 


0 0 0 


3 j 


(2.33) 


(2.34) 


with  eigenvalues 


a.  = 2 - = 0 . 585786 


a2  = 2 +Y2  = 3.414214 


(2.35) 


and  eigenvectors: 


(2.36) 
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We  may  test  Temple's  formula.  For  this  we  need  A2.  A lower 


bound  to  a^  is: 


(G1|a2|u1>-<G1|a|u1)2 


a1>  ax 


a2  - al 


(2.37) 


as  a trivial  application  of  Eq.  (2.6).  This  gives 


ax  > 0.58786  - 0.052583  = 0.533203 


which  is  not  a very  good  lower  bound  and  in  fact  is  suggesting  that  we 
should  not  expect  very  good  results  from  the  Eckart  criterion  in  this 
example . 

We  can  now  compute  (u-Ji^)  and  (u^Ug)  exactly: 

ail  = <U1  1^1 > = 0-999835;  a22  = (u2  | u2 ) = O.998878 


which  means  then: 


1 - a' 


11 


0.000331,  1 - a22  = 0.00224 


Eckart's  criterion,  Eq.  (2.I3),  applied  to  Eq.  (2.30)  gives: 

a,  -a. 


1 . a2  A 

11  a2~ai 


= 0.00237  . 


The  bounds  for  1 - a^  are: 


(2.38) 


(2.39) 


1 ~ a22—  O.OI99  from  Eq.  (2.26) 

1 " a22  ~ °-0197  from  Eq.  (2.27)  (2.40) 

1 - a22  - 0.00228  from  Eq.  (2.29). 

We  see  that  Eckart's  criterion  gives  considerably  more  error 
than  Temple's  bound.  We  said  before  that  both  formula  are  sharp,  and 
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this  follows  from  considering  \|r  = a ¥ + a.^.  : la  |2  + |a  |2  = 1.  Then 

o o 11  o 1 

<H>  - |ao|2Eo  + UJ^  « Ex  + 
or 


which  is  Eckart's  Eq.  (2.13)  when  the  equality  holds.  Similarly  if  we 
consider  (H-Eq)(H-E^),  its  expectation  value  over  \[r , is  zero,  i.e.: 

(AH)2  + ((h)-Eq)((h)-E1)  = 0 

from  which  Temple  formula  follows  as  an  equality. 

The  bounds  for  1 - a2^,  from  Eqs.  (2.26)  and  (2.27)  have, 
roughly,  the  same  error  as  Temple's  formula,  but  Weinberger's  bound  is 
very  accurate.  If  this  is  true  in  practical  application,  it  remains  to 
be  seen  but  at  this  time  the  existence  of  these  formula  suggests 
further  work  in  connection  with  W(£). 


CHAPTER  III 


PARTITIONING  TECHNIQUE  AND  PERTURBATION  THEORY 

In  PTIV  and  PTX,  LBwdin  discussed  the  various  connections 
between  partitioning  technique  and  the  Br illouin-Wigner  (B-W)  and 
Raileigh-SchrBdinger  (R-S)  perturbation  expansions.  They  stem  out 
of  Eq.  (I.50)  by  using  the  relation  in  Eq.  (I.52)  for  T.  LBwdin. 
discussed  the  derivation  of  the  various  expansions  and  some  convergence 
criteria. 

In  this  chapter  we  intend  to  show  the  connection  with  another 

r r 1 

expansion,  the  Feenberg  perturbation  series  L . We  present  a deri- 
vation of  the  Feenberg  series  by  using  successive  partitioning  tech- 
nique, and  discuss  some  of  its  properties  as  well  as  methods  of  solution 

^.1.  Br illouin-Wigner  Perturbation  Theory 

If  in  Eq.  (1.1*8)  we  let  \|/  = where  ij/^  is  an  eigenfunction 

of  Hq,  and  apply  Eq.  (I.52)  we  get  a BW  type  series,  namely, 

00 

'Fr  = ; (T  VW  = \J r + T V*  + T VT  V\|r  + ... 
c l__i  ' o o o 0 0 000 

k=0 

(3-1) 

00 

- *0  Z (ToV>\ 

k=0 


1*6 
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4 ■ Eo  + Z ■ Eo  + <*olVl*0> 


k=0 


+ <*0|VT0V|*„>  + ■••  * Eo  + <*„IV  Z <ToV)k|»o>  • 


(3.2) 


k=0 

r4qi  r c=o  i 

This  series,  due  to  Brillouin  and  WignerL^  J has  been 
extensively  discussed  in  the  literature,  where  the  interest  has  been 
focused  in  the  points  £ “£.=*  E. 

If  we  restrict  ourselves  to  positive  definite  perturbations, 

V > 0,  and  consider  values  of  £.  < E°,  from  Eq.  (I.30)  it  follows  that 


T 


o £-eJ  d-EJ 


f:  0 


(3-3) 


-1 


hence  - T 0.  From  theorem  1.6,  V > 0 implies  V >0.  Applying 
these  results  to  Eq.  (1.46): 

t'1  = v"1  - T > V"1  > 0. 
o 


By  theorem  1.7  again. 


0 ^ t ^ v 


(3*4) 


Therefore,  the  value  of  is  bracketed  by: 

E-E  + (\|r  1 1 1 \|r  ) = <5,  — E + (llf  I v|  \|/-  ) . 
o o To‘  ,Yo  1 o yo‘  ,Yo 


(3-5) 


On  the  other  hand,  the  series  expansion  may  or  may  not  converge. 
Examination  of  the  terms  appearing  in  Eq.  .(3*2)  shows  that  the  ones  of 
even  order  are  positive,  and  that  those  of  odd  order  are  negative. 
Therefore  the  series  is  alternating,  and  it  will  converge  if  and  only 
if  the  magnitude  of  the  n-th  term  goes  to  zero  when  n goes  to  infinity. 


k8 


We  can  write  this  condition  on  the  terms  of  Eq.  (3-2)  as  a 
condition  on  the  operator  t,  expressed  in  symmetric  form: 


= V2^  (VZIV^M 


k=0 


which  will  converge  if  and  only  if 


or  simply 


JL  _l 

-1  < V2T  V2  < 1 


V"1  > -T 

r r 


(3-6) 


( 3 -T) 


Using  the  truncation  Eq.  (I.30)  as  applied  to  T , it  follows: 


-T0  < -To(p) 


Hence  it  is  sufficient  to  impose 


V"1  > -To(p) 


(3.8) 


(3.9) 


to  assure  Eq.  (3*7).  Therefore  a sufficient  condition  for  the  conver- 
gence of  the  Br illouin-Wigner  series  is  that 


1 1 

1 > -V2T  (p)V2 


j j 

which,  for  p = 0,  To(p)  = — q j simplifies  to 


(3-10) 


£ - En 


1 1 


1 1 


V2PV2  = V - V20V2  < E Y-l 


(3-11) 


which  is  certainly  satisfied  if 


V < E^  - £ 


(3-12) 


The  preceding  analysis  is  very  rough;  we  may  have  some  additional 
insight  if  we  consider  the  following  theorem. 
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v 


Theorem  3 . 1 . 

For  our  operator  B,  a necessary  and  sufficient  condition  for 

niS  B = 0 is  that  there  exist  a positive  definite  Hermitian  operator  A, 
for  which  A - BAB  is  positive  definite^^  . 

The  well-known  result  is  that  Bn  = Q if  and  only  if.  the 

n-*oo  J ’ 

eigenvalues  of  B are  all  of  modulus  less  than  unity.  We  present  here 
the  proof  of  the  sufficiency  in  theorem  3.1,  as  it  is  the  result  we 
are  going  to  use. 

Let  A be  any  positive  definite,  Hermitian  matrix  for  which 


A - BAB  > 0 . 


(3.13) 


+ 


Since  A > 0,  A can  be  written:  A = DD  and  be  replaced  in 

Eq-  (3-13): 


+ + + 

DD  - BDD  B > 0 . 


+ \-l 


By  theorem  1.3,  setting  /l  = (D  ) 

1 - KK+>  0 (3-14) 

where  1 is  the  identity  operator  and  K = D ^"BD;  this  implies  that  the 
modulus  of  every  eigenvalue  of  K is  smaller  than  unity.  But  K is  a 
similarity  transformation  of  B.  Hence  all  the  eigenvalues  of  B are,  in 
modulus,  smaller  than  unity,  thus  Bn-»  0 when  n-*-°°. 

Consider  now  the  series: 


(A-B)"1  = A"1^(BA"1)n;  A+  = A > 0;  B+  = B < 0. 


(3.15) 


n=0 


The  power  series  converges  iff  BA  1 = 0.  According  to  the  theorem. 


n-*°° 
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the  series  will  converge  iff  we  can  find  an  operator  G,  Hermitian  and 
positive  definite  such  that 


c > ba"1ca“1b 


(3.16) 


Let's  set  G - A > 0.  Then  the  condition  for  convergence  of 

Eq.  (3.15)  is 

A > BA_1B  . (3.17) 

In  fact  we  have  C entirely  at  our  disposal.  We  may  make  other 


choices,  for  instance  C = A which  yields 


2 2 
Ac  > D 


(3.18) 


as  a sufficient  condition  for  convergence  of  Eq.  (3.I5). 

Notice  that  Eq.  (3.I8)  is  not  a trivial  consequence  of  Eq.  (3.I7). 

Only  if  A and  B commute,  Eq.  (3.18)  follows  from  Eq.  (3.I7). 

The  relevance  of  these  relations  in  BW  theory  follows  from 

identifying  A with  V and  B with  T . Then  sufficient  conditions 

o 

for  convergence  are: 


v_1>  TqVTo  or  1 > (AvM  (3.I9) 

which  is  the  same  condition  as  Eq.  (3.6). 

From  Eq.  ( 3. 18) , we  obtain: 

V"2>  T02  or  1>VTo2v  . (3.20) 

Now  T and  T (p)  commute,  then  from  T > T (p)  it  follows 
o o 00 

2 2 

Tq  4 (p),  so  it  is  sufficient  to  impose 

1 > vto2(p)v 


(3.21) 


which  for  p - 0 simplifies  to 
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(E°-£)2>VPV  (3.22) 

which  should  be  compared  with  Eq.  (3.H).  It  is  a better  bound,  though 
its  applicability  doesn't  go  beyond  the  relation  in  Eq.  (3-12)  because 
the  best  bound  for  VPV  is  VPV  < V2,  hence  Eq.  (3.22)  is  certainly 
fulfilled,  if  (E^  - £)  > V which  is  obtained  by  squaring  relation 

Eq.  (3.12). 

3.2.  Upper  and  Lower  Bounds  in  BW  Perturbation  Theory 

This  was  discussed  by  LBwdin  in  PTX.  The  starting  point  is  the 
symmetric  expansion  of  t: 

t = V + VTo  + VTotToV  (3.23) 

which  leads  to 

2n-l 

c = V Z (ToV)k  + (VTD)nt(ToV)n  . (3.24) 

k=0 

If  V > 0,  £,  < E°,  we  have  seen  that  0 < t < V.  But  this 

means  then: 

(VTo)nt(ToV)n<  (VTo)nV(ToV)n  (3.25) 

and  therefore 

2n-l  2n 

V Z (T0v)k^ToV)k<:  C < V E (3-26) 

k=0  k=0 

which  is  a bracketing  of  the  operator  t.  In  addition,  it  is  clear  that 
partial  sums  to  even  order  in  the  BW  expansion  yield  lower  bounds  to 


52 


the  reaction  operator,  and  odd  orders  provide  upper  bounds.  If  we 
define  the  functions  X^  and  the  matrix  elements  by: 

\ ' <T„v>kV  ek+i ' <*0MV  ■ (3.27) 

The  series  in  Eq.  (3*2)  may  be  written: 

00 

£1  - Eo  + E fk  (3.28) 

k=l 

where  the  functions  X^  are  defined  through  inhomogeneous  equations. 

If  EC  £ C E°,  Eq.  (3.26)  leads  to: 

2n 

^1>E0  +Z  «k  • (3.29) 

k=l 


Therefore,  the  partial  sums  to  even  orders  are  lower  bounds  to 

v V 

the  ground  state  energy,  provided  that  V > 0,  and  <£.  4 E°.  Several 
r g] 

authors1  - have  derived  variational  methods  to  obtain  upper  bounds  to 

2n 


the  partial  sums 


Unfortunately,  it  would  be  more  interesting 


to  obtain  lower  bounds  instead,  because  in  such  a case  the  series 
would  still  provide  a rigorous  lower  bound  to  the  ground  state  energy. 
The  need  for  variational  methods  comes  from  the  fact  that  Eq.  (3>27) 
cannot  in  most  cases  be  solved  exactly  for  X^.  We  shall  discuss  in 
the  next  chapter  the  variational  methods  employed  in  solving  equations 
of  the  type  in  Eq.  (3.27)  in  RS  and  BW  theory. 

Going  back  ;to  Eq.  (3-24),  LBwdin  suggested  that  the  remainder 
(VTo)nt(TQV)n  may  be  replaced  by  an  inner  projection  on  t,  t^  = t20t2, 
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in  which  case : 

2n-l  2n 

V (T0v)k  + (VTo)nt'(VTQ)n  < T < V Y (T0V)k  • (3-30) 

k=0  k=0 

He  showed  that  a suitable  choice  of  the  inner  projection  leads  to  an 
estimate  of  the  remainder  in  terms  of  quantitites  already  evaluated. 

The  following  formula  is  obtained  in  PTX. 


2n 


£ l > E°  + ^ ^ + ®:+  A 


(3-31) 


k=l 


where  A^e  - 6 k+e+1  " ^k+e+2’  ( > e=0, 1,  2, . . . ,m-l)  and  is  a row  vector 

of  order  m:  « - (6„+1.  £„  J • 

For  a given  n,  m = 1,2,...  gives  an  increasing  series  of  lower 

bounds,  and  for  n = m they  involve  only  quantities  6,  up  to  f 

k 2n+l 

For  example,  if  n = m = 1 


> E°  + + £ + 


1 ° i 2 e1-e2 


(3.32) 


The  fact  remains  that  the  quantities  6^  must  be  evaluated. 


Equation  (3.27)  means: 


X = T VX  X = \|r 

k o k-1  o vo 


(3-33) 


Hence,  by  Eq.  ( 1.7) ; 


or  s imp ly 


P(£  - H )X,  = PVX, 

^ o'  k k-1 


(£  -H  )x  = VX  - € ^ 

' o k k-1  kYo 


(3.3*0 
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and  these  inhomogeneous  equations  cannot  be  solved  exactly.  At  this 
stage,  it  may  be  convenient  to  reconsider  Eq.  (1.46)  and  introduce  the 


resolution  of  T . 

o 


3.3.  Truncations  of  in  BW  Theory 

In  what  follows,  will  refer  to  an  unperturbed  eigenfunction 

to  H , not  necessarily  the  ground  state.  We  can  give  bounds  for  T , 
o o 

as  indicated  by  Eq.  (1.30).  Notice  that  the  0^.’  s are  identical  with 
the  0,  ’s  just  because  \|r  is  a single  eigenfunction  to  H . Then  we 


have,  for  t < E 


p+1 


s I 


P ' 0° 


k=0 


£-e{ 


> T 


‘I 


P ' 0° 


k=0 


£-E 


T (p) 


P ' 0° 


k=0 


(3.35) 


where  the  prime  in  the  summations  indicates  exclusion  of  the  term 
involving  0g  = |i|r°)  (tyg  J where  \J/°  is  equal  to  i|r  . 

Instead  of  solving  Eq.  (3.2),  we  may  be  tempted  to  consider 


instead : 


- E° + <t0|[v'1-Typ)]‘1|*0> 


■ Eo  + E b0lv[i;(p)v:klt0> 


(3.36) 


k=0 


Eo  + ^olVI^o)  + <^0lVT0(p)Vl^0>  + 


/ 
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The  individual  terms  are  easily  evaluated,  i.e.,  for 


<*0|VT>)V|*0>  * £ FI? 

k=b  t_Ek 


V V 

ok  ko  _ r 1 

C2 


> e 


and  similarly  for  the  higher  orders.  Again,  we  see  that  we  get  upper 
bounds  for  the  second  order,  when  it  would  be  desirable  to  have  a lower 
bound  instead.  Further,  we  do  not  have  bounds  for  the  higher  orders. 

The  interest  of  considering  Eq.  (3-36)  at  all  is  that  we  can  evaluate 
to  all  orders  by  just  solving  an  inhomogeneous  equation.  Consider 
Eq.  (3-36)  written  in  the  form: 


*eo  + 

Define  a function  X by  the  inhomogeneous  equation: 

[v~1-t'(p)]~ V,  = X , 


which  may  be  written: 


-1  ' It°  )(t°lx) 

*o=V 


k=0  £ “ Ek 


which  is  equivalent  to  the  system  of  equations: 


(3.37) 


(3.38) 


V 

eo 


•.-E 


Vl 

£ -E? 


(3.39) 


where  V , 
ek 


ck  - 4°|x). 


We  are  interested  in  C which  is 

o 


given  by 
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c 

o 


= v + y 

oo 


VQkCk 

£-e° 


(3.40) 


For  a given  value  of E,  one  solves  the  system  of  linear  equations  in 
Eq.  (3-39)>  and  replaces  the  values  of  the  C^’s  in  Eq.  (3*40),  thus 
solving  Eq.  (3.37)  exactly.  We  should  realize  that  we  are  not  neces- 
sarily bracketing  an  eigenvalue  of  H.  In  fact,  as 


t_1  = V’1_To  ^V_1-T’(P)  = t-1(p) 


if  t 1(p)  > 0, 


by  theorem  1.7, 


(3-41) 


t'(p)  > t (3.42) 

and  £ | is  an  upper  bound  to  Therefore  if  £ is  an  upper  bound 

to  E , it  does  not  follow  that  S'  is  a lower  bound.  In  this  connection, 

it  is  interesting  to  study  the  point  £.'=£.=  e'.  In  order  to  do  so, 

1 o 

let's  rewrite  Eq.  (3*37)  in  a way-' that  makes  explicit  the  dependence 
on  specific  states.  For  the  state  g: 


£|(g)  = E°  + (^rtv-^p)]-1!^) 

where 


*:(P)  = 


6-H 


P = 1 


WI 


Then,  we  may  rewrite  Eq.  ( 3 - 39 ) : 


V = C 
eg  eg 


r 

I 


V . C. 
ek  kg 

e-Eu 


(3.43) 


(3-44) 
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where  C - (i|r°|X  };  <f'(g)  = E°  + C . 

eg  ye‘  g lvs'  g gg 

^ = £-^(g)  we  obtain  finally: 


If  we  impose  the  condition 


The  system  of  equations  in  Eq.  (3.45)  has  a solution  if: 


V 


ek 


|A|  --  Det[£l(lK  ' 8ekl ' °' 

By  multiplying  the  k1"*1  column  of  A by  £ '(g)-E°  we  obtain: 


7T  (^(gK)'1  Det[V ek-(£;(g)-E“)5ek]-  0. 
k=0 


(3.45) 


(3.46) 


But  (i|r°|H|\|r°)  = E°5gk  + Vek,  hence  Eq.  (3-46)  is  equivalent  to: 

Det[Hek-^(g)]  = 0 (3.47) 

which  means  that  the  points  € = £.j(g)  are  nothing  but  the  roots  of  the 

secular  equation  constructed  with  the  basis  i|/-°, \|r°,  . . , \|/°,  and  by  the 

variational  theorem  they  are  upper  bounds  to  the  eigenvalues  Ek<  This 

shows  that  the  truncation  on  the  left  hand  side  of  Eq.  (3.35)  does  not 

lead  to  very  useful  results  for  lower  bounds,  but  it  provides  the 

explicit  solution,  in  terms  of  partitioning  of  the  secular  equation  in 

a basis  of  eigenfunctions  to  Hq. 

This  should  not  be  at  all  surprising,  because  Eq.  (3*47)  is  the 

secular  equation  arising  from  consideration  of  the  truncated  Hamiltonian 

H'(p)  = H'(p)  + V,  where  H'(;p)  = V E^O.0  > H . As  H'(p)  > H,  the  result 

k=0 
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is  k-0  compatible  with  the  ordering  theorem.  The  truncation  T^(p)  or 
HQ(p))  restricting  the  operators  to  one  subspace,  has  as  a consequence 
the  loss  of  bracketing  of  the  exact  eigenvalues,  unless  we  feed  in  a i 
which  is  a lower  bound  itself,  but  this  is  not  a useful  approach.  On 
the  other  hand,  the  study  of  Eq.  (3.44)  provides  an  interesting  test  case 
for  the  partitioning  method  as  applied  to  the  solution  of  a secular 
equation.  This  was  done  by  Meath,  Sando,  Goscinski  and  Hirschf elder 
and  the  preceding  discussion  places  their  work  in  the  general  framework 
of  truncations  of  Tq.  In  the  notation  used  in  this  paper,  what  they 
have  done  is  to  consider  Eq.  (3.44),  and  solve  it  by  an  iterative  pro- 
cedure : 


C = V + 
gg  gg 


E E° 


v i ci 

-gk-hg- . • 


(3-48) 


— ~ +C  -E. 
k^g  g gg  4 


By  choosing  a starting  value  for  C , the  system  of  Eq.  (3-45) 

oo 

can  be  solved  for  , which  when  replaced  in  the  r.h.s.  of  Eq.  (3.48) 

yields  a new  value  for  C and  the  procedure  is  repeated  till  conver- 

gg  K 

gence.  In  Appendix  I we  present  some  numerical  results  which  illustrate 
the  convergence  in  terms  of  the  derivative  of  an  appropriate  bracketing 
function. 

A final  comment  should  be  made,  and  it  is  in  regard  to  Eq.  (3.42) 
We  stated  that  if  t \p)  > 0,  Eq.  (3-42)  followed.  The  way  to  insure 
that  t ^(p)  > 0,  was  to  require  V > 0 and  T^(p)  < 0,  hence  we  should 
impose  the  condition  d<E°,  yet  by  virtue  of  the  connection  with  the 
variational  principle,  we  are  able  to  arrive  at  the  conclusion 
£^(g)  — Eo,  in  order.  This  points  to  the  fact  that  the  condition 
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t (p)  > 0 is  too  strong,  and  we  were  able  to  bypass  it  by  solving  the 

equations  and  relating  them  to  a variational  solution.  We  will  return 

to  this  point  in  Chapter  V,  where  we  discuss  the  more  interesting 

truncation  T > T (p) 
o o' 


3-4.  Feenberg  Perturbation  Theory 

As  we  have  seen  in  section  3-2,  the  convergence  properties  of 

Br illouin-Wigner  perturbation  theory  need  to  be  studied  in  considerable 

f53l 

more  detail.  Morse  and  Feshbach  carry  out  a thorough  analysis. 

For  the  present  purposes  we  would  like  to  point  out  that  the  conver- 
gence peculiarities  derive  from  the  power  series  expansion.  A simple 
way  to  look  at  this  problem  is  to  consider  Eq.  (3-37)-  It's  solution 

is  given  by  Eq.  (3*39)  an<3  Eq-  (3-40).  We  may  try  instead  a BW  type 

expansion  for  Eq.  (3.37).  We  obtain  the  infinite  series  Eq.  (3.36). 

If  we  consider  now  the  individual  orders,  we  have  for  second  order. 

(3-k9) 

k=l  k 

which  blows  up  for  £ = E°,  k = l,2>...,p.  On  the  other  hand,  the  exact 

solution  of  Eq.  (3-39)  is  well  behaved  at  these  points,  as  it  can  be 

seen  by  rearranging  it  in  the  following  way: 


V (£  - E°)(  £.-  E°)  C - V VokCk 
eo  m' ' nr  e 

k^O  (£  " Ek) 

k^m 


V C , 
om  m * 


Let  £ - E°;  then  V ,C  - 0,  and  if  V ^ 0,  C =0  excluding  the 
m om  m om  m 6 
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pathological  term  in  Eq.  (3-39) • This  shows  that  the  analytical 
behavior  of  £|(p)  is  not  reproduced  accurately  by  the  different  terms 
of  the  power  series  expansion,  and  this  curtails  the  convergence  prop- 
erties of  the  expansion.  These  remarks  apply  just  as  well  to  the  BW 
series  with  Tq  instead  of  T^.  Let's  see  what  happens  with  more  detail. 
The  third  order  energy  in  the  BW  series  is,  for  the  state  g: 

€3(S)  * (t°|VToVTo|*°)  (3.5O) 

and  in  terms  of  the  spectral  resolution: 


<h(s) 


-E 


k^g 


V V V 
gk  ke  eg 

(£-  E°)(£-E°) 


(3.51) 


where  there  is  no  restriction  k^e.  But  this  means  that  the  terms 
^gk^kg^kk^^ ~^k^  are  aH°wed*  And  this  type  of  term  is  unwanted. 

If  we  consider  the  characteristic  equation  of  H in  a basis  of  eigen- 
functions to  Hq,  for  a finite  or  an  infinite  basis,  there  should  not 
appear  any  coefficient  which  involves  products  of  matrix  elements  of  the 
same  row  or  column  of  the  secular  determinant: 
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■ V 
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El+Vll' 
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21 


12 


E“+V  - £ 
2 22  c< 


= 0.  (3.52) 


In  Eq.  (3.52),  a term  involving  Vq^v^qV11  can  not  aPPear-  Yet>  as 
Eq.  (3.51)  indicates,  this  term  will  appear  in  the  solution  of 

£x(g) = Eg  + vgg  + ^g)  + S(s)- 


6 1 


[ 5 ] 

As  Feenberg  ' pointed  out,  these  unwanted  terms  may  be  removed 
by  what  he  called  the  rearrangement  formalism.  A typical  term  may  be 
indicated  by  a diagram  - see  figure 

The  diagram  represents  the  product 


V V V V V V 
lr  rt  ts  s2  22  51 


and  the  "intermediate  states"  2 and  s occur  twice.  Feenberg  showed 

that  it  is  possible  to  eliminate  the  repeated  terms  if  one  modifies 

the  denominators  in  the  BW  series.  There  are  several  ways  to  go  about 

[R4I  r^i 

this.  Feenberg  s rearrangement  formalism  y , Bisshopp's  asymptotic 
approximation  to  a Fredholm  type  solution  and  Morse  and  Feshbach's 
treatment,  based  on  the  expansion  of  the  characteristic  equation. 

[19] 

KumarL  7 introduces  a rearrangement  operator  and  a "Feenberg  potential" 
in  order  to  derive  the  equations  which  seems  a bit  too  dangerous, 
considering  that  the  equations  are  nothing  but  a rearrangement.  In 
fact,  we  intend  to  show  that  they  can  be  derived  by  the  success ive 
partitioning  technique  discussed  by  NordlingL'?°-' . The  significance 
of  this  is  important.  The  partitioning  technique,  PT,  pursues  the 
solution  of  the  secular  equation  via  a successive  partitioning  technique, 
SPT.  A complete  solution  of  the  system  of  equations  is  not  performed, 
only  a first  triangular ization  of.  the  solution  being  needed.  We  shall 
show  this  in  the  next  section  where  we  will  prove  that  SPT  applied  to 
Eq.  (3-44)  yields  the  pt4  approximant  to  Feenberg1 s perturbation  scheme, 
to  which  it  is  entirely,  equivalent. 

We  shall  review  several  derivation  of  Feenberg1 s formalism 
and  then  present  the  successive  partitioning  derivation  which  illuminates 
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some  properties  of  the  method,  in  particular  the  connection  with  the 

continued  fraction  representation  of  the  solutions. 

[I9I 

Kumar  w defines  Pp,  the  Feenberg  operator,  to  be  such,  that  in 
the  expansion  of  (\j/^  j P^tp  | ) , the  state  \|r  is  not  equal  to  any  inter- 

mediate state,  and  where  tp  is  a modified  reaction  operator,,  involving 

the  Feenberg  potential,  which  he  defines  as  = Y*  (t„),  , 0°  So,  the 

F 4_j  F kk  k 

operator  being  considered  is,  by  making  the  formal  translation  into 
the  language  employed  in  this  dissertation 


t = V + V 7: t_ 

F 6 - H -V  F 

o F 


(3-53) 


This  definition,  which  should  be  compared  with  Eq.  ( 1.44a) 
implies  the  following  expansion: 


E = v + v£-h  _v; 


p p p 

F F F 

V + V v " „ „ V 7 ■„  „ V + ...  . 


o F 


£-H0-Vp  £-Hrt-V. 


o F 


If  we  use  the  spectral  resolution  of  the  operator 


and  noticing  that  the  definition  of  V_,  entails  that 

F 


£-Vvf 


H + V 
o F 


I K' 


we  obtain 


V0°V 

t = V + > 

F L 1-.  E°(.> 
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The  effect  of  is  to  impose  the  injunctions  on  the  indices  of 
summation.  If  we  use  the  following  notations: 


U2) 


pqr...,n  = En  + (tF}nn’ 

P^q,r,  . . 
q^r, s . . . 


(3-55) 


Eq.  (35*0  may  be  rewritten: 


tF  = V + 


E 


vo,°v 

k 


+ 
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vo.°vo°v 

k 


ku  ^O  [£-(k2)  ,][£-(x.2)  ] 

k^°  o5k  o,k  okj 


(3.56) 


vo°vo°vo°v 

ki  m 


i^o,k 

m^o,  k,  Jj 

Feenberg's  perturbation  series  follows  from  Eq.  (3.5 6)  by- 
setting  £ = <£,  = E in: 

fi,  - E°  + (O  - E°  + . (3-57) 

1 o F 00  o Yo*  F'Yo 

The  burden  of  the  proof  now  is  to  establish  that  (t  } = (t) 

F 00  00 

One  should  realise  that  this  is  very  formal,  and  the  Feenberg  operator 
may  not  be  definable  for  all  states  appearing  in  the  expansion  Eq.  (3.57). 
So  instead  of  reproducing  Kumar's  proof,  which  is  subject  to  criticism, 
we  shall  point  out  that  Eq.  (3-57)  may  Ee  rewritten  in  the  form: 


£1- 


E + V + 
o 00 


VQkVko 

,,  £-(^2)0k 

k^o  °> K 
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(3.58) 
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where 


(k  ) , - E°  + V.  . + 

' o,k  k kk 
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£?o,k  e'(X,2)ok^ 


(3.59a) 
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VkijVfeVmk 
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and  similarly; 


^2)ok,e  - El  + vu  + 


lm Vmjj 


m/o, k,£ 


£-(?) 


oki,m 


m^o,k,i 
n^o, k,$,m 


V,  V v . 
H m ma  ajj 


lkk,n] 


(3.59b) 


and  so  on. 

It  is  quite  clear  that  the  expansion  in  Eq.  (3.58)  with  the 
denominators  in  Eqs.  (3,59a)  and  (3. 59b)  seems  to  fulfill  the  require- 
ment that  no  coefficient  has  products  of  terms  V. . with  repeated 
indices,  i.e.,  terms  that  are  forbidden  in  the  expansion  of  a secular 
equation.  In  fact,  as  Feshbach^^  has  shown,  the  Eqs.  (3.58)  and 
(3.59)  represent  the  exact  solution,  by  elimination,  of  the  system  of 
linear  equations  whose  solution  is  determined  by  the  secular  Eq.  (3.52). 
In  other  words,  given  an  infinite  secular  equation,  or  a finite  one, the 
eigenvalues  are  given  by  Eqs.  (3.58)  and  (3-59).  If  the  basis  is  finite 
the  injunctions  on  the  indices  in  Eq.  (3.5k)  imply  that  there  will  be 
a finite  number  of  terms  in  such  an  expansion,  in  fact  it  is  easy  to 


see  that  for  an  N dimensional  bas 


6 


PZ 


is,  the  zth  order  in  the  expansion  of 

Eq-  (3*54)  has  (3(N,  z)  = (N-l)(N-2)  . . . (N-z+l)  terms  in  the  sum,  hence  in 

the  expansion  Eq.  (3*54),  the  orders  N+l  and  higher  are  vanishing.  The 

fact  that  Eq.  (3*58)  can  he  obtained  from  the  characteristic  equation 

constitutes  a proof  of  Feenberg's  expansion  and  of  the  formal  validity 

of  the  expansion  obtained  by  Kumar.  If  we  apply  this  elimination  to 

Eq*  (3*37).  we  will  obtain  a Feenberg  type  series  with  contributions 
t tl 

through  p order.  That  elimination  can  be  carried  out  explicitly  and 
yields  an  expression  not  too  wieldy  for  computational  purposes,  and 
Eq*  (3*37)  seems  preferable.  The  advantages  of  Feenberg's  scheme  are 
apparent  for  certain  cases  that  we  will  discuss  in  the  next  section. 

The  point  is  that  any  finite  approximant  Mp  to  the  Feenberg  series 
in  Eq.  (3*54),  which  consists  of  summations  over  k,U,m,  etc.,  with 
0 k — P>  an<3  injunctions  on  the  indices,  has  N orders,  each  of  which 
consists  of  a finite  sum  of  terms,  and  it  rigorously  yields  a solution 
to  a variational  secular  equation,  thus,  it  gives  an  upper  bound.  So 
far  the  index  "o"  has  not  been  specified  and  this  means  that  we  can 
obtain  any  of  the  first  p+1  states  by  setting  "o"  equal  to  the  ground 
state,  first  excited  state,  etc.  The  approximants  M were  obtained  by 
Bisshopp  in  a different  way,  but  we  want  to  emphasize  their  varia- 

tional connotations,  and  absence  of  convergence  problems  in  this  case. 

There  is  an  entirely  different  derivation,  given  by  Feenberg^] 
based  on  transformations  of  the  BW  series  which  amounts  to  a rearrange- 
ment of  the  series.  He  employed  a set  of  algebraic  identities,  origin- 

r cr  Q 1 

ally  due  to  Richards  and  by  very  interesting  manipulations,  the 
series  in  Eq.  (3*54)  is  obtained. 
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The  wave  function  itself  is  obtained  too,  with  e xpansions  for 
the  coefficients,  as  given  by  Morse  and  Feshbach.  We  shall  omit  a 
discussion  of  the  convergence  properties  of  the  series,  to  concentrate 
on  the  features  of  the  elimination  approach  that  led  to  Eq.  (3.58).  A 
secular  equation,  finite  or  infinite,  is  a condition  for  the  existence 
of  solutions  of  a system  of  homogeneous  equations.  The  elimination 
procedure  can  be  related  to  partitioning,  and  that  is  what  we  discuss 
in  the  next  section. 

3.5.  Successive  Partitioning  Technique 
and  Perturbation  Schemes 

For  simplicity  we  shall  consider  first  the  SchrBdinger  equation 

U C = E <£  (3.60) 

which  can  be  written 

1D(C  = 0,  Ml  = HI  -El  . (3.61) 

Conventional  partitioning  technique,  as  discussed  in  the  PT 
series  implies  that  we  can  partition  Eq.  (3.61)  into  four  blocks, 
corresponding  to  two  subspaces,  "a"  and  "b",  thus  obtaining  the 
equations : 

»aa  ' “aA'ka^a  " 0 (3-&0 

Cb  " - “L1  V ®a  • (3-63) 

Conventional  partitioning  technique  corresponds  to  choosing 
subspace  "a"  as  one  dimensional,  in  which  case  we  obtain: 

[Mn  - “ib  "tb  >Wci  “ 0 


which  for  Ch  f 0 implies: 


(3.6b) 


This  equation  should  be  compared  with  Eq.  (1.12).  In  fact,  Eq.  (1.12) 
is  nothing  but  the  same  equation  in  operator  form  when  the  subspace 
"b"  is  the  orthogonal  complement  of  "a".  That  is,  when  we  deal,-  in 
principle,  with  an  infinite-dimension  secular  equation.  But  going 
back  to  Eq,  (3*64),  we  can  establish  the  bracketing  theorem  for 

= f (£•) . See  PTI  and  PTII.  This  has  been  thoroughly  discussed  in 
Chapter  I.  A very  interesting  alternative  is  given  by  choosing  sub- 
space "b"  as  being  one-dimensional.  In  such  a case,  Mj^  = and 

= C^T.  Equation  (3.62)  is  now: 


r mj  T 

MI  - aN^£ 
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^NN 


<E  =0 
a 


(3.65) 


We  have  reduced  the  order  of  Eq.  (3.61)  by  one,  and  this  may  be  done 
again.  In  fact,  we  may  define  a series  of  reductions,  by  successive 
partitioning,  until  we  obtain  a one-dimensional  problem.  We  will  use 
the  notation  introduced  by  Nordling1^  .and  define  the  following: 


Ml 


(i"1)  = M(0 


aa 


,(N)  = 

M)(1)mj(1) 
ai  1a 

XX 


(3.66) 


i = N,N-1, ...,£. 


where 


Ml 


(i) 


Ml 


(i) 
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m/  ^ 
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M]^1^  m/^ 

ia  ii 


(3.67) 
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That  is,  is  a square  matrix  of  order  i,  MJ^1^  is  its  upper  left 


square  part,  of  order  i-1,  etc.  The  vectors  <E^  are  defined  as 
follows : 

= d 


d^1"^  = d^1^ 

a 


i = N,N-1, . . .,2 


(3-68) 


«<*> 


.(i) 


c(i) 

1 


where  is  the  column  vector  formed  by  the  first  i-1  components  of 

the  column  vector  d^^. 

The  successive  partitioning  amounts  to  replace  Eq.  (3.61)  by 
the  sequence  of  eigenvalue  problems: 


= 0 


i = N, N-l, . . ., 1 . 


We  are  interested  in  the  last  equation  of  the  sequence: 

M(1)c( 


^U1)  = 0 


(3.69) 


(3.70) 


,(i) 


[E-H^1)]C(1)  = ° (3-71) 

which,  for  ^ 0 is  equivalent  to: 

E = = ti(E)  (3.72) 

where  E is  an  eigenvalue  of  Eq.  (3.60).  Actually,  we  should  consider 
the  equation 

= h(£.) 

and  compare  it  with  = f(£,)  given  by  Eq.  (3.64). 


(3-73) 
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We  want  to  prove  that  h(£)  and  f(£)  are  one  and  the  same 
function,  and  that  the  only  difference  between  Eqs.  (3.73)  and  (3.64) 
is  that  the  former  corresponds  to  a triangularization  of  the  system  of 
equations  (3.63)  with  = C^,  necessary  to  compute  Eq.  (3.64).  We 
have  met  an  example  of  a system  of  linear  equations  of  this  type  before, 
namely,  Eq.  (3.44). 

Before  proving  8±  = f(£)  = h(£) . we  should  remark  that  once 
the  eigenvalue  E has  been  found,  the  eigenvector  can  be  obtained  by 
setting  - 1,  or  any  other  constant,  and  performing  the  sequence  of 
operations : 


^(i)  _ c(i-1) 


*(i) 


ta  a 
*(*■) 


M: 
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a 


The  proof  follows  from  the  fact  that  the  successive  partitioning 
as  defined  by  Eq.  (3.68)  corresponds  to  a triangularization  of  a 
of  equations.  Consider,  for  instance,  the  system  of  equations 


scheme, 

system 


MIC  = <£ 


(3.75) 


A lower  triangularization  of  this  system  amounts  to  transform 

it  into: 


Ml'c  = d7 ' 


(3.76) 


where  M 1 is  a lower  triangular  matrix,  i.e., 


all  its  elements  above 


70 


the  diagonal  are  zero.  The  connection  of  Eq.  (3.76)  with  Eq.  (3.69) 
follows  from  the  fact  that  the  ith  row  of  M » is  defined  by  the  last 
row  of  and  d \ , the  i element  of  the  column  vector  dJ  ' is  given 

by  the  last  element  of  the  vector  dl^\  defined  as  follows: 
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(3.77) 


Once  Eq.  (3.70)  is  obtained,  the  solution  of  Eq.  (3.7b)  follows : from: 
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(3.78) 


Consider  now  the  problem  MKC  = 0.  It  is  of  the  type  in 
Ecl>  (3-75)»  with  d)  = 0.  If  we  carry  the  first  N-2  of  the  N-l  steps 
indicated  by  Eq.  (3-7l)>  we  obtain: 


Ml  C - 0 (3.79) 

where  the  last  N-l  rows  and  last  N-l  columns  of  M]  are  in  lower 
triangular  form.  In  fact,  we  can  write  the  last  N-l  rows  of  Eq.  (3*79) 


as 
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“bbh, " - “bA 


where  is  lower  triangular. 


The  first  row  of  Eq.  (3.61)  can  be  written: 


Mnci  + “ibS  ‘ 0 


UQder  the  tursnsfoirmstions  thst  1 g cid  to  Eq#  (3#T9)  it  boconios! 
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For  C - 1,  it  is  simply: 


E = u(2)  + 

11  M12  C2 


From  Eq.  (3.78) 
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and  then 
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= h(E) 


(3.80) 


(3.81) 


(3.82) 


(3.83) 


This  is  clearly  the  result  from  successive  partitioning  tech- 
nique, as  determined  by  Eq.  (3-72).  Thus,  the  elimination  technique 
leading  to  a triangularization  of  Eq.  (3.61)  is  identical  to  the 
successive  partitioning  scheme.  Furthermore,  we  may  consider  the  par- 
titioning equation  (3-64)  written  in  the  following  way: 


E - Hn  + “ibS  - £<E> 

with  given  by  the  system  of  N-l  linear  equations: 


(3-84) 


M3  <E  = - Ml 
bb  b bl 


(3.85) 
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The  N-l  transformations  leading  to  Eq.  ( 3 - T9 ) » and  hence  to  Eq.  (3.83), 
do  nothing  but  carry  out  a triangular ization  of  Eq.  (3.85);  no  change 
in  the  nature  of  the  solution,  but  a somehow  simplified  form  of  the 
equations  is  obtained,  which  when  solved  and  substituted  in  Eq.  (3.84) 
yields  the  partitioning  solution  f(E).  Now,  not  only  f(E)  = h(E),  but 
in  general,  f(£.)  = h(£)  holds.  One  way  of  seeing  this  is  to  expand 
the  denominators  in  Eq.  (3.58),  obtaining  the  Brillouin  series,  which 
may  be  summed  up  again  to  yield  the  conventional  partitioning  expression. 
Of  course,  for  finite  N the  expansion  of  h(£)  given  by  Eq.  (3.58)  can 
be  transformed  back  into  a system  of  linear  equations  of  the  type  in 
Eq.  (3.44).  In  the  case  of  an  infinite  basis,  the  elimination  proce- 
dure can  be  considered  as  being  possible  in  a formal  sense.  The  equa- 
tions are  derived  by  mathematical  induction,  as  shown  by  Morse  and 
Feshbach  and  the  same  considerations  apply  to  the  successive  parti- 
tioning approach.  This  last  approach,  as  well  as  the  observation  that 
no  term  involving  Vgs  appears  in  the  sums  in  Eq.  (3.58)  implies  that 
the  matrix  elements  of  V may  be  replaced  by  matrix  elements  of  H. 

We  went  into  the  successive  partitioning  scheme  with  consid- 
erable detail  because  it  illuminates  some  aspects  of  the  Feenberg 
approach.  From 

f(£)  = h(£)  ' (3-86) 

it  follows 

0 >f(E.)  = h-(e)  . (3.87) 

The  importance  of  this  result  is  quite  apparent  . The 
Feenberg  perturbation  expansion,  given  by  Eqs . (3.58)  and  (3.59)  is  a 
bracketing  function,  namely  that  and  £ bracket  one  eigenvalue 
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of  H.  If  one  has  a finite  approximant,  i.e.,  with  a finite  number  of 
orders  included,  the  solution  to  E = f(E)  is  bracketed,  but  this  is  not 
the  exact  eigenvalue,  but  an  upper  bound  to  it.  Actually  one  is  not 
interested  in  finite  approximants  to  the  Feenberg  series,  but  in  the 
infinite  expansions.  The  reason  for  it  is  that  they  assume  particularly 
useful  forms  under  special  circumstances,  and  these  can  be  discussed  in 
terms  of  the  properties  of  successive  partitioning.  As  pointed  out  by 
Nordling  ^ , successive  partitioning  of  a triple  diagonal  matrix  has 

a very  simple  form.  By  triple  diagonal  matrix,  we  mean  a matrix 


«1  P2  ° 

P2  oc2  P3  0 


° PN 


(3.88) 


For  any  symmetric  real  matrix  it  is  possible  to  define  a unitary 
transformation  which  brings  it  into  triple  diagonal  form.  Many  matrices 
corresponding  to  physical  problems  are  already  in  this  form:  systems 

with  nearest  neighbour  interaction.  There  are  equations  where  the 
nature  of  the  perturbation  is  such  that  it  has  a triple  diagonal  form 
automatically,  i.e.,  the  Mathieu  prblem,  the  anharmonic  oscillator 
with  V = X.X.  Actually  the  matrix  corresponding  to  the  Mathieu  problem 
is  triple  diagonal"  but  with  a string  of  zeroes  between  the  main 


diagonal  and  the  two  parallel  strings  of  non-vanishing  elements ; this 


is  due  to  the  fact  that  ^nm=Q  unless  m - n+2  or  m — n.  Consider  then, 
in  general,  a triple  diagonal  matrix  with  k strings  of  zeroes  separating 


the  strings  of  non-zero  components.  We  can  generalize  the  scheme 


partitioning,  in  terms  of  a continued  fraction  for  an  infinite  matrix, 
or  in  terms  of  a truncated  continued  fraction  for  any  finite  problem. 
This  corresponds  to  an  extension  of  Nordling's  scheme.  The  observation 
of  the  continued  fraction  representation  leads  torcompact  derivation 
of  several  properties  of  this  type  of  solution. 


strings  of  zeroes  between  the  diagonal  and  the  two  parallels  of  non- 


defined  by  Eq.  (3*66)  and  obtain  an  explicit  solution  via  successive 


In  terms  of  Eq.  (3.66)  define: 


(3.89) 


For  any  triple  diagonal  matrix  of  the  types  discussed  above, 
the  matrix  ^ has  only  one  non-vanishing  element.  If  there  are  k 


zero  matrix  elements,  the  non-vanishing  components  of JD^ ^ will  be: 


i = N-l,N-2 


> • • • > 


k+1.  (3.90) 


To  make  this  clear,  we  shall  use  the  following  notation 


T = 


ce1  0 0 ^ 0 0 

0 «2  0 0 0 0.. 


(3.91) 


0 0 a3  0 0 0g  . . . 

0 0 a,  0 0 . . . 


where  we  have  depicted  the  case  k = 2. 


75 


Then  Eq.  (3*90)  can  be  written: 

D(i)  _ ^i+1 

1+1 

where  c/P  is  given  by: 


i = N-l,N-2, 


k+1 


(3.92) 


cfc-V  -a^  = 

i-k  i-k 


a. 

1- 


-D(0 

i-k, i-k’ 


i = N-l.N-2, . . .,k+l. 


(3-93) 


Equation  (3-93)  states  simply  that  when  we  carry  the  successive 
partitioning  corresponding  to  Eq.  (3*66),  if  the  matrix  is  triple  diag- 
onal, only  one  diagonal  element  in  the  subspace ja_  is  going  to  be  changed. 
The  process  indicated  by  Eq.  (3*93)  is  carried  until  we  obtain 


E = h(E)  = . 


(3.9*0 


What  we  want  to  stress  is  that  Eq.  (3*9*0  can  be  used,  in  con- 
junction with  Eqs . (3*93)  and  (3-92)  to  construct  an  explicit  solution 
by  "back"  iteration: 

From  Eqs.  (3-94)  and  (3.93): 

e = h(E)  = a[k+lS>  = ax  - D^k+1) 
and  using  Eq.  (3-92) : 


E = 


a. 


&2)-e 


= a,  + 


fi2 
Pk+2 


h +IZ  + n(2k+2) 

E °k+2  °k+2,  k+2 


= a,  + 


q2 

Pk+2 


'1  ' 2 

E-a  + — 

^+2  + (2k+3) 

a2k+3  E 


(3.95) 


Let  s consider,  for  example,  the  case  k — 0.  This  is  the  conventional 
triple  diagonal  matrix  case: 


0 . • . 


Pg  Qg  * " * 

0 p3  a3 


(3.96) 


The  first  eigenvalue  is  given  by: 


E = 


°1  + 


^3 

E-a2 L 


(3.97) 


e-ol 


e - a,  - 


As  pointed  out  by  Nordling^^ , any  NxN  symmetric  real  matrix 
H can  be  transformed  into  the  form  of  Eq.  (3*9^)  by  a unitary  trans- 
formation. Hence  Eq.  (3-97)  is  the  explicit  representation  of  the  first 
eigenvalue.  Of  course,  suitable  modifications  of  Eq.  (3.97)  hold  for 
other  eigenvalues.  Equation  (3*97)  is  valid  for  an  infinite  matrix  of 
the  form  in  Eq.  (3-9^)  in  which  case  it  provides  a continued  fraction 
type  of  solution.  In  either  case,  we  want  to  stress  the  importance  of 
the  continued  fraction  solution,  obtained  by  the  "back"  iteration  in 
Eq.  (3-95).  b ecause  it  provides  fast  numerical  convergence  and  because 
the  formal  properties  of  continued  fractions,  like  recurrence  relations, 
etc.,  give  a simple  explanation  of  the  properties  of  triple  diagonal 
matrices,  their  determinants  and  eigenvalues,  and  this  is  certainly 
an  alternate  path  to  the  one  discussed  by  Nordling. 
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Applications  of  the  Feenberg  scheme  are  interesting  when  the 
matrices  are  triple  diagonal,  but  in  these  cases  the  continued  frac- 
tion solutions  are  a natural  tool,  which  follows  from  successive 
partitioning.  When  the  matrices  are  not  triple  diagonal,  the  possi- 
bility of  implementing  Feenberg's  scheme  relies  on  finite  approximants , 
which  are  related  to  the  variational  principle,  suggesting  that  other 
schemes  may  be  equally  or  more  useful,  in  fact,  the  truncations  of 
Tq  discussed  in  Chapter  V and  the  applications  of  the  techniques  of 
PTX  are  more  interesting  in  terms  of  usefulness.  On  the  other  hand, 
the  continued  fraction  solutions,  for  problems  with  nearest  neighbour 
interactions  should  prove  very  useful. 


CHAPTER  IV 


RECIPROCAL  VARIATIONAL  PROBLEMS 
AND  PERTURBATION  THEORY 

The  implementation  of  conventional  perturbation  theory  requires 
the  computation  of  the  individual  orders,  both  in  BW  and  RS  theory. 

Exact  evaluation  involves  an  infinite  sum  and  an  integration  over  the 
continuum  when  using  a resolvent  technique;  solution  of  an  inhomoge- 
neous equation  is  the  alternative,  not  easily  attainable  in  most  cases. 
Thus  the  need  for  approximate  solutions,  which  are  of  variational  type. 
These  may  be  very  successful,  but  we  should  emphasize  that  only  through 
computation  of  upper  and  lower  bounds  to  the  individual  orders,  some- 
thing may  be  said  about  the  nature  of  the  approximation. 

[59] 

Hylleraas  was  the  first  to  derive  a variational  expression 
for  the  second  order  energy  in  RS  theory,  giving  an  upper  bound  to  it. 

In  fact,  as  the  inhomogeneous  equations  that  characterize  the  different 
orders  in  RS  and  BW  theory  can  be  related  to  variational  principles, 
it  turns  out  that  the  even  order  energies  have  upper  bounds  given  by 
functionals.  The  minimization. yields  an  upper  bound  to  the  even  order 
energy  and  a function  that  can  be  used  to  compute  approximations  to 
the  other  orders.  Knight  and  Schefcr^60^  exploited  this  in  the  He  series. 
Midtal  discussed  the  orthogonality  properties  that  explain  their 
success.  Prager  andeHirschfelder obtained  an  expression  for  a lower 
bound  to  the  ground  state  second  order  perturbation  energy  in  the  RS 
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scheme  and  this  was  generalized  by  Sando  and  Hirschfelder^*^  . Meath 

r i 

and  Hirschfelder L obtained  variational  upper  bounds  in  the  BW  scheme, 
similar  to  the  Hylleraas  principle.  Sinanoglu^^  stated  an  orthogo- 
nality requirement  towards  the  lower  unperturbed  eigenfunctions  which 
allows  extension  of  the  Hylleraas  principle  to  excited  states,  and 
this  was  carried  out  for  the  BW  case  in  reference  63. 

Hirschfelder  and  Prager  derived  their  formula  by  analogy 
with  electromagnetic  theory.  We  shall  introduce  the  geometric  lan- 

fO-l 

guage  of  Synge1  , in  the  case  of  positive  metric,  and  then  show  that 
special  choices  lead  to  the  Hylleraas  principle  and  to  the  Hirschfelder 
and  Prager  formula  in  a natural  way. 

4.1.  Reciprocal  Variational  Problems 

Consider  a linear  vector  space  -A,  with  vectors  £.  We  don't 
specify  their  nature  at  all;  they  may  be  functions,  or  sets  of 
functions.  Consider  a Hermitian  positive  metric  in  the  space  -A-, 
defined  by 

Q (£>£.)  > 0 (4.1) 

where  the  equality  sign  holds  if  and  only  if  £ = 0. 

Consider  two  linear  subspaces  of -A.  ; Z and  JX.  We  will  say 
that  they  are  orthogonal  if  and  only  if  for  all  £ e Z and  all  cn  e IX 

Q (0,05)  = 0 . (4.2) 

If  we  interpret  Q (£»£.)  as  the  square  of  the  norm  of  £,  a "Pythagorean 
theorem"  is  fulfilled: 

Q(cn-£,cn-£)  = Q(ajb£,ajfa)  = Q(to>cn)  + Q(£>£)  . (^*3) 


8o 


Within  this  context  define  two  variational  problems: 
i)  Given  two  vectors  and  c^.  a linear  subspace  -/"l.  and  its 
orthogonal  complement  Z,  minimize 

Ve.)  = Q(£-£0>£.-£0)  (4.4) 

with  the  constraint  that  £ - £q  should  belong  to  -0-  . 

Il)  For  the  same  PQ» » -O-  and  Z as  before,  minimize 

d2(E.)  “ Q(£-£0»£-£0)  (4.5) 

subject  to  the  constraint  that  £ - £q  should  belong  to  Z. 

These  two  variational  problems  have  some  properties  of  interest, 
easily  derived.  Let  £ be  the  solution  of  Problem  I,  i.e.,  it  minimizes 
Eq.  (4.4).  Then  consider  a specific  £,  of  the  form: 

£L  = u + (4.6) 

for  real  e and  any  e £1.  Substitute  Eq.  (4.6)  into  Eq.  (4.4),  and 
set  the  first  variation  of  e equal  to  zero. 

dq(£q)  = Q(Ei’Ei)  = Q(E+e0"3o»E+e0"2o) 

Sdl^l^  = ‘ 2Q(cu,£q)  = 0 

which  may  be  written 

Q(£>£0>^)  = 0.  (4.7) 

This  consists 
Theorem  4.1: 

The  necessary  condition  for  Problem  I to  have  a solution  is  that 
£ “ Hq  should  belong  to  Z.  Notice  that  this  is  precisely  the  constraint 


of  Problem  II. 
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Theorem  4.2: 

_u  is  unique.  The  proof  consists  in  assuming  two  solutions  u^ 

and  Ug.  Then  u^c^  e 2,  and  e 2.  Hence,  u^-Ug  e 2.  But  the 

constraint  of  Problem  I imposes  u -p  e ./l,  u -p  e jTX  and  hence 

“1  o —2  o 

e ^is  means  that  u^-ui^  belongs  to  the  intersection  of  2 and 

A , which  are  orthogonal  subspaces,  and  therefore  contains  only  the 
zero  vector.  Therefore  u^  = u^  . Q.E.D. 

The  following  two  theorems  have  entirely  similar  proofs: 

Theorem  4. 3 : 

The  necessary  condition  for  Problem  II  to  have  a solution  is 
that  v - £o  is  in  A , where  v is  the  vector  that  minimizes  Eq.  (4.5). 
Notice  that  this  is  the  constraint  of  Problem  I. 

Theorem  4.4: 

v is  unique. 

The  significant  result  in  this  context  is: 

Theorem  4 : 

_u  = v . 

Proof: 

u - e 2 (Theorem  4.1). 

v - e 2 (Constraint  of  Problem  II). 

Hence 

u - v e 2. 

S imilarly 

v - ^ e A (Theorem  4.3) 

H - £q  e jT\-  (Constraint  of  Problem  i). 


82 


Then 

_u  - v e 

and  it  follows  that  u - v must  be  the  zero  vector.  Q.E.D. 


Theorem  4.6: 


Let  be  the  minima  of  Problems  1 and  IX,  respectively. 

Then:  ^ + d.,  - 

Proof : 


dl  + d2  - QCii-^iL-c^)  + • 

By  the  Pythagorean  relation  in  Eq.  (4.3): 


dl+  d2  - Q(u-£o-u+£o,u-£o-u+£o)  = Q^-c^-c^).  Q.E.D. 

The  sum  of  the  minima  of  the  two  reciprocal  problems,  I.  and  II,  is 
then  exclusively  determined  by  £^  and  £q.  We  say  reciprocal  because 
the  constraint  of  one  problem  is  the  Euler  condition  of  the  other  and 
vice  versa.  The  following  theorem  also  holds,  and  is  the  key  to  the 
applications : 


Theorem  4.7: 

The  sphere  of  radius  I £-£ 

l~ 


about  the  point 


cal  locus  of  £. 

The  proof  follows  from  consideration  of 


£+£  is  a geometri- 

“5“ 


£+£  £+£ 

4Q(n 2~’- 2~ ) = Q (£"£+£"£>  £"£+£"£)  = Q (£"£>£-£.) 

+ Q(£-£»£-£)  = Q(£-£-(£-£)  »£-£-(  u-£)  = Q(£-£,£-£) 


C^-8) 


The  third  line  in  Eq.  (4.8)  follows  from  the  orthogonality  of  £-£  and 
£-£,  and  the  fourth  line  too,  bu  successive  application  of  the 
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"Pythagorean  Theorem"  in  Eq.  (4.3).  But  Eq.  (4.8)  is  precisely 
theorem  4-7>  as  it  can  be  seen  by  dividing  through  by  four: 


£+£ 

£-£ 

U " T 

2 

(4.9) 


We  can  restate  the  theorem.  The  common  solution  to  the  recip- 
rocal Problems  I and  II  deviates  from  the  arithmetic  mean  of  two  trial 
vectors  £ and  £,  satisfying  the  constraints  for  Problems  I and  II, 
respectively,  by  an  amount  equal  to  half  the  distance  between  the  two 
vectors  £ and  £.  This  is  nothing  but  a version  of  the  theorem  of  Thales, 
as  it  can  be  seen  from  Figure  .5. 

Equation  (4-9)  constitutes  an  estimate  of  the  norm  accuracy 
of  the  trial  vectors,  p and  g;  for  the  reciprocal  Problems  I and  II. 
Theorem  4.6,  has  the  important  consequence  that  an  upper  bound  for  one 
problem  constitutes  a lower  bound  to  it's  reciprocal,  since  the  sum  of 
the  two  minima,  d^  + d , is  given  in  the  specification  of  the  problems: 

di  + d2  = = a • (4a°) 


Let  d^>  d , then 


. > -V/ 


d 


2 


a - dx  > a - dL  = d2 


Thus,  an  upper  bound  to  the  direct  problem,  d^,  leads  to  a lower  bound 
to  the  reciprocal  problem,  d^. 

Instead  of  exploring  the  many  other  cases  discussed  by  Synge 
in  his  book  [8],  we  will  consider  a quantum  mechanical  application  of 


the  present  formulation. 
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4.2.  Inhomogeneous  Equations  in  Rayleigh-SchrBdinger 
Perturbation  Theory 

It  is  very  well-known  that  the  different  orders  in  RS  pertur- 
bation theory  are  defined  by  a series  of  inhomogeneous  equations.  See 
PTIX,  and  reference  6. 

n 

(Ho"Eo^n+l  = ^ £ n+1  -k^k  “ n=0,  1’  * * * (^-H) 

k=0 


with 


:k+!  = (^olVKk> 


oo 

-.♦I 


k=l 


- I 


(4.12) 


k=0 


where  we  have  split  the  Hamiltonian  H into  an  unperturbed  part,  H , 

and  a perturbation  V.  Further  we  assume  H \|r  = E \|r  . 

o o o o 

Solutions  of  the  "low  orders"  in  Eq.  (4.11)  and  Eq.  (4.12) 
have  been  obtained  in  several  problems  of  interest  by  the  Wisconsin 
group1  5 . One  very  important  property  is  the  result  originally  due 
to  Wigner^^,  if  we  know  the  approximate  wave  function  through  order 
n,  the  energy  can  be  computed  through  order  2n+l. 

Let  us  consider  the  first  equation  in  Eq.  (4.1l),  and  assume 
that  we  deal  with  the  ground  state  energy 


(VEo>h  - - (v-^)*0 


(4.13) 
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where  the  first  order  energy  is  given  by 

ei  ‘ <*0|v|*o>  - v 


oo 


(4.14a) 


The  second  order  energy  is 


^2  " 


(4.14b) 


We  choose  \|r^  orthogonal  to  \Jr  , to  satisfy  an  intermediate 

normalization  condition.  We  can  construct  a variational  principle 

for  Eq.  (4.14)  if  we  notice  that  H - E 2.  0,  i.e.,  is  a positive 

o o r 

operator.  Hence  it  can  be  considered  as  defining  a positive  metric  in 
a space  A.  In  fact,  let  us  identify  the  vectors  £ with  functions  \|/; 

£q  = t0  and  £o  = 

The  direct  Problem  I,  of  Eq.  (4.4)  now  reads: 

d]_(^)  = Q(^"^1»t“t1)  • (4.15) 

The  minimum  will  be  attained  by  % = \|r  , the  solution  of  Eq. 
(4.13).  What  will  be  d^,  the  minimum  of  d^(\J/)? 

This  is  seen  by  writing  Eq.  (4.15)  explicitly:: 


dx(T)  = (f-t1jHo-Eo|'J-t1) 

= <f|Ho-Eo|?)  +2(J|Ho-£1|to)  - 62 


(4.16) 


where  we  have  used  the  Hermitian  property  of  H^,  and  Eq.  (4.13).  The 
minimum,  attained  by  \| r = is  0,  as  it  can  be  seen  by  further  use  of 
Eq.  (4.13).  We  have  then: 

*2  ~ <?|Ho-EJt$  + 2<IjV-£1|t0> 


(4.17) 
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This  is  nothing  but  the  variational  principle  ■ for  the  ground 
state  second  order  energy  due  to  Hylleraas^^  > but  derived  within  the 
geometric  formulation  of  Synge. 

This  would  not  be  of  much  interest  if  it  were  not  for  the 
reciprocal  problem. 

Let  £ = 0.  Then  Eq.  (4.5)  reads: 

dp(0)  = Q(0-V$-O  = Q(0.0)  = <0lH  -E  |0)  (4.18) 

^ u u O O 

with  the  constraint  0-\J/^  e Z.  The  constraint  amounts  to  impose 
Q(w~i|/^>ty)  for  all  \|/.  This  would  amount  to  restrict  0 to  be  equal  to 

but  in  fact  we  have  more  freedom  than  that.  Consider  first  Eq.  (4.18). 
The  minimum  is  given  by  ' 


di  + d2  = Q(t0-t1»t0-^1)  = Q(t1»^1) 

~ ^JHo"Eo^l^  ""  -<^1 1 Hl"S  I^o>  ~ ~ ^“2 


Then 


and 


d2  " ^2 


^2  > -d2(0)  • 


(4.19) 


Therefore,  if  we  could  implement  Eq.  (4.l8),  we  would  have: 

-d2(0)  6 €2  £ ^(V)  • (4.20). 

The  Dalgarno-Lewis procedure  to  determine  solutions  to 
inhomogeneous  equations  of  the  type  in  Eq.  (4.13)  can  be  applied  in 


this  case. 
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Let  0 be  written  as  0 = f\|r  . Then  we  have, 

o 

( H0-E0 )0  = (K+V-EQ)fto  (4.21) 

where  we  have  set  = K + V;  K is  the  kinetic  energy  operator,  which 

we  are  going  to  represent  symbolically  by  a 3N-dimensional  Laplacian: 

1 2 

K - - — V . V is  the  potential  energy  in  the  unperturbed  Hamiltonian 

and  we  will  assume  it  is  a function  of  the  coordinates  of  the  N particles. 

As  shown  by  Hirschfelder L , we  can  rewrite  Eq . (4.2l)  by 

subtracting  the  expression  f(H  -E  W from  both  sides: 

v o o'  o 

(Ho-Eo)0  = (K+V-Eo)f^  - f(K+V-Eo)to 
= [K,f]^  = |[V2,f]t0. 

The  square  brackets  denote  a commutator.  Furthermore,  we  can  write 

0(VEo)8)  = 

- - . (4.22) 

This  follows  by  standard  vector  analysis: 

♦0[v2.f]t0  - * v2(£*0)  - <l<0fv\ 

= \|r  V'V(fty  ) - \lr  fV  \|/  = Air  V*(fV\l/-  Vf)  - il/-  fV2\|/ 

To  ' o o o o ' o o o o 

= 2\|/QVf  'V\J/o  + ^2f  = V(f^Vf)  . 

Therefore  Eq.  (4.18)  is  transformed  into: 

d2(0)  = d2(f)  = - | J fV.(t2Vf)(dv) 

= l/ *o^f(dv)  " |fv.(^2fVf)(dv)  . 
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The  divergence  term  vanishes  if  we  impose  suitable  boundary  conditions 
on  0.  This  means  that  Eq.  (4.18)  can  be  replaced  by: 

d2(^')  = 2 (4.23) 

where  q is  a vector  function  that  satisfies  an  equation  derivable 
from  the  inhomogeneous  Eq.  (4.13).  Let  0 = Y.  = f\|r  . 

(H  -E  )ft  = - (V-£\)t 
v o o o v l/Yo 

By  the  same  arguments  as  before: 

[V2,f]^  = 2(V-61)lifo 
or 


v-(t^f)  = 2^/o(v-61)to 


Let  Vf  = q.  Then  if  q satisfies: 


V-(^q)  = 2t0(V-^)t0 


then. 


d2C5)  - 


^q.q(dv) 


gives  a lower  bound  to  £^: 


- d2(q)  ' ^2  ~ diM 


(4.24) 


(4.25) 


(4.26) 


In  fact  we  have  shown  through  a first  integration  of  the 
inhomogeneous  equation  that  it  is  possible  to  formulate  the  reciprocal 
problem  and  hence  get^a  lower  bound  formula.  It  is  clear  that  an 
admissible  trial  vector  for  Problem  II  has  to  satisfy  Eq.  (4.24),  and 
this  is  by  no  means  easy  to  do  in  general.  Hirschfelder  and  Prager 
used  it  for  the  case  of  the  polarizability  of  atomic  hydrogen  with 
very  good  results  as  compared  with  other  approxinate  methods.  Its 
applicability  in  the  future  dependes  then  in  solving  Eq.  (4.25)  for 
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complicated  potentials  V.  The  extensions  to  higher  even  orders  have 
been  given  by  Sando  and  Hirschfelder^2^ . 

Our  interest  is  rather  in  studying  the  forcefulness  of  this 
geometric  approach,  and  we  can  consider  the  direct  problem  again, 
within  the  context  of  perturbation  theory,  but  for  the  excited  states. 
The  pertinent  inhomogeneous  equation  is  Eq.  (4.13)  again,  but  modified 
as  follows: 


(vEi“Mn)  - - (v-4n))i'in)  • 

Xn  order  to  construct  an  upper  bound  to  the  second  order  energy^ 
we  have  to  insure  that  the  metric  defined  by  the  kernel  (H  -E^nb 


(4.27) 

(n) 

2 * 

is 


positive.  This  is  achieved  by  restricting  the  space  of  trial  func- 
tions to  be  orthogonal  to  the  \|/k),  k=0,  1,  . . . , n-1,  i.e.,  to  the 

lower  unperturbed  eigenfunctions.  Under  this  conditions,  the  direct 
Problem  I applies  and  the  extension  of  the  Hylleraas  principle  to  the 
R-S  excited  state  second  order  energy  follows  trivially,  thus  proving 

the  Sinanoglu  results:  \Jr  = V*  [V.  /(E^-E^bli/1^  + 

n .‘-J  kn  v o o /JTo  v 

k=o 

This  shows  that  the  initial  care  in  introducing  the  geometric 
approach  has  the  consequence  of  providing  unified  and  simple  derivations 
of  the  lower  and  upper  bound  formulae  in  R-S  theory.  It  is  certainly 
of  great  interest  to  extend  it  to  other  problems,  and  by  no  means,  this 
has  been  done  to  any  extent. 


4.3.  Inhomogeneous  Equations  in  Br illouin-Wigner 
Perturbation  Theory 

There  is  an  entirely  analogous  derivation  for  upper  bounds  to 
the  B-W  even  order  energies.  In  this  case,  for  £ < E°,  the  operator 
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Hq  - £ is  positive  definite  in  the  orthogonal  complement  to  \Jr  , and 
the  subspace  of  £ vectors,  i.e.,  the  ^ functions  are  to  be  chosen 
orthogonal  to  i|f  . 

Consider  the  inhomogeneous  equation  for  the  ground  state  first 
order  wave  function  in  the  B-W  scheme 


(H0-&)t1  + (V-^)^  = 0 . (4.28) 

The  second  order  energy  is  defined  by 

^2  = <^ofV|V  *'  (4-29) 


As  before,  let  £ - \jr,  £^  - and  £ = \|/  . The  direct  problem 

in  Eq.  (4.4)  yields: 

% 

d1(t)  = Q(t-t1»t_t1)  = H0-f| 


By  use  of  Eq.  (4.28) : 

dx(t)  = <^/ho-£|T>  + 2<t|v-  v-C1h0> 

= 4^fHo-£|^  + 2(tlv|to>  - 62  . 


(4 .30) 


The  minimum  of  d^(\|r)  is  d^  = 0,  for  \|r  = i(r-.  Therefore, 

6 2 = OlH0-£|V>+  2<\|r|Vl\|ro>  (4.31) 

r>S 

for  all  \|r  orthogonal  to  , and  this  is  the  desired  variational  prin- 
ciple, whose  generalizations  to  higher  orders  and  to  excited  states 
follow  the  same  lines  as  the  previous  section.  On  the  other  hand,  the 
reciprocal  problem  cannot  be  implemented  in  a simple  fashion,  because 
a first  integration,  by  quadratures,  of  the  B-W  equations  is  not 

ki  [67] 

available 


CHAPTER  V 


« PARTITIONING  TECHNIQUE  AND  THE  METHOD 
Vi  0F  INTERMEDIATE  PROBLEMS 

la,  Chapter  III  we  have  discussed  a truncation  of  T , denoted 
; o’ 

by  Tj(p).  We  did  not  discuss  at  all  the  truncation  T (p),  which  for 
\ 

our  purposes  £n  obtaining  lower  bounds  is  far  more  reaching.  We'll 

i . . X 

relate  it  to  the  method  of  intermediate  problems  and  point  out  some  of 

i 

the  limitations  of  the  latter.  Finally,  we'll  give  some  by-products 

><  4 

of  the  partitioning  approach  to  the  intermediate  Hamiltonians,  like 
\ 

lower  bounds  to  the  E^’s  and  improved  energy  estimates. 

5.1.  The  Truncation  T^fp) 

We  will  start  with  a single  reference  function,  an  eigenfunc- 
tion to  H^,  the  unperturbed  Hamiltonian.  Then  the  following  formulae 
are  valid,  as  shown  in  Chapter  I: 


H i|f°  = e!V 
ok  kvk 

(5-1) 

®k  - \\ 

(5-2) 

£i,k " Ek  + ^kMO 

(5-3) 

) 

t = (V_1-T  )_1 
? x o 

(5*4) 

/ 

T 

o 


(5-5) 


and  if  E^<  £ , ^ < E^  follows. 
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Assume  that  £ < E°+^.  Then,  from  Eq.  (I.30): 


T = 
o 


°o  r~> 

Zz±+L 


o°  p 


0 


£ 


l 


, , f _ rO 

Jh^k  £ U=p+1  £ 


o / °n 

£ - E9  t -E°+1  L-j  £ 

ik  £ P+1  i-p+i 


00 

E 


•z 

£*k 


t-n 


s + d. 


^rl1  - 1 °i]  ■ T°(p) 

£= o 


(5-6) 


Let's  assume  that  £ < E°  and  that  we  consider  the  ground 


state.  Then  0 ? Tq  TQ(p).  For  positive  V,  V ^ > 0 and 
0 < V_1-t  h t"1  & V-1-T  ( p)  - t'L(p)  . 


(5-7) 


According  to  theorem  1.6,  it  follows  that 


0 < t( p)  < t 


(5.8) 


and  Eq.  (5-3)  way  be  replaced  by 


t-l  - Eo  + <*„l|:v'i-T0(p)l'1|*0> 


-1 


I -1 1 


(5-9) 


where  we  have  a bracketing  of  the  ground  state  eigenvalue: 

2 Eo  < <£.  . (5-10) 

It  is  interesting  to  see  that  there  is  no  need  to  expand  the 
inverse  in  Eq.  (5.9)  io  a power  series  to  obtain  a Brillouin-Wigner 
solution.  We  can  solve  Eq.  (5*9)  to  a^l  orders  by  introducing  an 
inhomogeneous  equation  approach  entirely  similar  to  what  was  done  in 
Chapter  III.  Let  X be  a function  defined  by 

[v"1-To(p)]”1\Jfo_  = X 
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or 

[V^-TjpJJX  - tD  . (5.11) 

We  just  have  to  replace  Eq.  (5 .6)  into  Eq.  (5.H)  and  perform 
some  algebraic  transformations:  ' 

P 

+ Z +^(e-E°+1)-1-(t-E°)-1(i-6oi)]<1|,°|x)  . 
l~  O 

We  can  multiply  by  i|f^[V  -(£-E°+^)  ^ s \|/-°v,  and  integrate,  thus 

obtaining 

<*k|v|*o> 

- <(,°|x>+  £ (*“|v|*°>[(£.-E“+1)-1--(£.-E°)1l-Soi:)]<t°|x>  . 

!>- O 

Let  us  denote  <\|r£  I v |i|rj  > by  v^,  (^°|x)  by  C^,  and 
[(£  - Ep+l)  “ (^"Ep  by  d^*  Then  Eq.  (5.1l)  is  finally 

transformed  into  the  system  of  equations: 


P 

vko  vk JJ  djl  ^Ji*  _ °»^»,,,»P»  (5  *12) 

i= o 

We  can  express  Eq.  (5. 12)  in  more  compact  form  by  using  matrix 

notation.  Let  V be  a column  vector  with  the  v,  ’ s as  components;  let 
o ko 

<C  be  a column  vector  with  the  C ’ s as  components  and  let  A be  a square 
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matrix,  with  (A  ) ^ v^d^.  Then  Eq.  (5. 12)  can  be  written: 

= C + AC  = (1+A)C  . (5.13a) 

This  means  then: 

a = (n-A)_1vo  , (5.13b) 

and 

P 

(0X)  3 [(1+A)  \]0  = Yj  C(1+A)"1]oiVio  (5-14) 

i=o 

is  the  explicit  form  of  (\|/Q  | t(p)  j ijr^  }.  It  should  be  pointed  out  that 
the  inverse  matrix  in  Eqs . (5.I3)  and  (5-14)  does  not  have  to  be  com- 
puted, because  all  we  need  is  the  solution  of  the  system  of  Eq.  (5. 12). 

The  discussion  of  the  solution  of  this  system  of  equations  parallels 
the  discussion  of  the  system  of  equations  that  arose  in  Chapter  III 
and  led  to  Feenberg's  perturbation  theory.  The  coefficients  are  slightly 
different  and  there  are  new  integrals,  the 

VM  ■ <<ltv'1-(«-E°+1)-1]-1l^>  (5.15) 

to  be  evaluated.  The  compensation  of  this  additional  complication  is 

that  for  <£  ^ E^  we  obtain  a lower  bound  to  the  ground  state  energy. 

The  fact  that  Eq.  (5-9)  itself  is  a bracketing  function,  i.e.,  it  has 

negative  slope,  suggests  the  examination  of  the  point  = £ = E1. 

It  is  quite  clear  that  E'  is  bracketed  by  the  pair  t J,  £ ; but  this 

does  not  imply  a simultaneous  bracketing  of  E . Is  E ' a lower  bound 

O o 

to  Eo?  The  answer  to  this  is  yes,  and  it  is  given  by  the  theory  of 
intermediate  problems,  as  we  will  see  in  the  next  section. 
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Could  we  apply  the  same  treatment  to  excited  states?  It  is 
easy  to  see  that  a system  of  equations  entirely  similar  to  Eq.  (5. 12) 
could  be  obtained,  but  we  have  not  proved  that  Eq.  (5.8)  applies.  In 
fact  we  will  see  that  Eq.  (5*8)  is  too  strong  a condition  and  that  it 
will  have  to  be  modified.  So,  even  though  we  have  not  proved  that 
(\!c(p)  |\|r£>  is  smaller  than  1 1 ) , we  will  get  an  explicit  form 

for  it:  it  is  just  a modification  of  the  system  of  Eq.  (5. 12): 

P 

Vmk  = Gmk  + ^ Vm£d^k)CJk'  U>m>  k=°>  l>  • • • > P ) (5-16) 

1,-0 

where  we  defined  C , = (#° I X ) 

mk  Ym  k 

and  d^k)  - (C-Ep’V^)]  ; 

the  v s are  the  same  as  in  Eq . (5.I5). 

The  solution  of  this  system  of  equations  leads  to  an  explicit 
form  for  Eq.  (5.3): 

e.]  'Ek  + Ckk-  (5-l7) 

1,  k 

In  Eq.  (5. I7)  we  can  impose  the  condition  £,,=£=  E ' . In 

1,  k k 

Chapter  III,  we  showed  that  this  condition  led  to  a secular  equation. 

It  can  be  _shown  that  the  same  is  true,  and  we  examine  this  question 
in  the  next  section;  the  only  point  we  want  to  emphasize  here  is  that 
the  truncation  TQ(p)  implies  that  there  is  no  need  for  a power  series 
expansion  of  the  inverse  operator  [V  ^-TQ(p)]  ^ in  order  to  compute  its 
matrix  elements,  thus  avoiding  the  convergence  and  further  truncations 


of  the  BW  series. 
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9.2.  Intermediate  Problems 

The  method  of  intermediate  problems  was  developed  by  the 

[31 

Weinstein,  Aronszajn  school  - In  particular  Bazley  and  Fox  made 
extensive  contributions  to  it,  within  the  context  of  quantum  mechan- 

r 34 1 

ical  applications . , In  terms  of  the  theorems  introduced  in  Chapter  I 

it  is  simple  to  define  a sequence  of  intermediate  problems  corresponding 
to  a specific  eigenvalue  problem.  Let  A be  a self-adjoint  operator 
bounded  from  below.  Then 


Auk 


3kUk 


(5-18) 


is  the  eigenvalue  problem  under  consideration.  Suppose  that  we  can 
construct  a scries  of  operators  A^p^that  satisfy: 

A(0)^  A<lJ  - ...  £ A(P)  a A(P+1>^  ...  A . (5.19) 

For  each  of  them  there  is  an  eigenvalue  problem 


A 


(P)U(P)  = a(P>u(p) 


(5-20) 


Then  theorem  (1.2)  states  that 


(o) 

ak  * 


*4p)^4p+1) 


a^,  in  order 


(5-21) 


Thus,  to  obtain  lower  bounds  to  the  eigenvalues  of  Eq.  (5.18) 

we  need  the  requirements:  we  should  be  able  to  construct  the  operators 

A^^;  they  should  be  solvable,  and  we  should  keep  track  of  the  ordering 

Eq.  (5.2I).  The  last  requirement  is  important  in  connection  with  lost 

eigenvalues.  As  the  kC^  eigenvalue  of  A^^  is  a lower  bound  to  the 
t tl 

k eigenvalue  of  A,  it  is  important  to  make  sure  that  all  the 


97 


eigenvalues  of  A^P^  are  obtained,  otherwise  the  ordering  can  not  be 
implemented. 


If  A - A + B,  where  A has  a known  spectral  resolution, 
o o 

can  define  a series  of  truncations  A^P^: 

o 


we 


>(p) 


k=o 


k=o 


(5-22) 


u rs°  0 \ / O . O OO 

where  0.  = I u,  Mu,  ; A u,  = a,  u,  . 
k k k o k k k 

If  A is  bounded  from  below,  and  a°  £ a°  £ ...  6 a°i  ° 


o 1 


it  follows  that 


— a £ • • • , 

p p+1 


± A^  . . . £ A^P^  A^P+1^  . . . £ A . (5.23) 

o o o o o 

Then  we  can  define  a sequence  of  operators  A^  = A^  + B which 
satisfies  the  relations: 


A(o)*  A<l)  S ...  - A<pM  a{p+1M  ...  - A . 


(5-24) 


This  means  that  if  we  can  solve  the  eigenvalue  problem  for 
A^  > a set  of  lower  bounds  to  the  eigenvalues  of  Eq.  (5.18)  is  obtained. 
This,  so  far,  does  not  require  B to  be  positive  definite.  But  if  the 
eigenvalue  problem  corresponding  to  the  A^P^’s  is  not  solvable  or  not 
easily  solvable,  in  the  case  of  B positive  definite,  a new  sequence  may 
be  defined 


A'(°)  — A'^1^  . . . ^ A'(p)  ^ A'(p+1) 


— A 
• • • 


(5-25) 


where 


(x)  = + B20B2  = A^  + B 

o o 


(5.26) 
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By  theorem  1.9  B 1 — B and  then  Eq.  (5.25)  follows.  The  important  thing 
in  this  connection  is  that  the  inner  projection  is  entirely  arbitrary. 
Consider  for  instance  the  Bazley  projection  defined  in  Eq.  (I.57): 


B'  = | Gi)  (th  | B _ 1 J Ch  )<jh|  . 


The 


operator 


is  now: 


(5.27) 


.(p) 


k=o 


+ a 


P+1 


(i-£  k>kd 


k=o 


(5-28) 


+ |lh)<ffi|  B-1|lh><lh|  . 

The  solutions  of  A^^uj^  = a^^u^^  will  be  spanned  by  the  functions 

u°,  k = o,...,p  and  by  the  linear  manifold  Dn  . We  may  choose  it  to 
consist  precisely  of  the  same  functions  u°.  Then  as  Bazley  and  Fox 

[34] 

have  shown  , the  solutions  may  be  shown  by  solving  a secular  problem. 

In  factj  the  common  feature  of  the  intermediate  problems  is  that  they 

lead  to  a finite  secular  equation  for  the  lower  bounds,  * ^ . 

We  want  to  emphasize  that  partitioning  technique  can  give 

explicit  solutions  to  the  secular  equations  of  the  intermediate  problems, 

the  only  requirement  being  that  equivalent  truncations  are  made  for 

or  A1^^  and  the  appropriate  reduced  resolvents.  Consider  for 

instance,  Miller 1 s work  on  the  excited  states  of  He.  He  applied 

Bazley's  Eq.  (5.23)  to  Hq.  More  explicitly,  given  Hq  + V for  the  He 

atom,  where  V is  just  r^,  t'le  sPectra^  resolution  of  H^  is  known. 

The  truncation  H^^  leads  then  to  a secular  equation.  Instead  of 

deriving  the  secular  equation  in  the  same  way  as  Miller,  we  point  out 

that  the  truncation  T (p)  is  just  the  inverse  of  H in  the 

o' r 7 J o 
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orthogonal  complement  to  \|/°.  Therefore  it  follows  that,  the  solution 
of  Eq.  (5-17),  given  by  the  system  of  Eq.  (5.16),  for  £ = £ should 

-L  IV 

coincide  with  one  of  the  roots  of  Miller's  secular  equation,  just 

because  P(£  - )Tq( p)  = P.  Equation  (5.16)  is  one  solution  of 

HqP)  + V,  when  £q>k  =6..  We  can  verify  this  explicitly,  and  this 
constitutes  an  alternate  proof.  Of  course  the  proof  is  similar  to  the 
transformation  done  in  Chapter  III. 

Consider  Eq.  (5.16).  If  we  replace  d^(k)  by  d^,  defined  by 


o 1 1 E°  - E° 

di  ■ «£-e;+i>  - (£-EP  ] 211 1 


<£-Ep+l>(£-Ep 


(5.29) 


The  only  term  that  changes  in  Eq.  (5.I5)  is  the  kth  term  in  the 
summation: 

v d (k)C  = v d.  C.  . + v . (6  -E.°)  "LC.  , 
mk  kk  mk  k kk  mkv  k'  kk 

But  if  we  impose  - E°  = C , Eq.  (5.16)  may  be  written: 


° Yj  + 6r^)Cik 

i= o 


(5.30) 


The  condition  for  existence  of  a non-trivial  solution  is 


IVA  + Sn^l  ‘ 0 


(5-31) 


Note  that: 


mi 


V( l - E°  ) 

<*°l  „p+1  Iff 

l - E°  -V  * 

P+1 


m1 


(£-Ep+i)8rf  • <0(Ep+r£>v)'1lff 


(5-32) 
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Then  Eq.  (5.31)  can  be  written  in  the  following  way: 


(1-  T^-f)  - — ~Vi)(VVi 


f 1 pO'  'vml''“p+l 

C - 


th 


& - EJ+1> 


= 0. 


If  we  multiply  the  !b  column  by 


o'  - EC 

c- 


(£'-bp°+1)(e»-e“+1) 


we  obtain 


mj? 

E°-E° 

& P+1 


• <<l<Ep°+1-WlK> 


= 0 


(5-33) 


which  is  indeed  the  secular  determinant  solved  by  Miller. 

We  haven't  derived  Eq.  (5.33)  because  it  is  in  any  sense  a 
better  equation  to  work  with.  We  wanted  to  establish  a formal 
connection  between  the  partitioning  formula  in  Eq.  (5.17)  and  the 
solutions  of  the  intermediate  problem  given  by  Eq.  (5.33).  The  fact  is 
that  Eq.  (5.I7)  provides  an  iterative  method  of  solution  to  Eq.  (5.33), 
and  any  pair  bracket  at  least  one  eigenvalue  of  the  intermed- 

iate problem. 

It  would  seem  then  that  the  partitioning  method  does  not  offer 
anything  but  a convenient  way  of  solving  the  intermediate  problem. 

For  the  ground  state,  provided  £ is  an  upper  bound  to  Eq,  the  bracketing 
function  yielded  a lower  bound  which  was  worse  than  the  ground  state  of 
the  intermediate  problem.  For  the  excited  states  we  were  not  able  to 
verify  Eq.  (5.8),  and  hence  it  seems  that  the  only  useful  information 
comes  from  the  intermediate  problems.  The  difficulty  here  lies  in  the 
use  of  a single  reference  function  as  we  will  see  in  the  next  section. 
Furthermore,  the  truncation  of  T (p)  is  by  no  means  the  only  way  of 
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evaluating  the  bracketing  function  Eq.  (1.50),  as  shown  in  PTX,  PTXI, 

r7n 

and  in  the  book  by  ChoiLI  J on  He-like  ions. 


^•3.  The  Bracketing  Function  Associated  with 
Consider 


6U  - K + <<l(V-1-T0(pJ71|t°> 

with  T (p)  defined  in  the  orthogonal  complement  to  \j/°.  Clearly, 


(5-34) 


d£ ' 


d£ 


^ i 0 


(5-35) 


Then  Eq.  (5-34)  is  a bracketing  function.  The  asymptotes  of  Eq.  (5-34) 
are  different  from  the  eigenvalues  of  PHP,  which  are  the  asymptotes  of 


l,k 


= Ek  + 


<t£l(v 


-1 


-1 


10 


(5-36) 


Butjfrom  the  previous  section  we  may  conclude  that  the  asymp- 
totes of  Eq.  (5.34)  are  given  by  the  eigenvalues  of 


PH^P  = P(H^p)+v)P^  PHP  . 

By  theorem  1.2  we  have:  E^P^  ^ E^,  in  order.  So  if  we  use  the 

ordering  theorem,  and  examine  under  what  conditions  Eq.  (5-34)  blows 
up,  lower  bounds  to  the  E^’s  may  be  found. 

There  are  two  consequences  of  this  fact.  One,  that  these  lower 
bounds  may  be  used  in  conjunction  with  the  Eckart-type  formulae  of 
Chapter  II,  and  with  the  Temple  formula,  as  mentioned  in  the  discussion 
following  Eq.  (2.I5). 
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The  second  consequence  is  in  connection  with  the  inequality 
c(p)  - which  holds  under  restricted  conditions.  From  the  preceding 
i-  foi-ows  that  Eq.  (5-35)  holds  and  that  the  asymptotes  of  Eq.  (5.34) 
are  lower  bounds  to  t for  all  £.  . We  can  examine  the  situation  in 
a diagram.  See  figure  6.  It  follows  that  there  are  ranges  where 
■-vP)  — may  have  several  situations,  and  for  the  purpose  of 

w.iis  discussion  it  is  enough  to  consider  a well  behaved  case,  where 
we  have 


Eo  < Eo<  Ei  < E1<  Ei  ^ E1  * E2  • 

It  is  clear  that  in  the  intervals  (E^Eh),  (E^,E  ),  t(p)  is  larger 

— an  t instead  of  being  a lower  bound  to  it.  But  regardless  of  the 

complexity  oi  the  situation,  we  may  find  the  asymptotes,  and  hence  the 

eigenvalues  of  PH'P,  which  are  lower  bounds  to  the  E.’s  of  Chapter  II. 

Consider  Eq.  (5.14),  the  explicit  formula  for  (\jr°  | t(p)  j\J/°) . It  blows 

c,  whenever  ( 1+  .4  ) SL-o,  1,  . . . , p,  blows  up.  Consider  for  simplicity 

she  case  p=o.  Then  it  blows  up  for  (l+A)  = o,  or- 

oo 


1 + v d =0 
oo  o 


(5-37) 


which  may  be  written,  using  Eqs.  (5-32)  and  (5.12) 


1 - t(£-E°)  - <+°!(E°-£  + v)-1|„°>]  -t_  = 0 

L - Ei 

and  then, 

(t°|(E°-£+v)-1k°)  “ 0 (5.38) 

■S  the  simplest  formula  that  gives  a lower  bound  to  E^  within  this 
scheme.  Or  course  Eq.  (5.14)  provides  a whole  family  of  equations  in  £ 
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that  when  solved,  give  lower  bounds  to  the  E^’s;  Eq.  (5.38)  does  not 
give  a very  good  lower  bound,  but  it  may  be  improved  by  considering 
higher  p’s. 

We  see  then  that  a byproduct  of  using  the  bracketing  function 
to  compute  the  eigenvalues  of  the  intermediate  problem  is  that  it 
provides  simultaneously  a scheme  to  compute  lower  bounds  to  the  eigen- 
values of  PHP  and  hence  quantities  useful  in  the  estimates  of  the 
Weinstein  function,  as  in  Eqs.(l.29),  (2.10),  and,  Temple's  formula 
Eq.  (2.9)-  These  estimates  may  prove  very  useful  in  connection  with 
the  full  resolvent  technique(Section  1.2).  The  integrals  appearing  in 
the  scheme  have  been  evaluated  numerically  by  Miller,  but  further 
research  should  prove  interesting,  especially  to  obtain  them  in  closed 
form  and  study  their  properties. 

The  drawoack  with  the  intermediate  problems  associated  with  a 

truncation  of  H , and  hence  with  the  truncation  T (p)  is  that  the 
o o' 

discrete  spectrum  is  by  no  means  complete.  To  expect  convergence  to 
the  exact  eigenvalue  one  should  be  able  to  compute  the  full  spectral 
resolutions.  The  integration  over  the  continuum  is  no  easy  task.  In 
the  next  chapter  we  discuss  this  further  and  make  an  attempt  to  cope 
with  this  difficulty.  Far  more  important  is  the  fact  that  the  trun- 
cation T^(p)  is  not  the  only  way  to  evaluate  the  bracketing  function. 

As  LBwdin  showed  in  PTX  and  PTXI,  the  inner  projection  on  V,  or  on  t, 
when  by  suitable  conditions  its  possitivity  is  insured,  yields  expres- 
sions where  the  operator  T^  has  been  transformed  away.  We  refer  to 
the  substitution  th  = (£  - H ) j in  Eq.  (I.57).  This  implies  the  use  of 
the  linear  manifold  jj,  called  LBwdin  Space  by  Choi^^  . This 
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substitution  removes  the  dependence  on  T , and  the  only  requirement  is 

one  on  the  domains  of  the  operators  involved.  Any  qualms  about  the 

accuracy  of  this  "variational"  estimate  of  T are  removed  by  the 

o J 

accuracy  of  Choi's  results,  and  they  go  beyond  the  intermediate  problems. 

The  multidimensional  reference  function,  discussed  in  [12]  and 
PTXI  has  the  further  flexibility  of  allowing  a treatment  without  any 
reference  to  an  intermediate  Hamiltonian,  i.e.,  provided  that 
i)  £ ^ E if)  the  reference  function  is  a linear  combination  of 

the  first  p+1  eigenfunctions  to  Hq,  iii)  the  Aronszajn  space  g, 

Eq.  (I.56),  contains  at  least  these  same  functions,  then  there  is  an 

1 A 

inner  projection  of  V,V  — V2QV2  associated  with  an  operator  t',  which 
fulfills  t'  < t,  and  then  = Eq  + (i)/-0 1 t ' | \|r°}  and  £.  bracket  one 
eigenvalue  of  H.  This  is  thoroughly  discussed  in  PTXI  but  we  mention 
it  for  completeness. 


5.4.  Improved  Energy  Estimates 

The  starting  point  is  Eq.  (l.l4),  and  by  use  of  Eq.  (I.I5)  and 
Eq.  (1.18)  we  obtain 


E ^ C C ” q ( E+9  €.)  £ " q(  E+96) 


which  may  be  transformed  into: 


E = £ - 


1 + ;q(  E+9  ) 


(5.39) 


We  shall  restrict  our  attention  to  the  gound  state,  and  assume  that  <£, 
is  an  upper  bound,  in  which  case  <5^  is  a lower  bound,  i.e.,  £ E - £. 
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Eq.  (1.17)  tells  us  that 

£ i E ± £*  ± £ 


(5.40) 


where 


l*  = e.  - 


£--£i 

1 + q(£) 


(5-41) 


In  fact  Eq.  (5.40)  follows  from  Eq.  (5.41)  and  Eq.  (5.39)  if  we  make 
use  of  theorem  1.1,  as  q(E+06)  q(  £,)  . 

Use  of  theorem  1.1  in  Eq.  (5.39)  leads  to: 


£- 


4 E 4 £- 


£-  £- 


1 + q(£1)  1 + q(<£) 

and  for  4.  £^,  where  is  a computable  lower  bound  to  the  brack- 

eting function: 


£'±  £-  — £ ~ ^ 

1 1 + q(£i) 


£ - € 


r,  4 E k £ - 


l 


1 + q(£) 


4 £ 


(5-42) 


The  inequalities  in  Eq.  (5-42)  are  very  interesting  because  they  indi- 
cate that  if  we  start  with  an  upper  bound  <f,  an  improved  upper  bound, 
and  an  improved  lower  bound  may  be  obtained  provided  the  slopes  may 
be  evaluated  at  the  points  £ and  £ . As  it  might  be  anticipated 
these  slopes  cannot,  in  general,  be  obtained  exactly,  yet  bounds  to 
them  may  be  used  in  Eq.  (5.I+0). 

If  we  recall  the  definition  of  q(g.)  : 


q(£)  = (^|HT2H|t)  (5.43) 

• * 2 

it  is  apparent  that  operator  inequalities  involving  T are  the  required 
tools  on  account  of  theorem  I.3.  Equation  (2.7)  gives  a bound  for  T, 
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provided  £ < E^: 


0 p T > 


As  T and 


7 — =-  commute,  it  follows  that 

C “ -J^ 


(5-44) 


Alternatively,  Eq.  (5*44)  may  be  derived  from  the  spectral  resolution 


This  is  a very  interesting  relation,  because  it  gives  a bound  to  the 
slope  of  the  bracketing  function  in  terms  of  the  mean  square  deviation 
of  the  Hamiltonian  with  respect  to  the  bracketing  function  and  of  the 
first  eigenvalue  of  PHP. 

Lower  bounds  to  E1  may  be  replaced  into  Eq.  (5.45)  without 
violating  the  inequality,  and  we  have  then: 


that  were  suggested  in  the  previous  section.  Substitution  of  Eq.  (5.45) 
in  Eq.  (5.42)  yields: 


2 

of  T . Substituting  into  Eq.  (5.43): 


q(£)  £ (^iHPHU)  = fell)2 

(Er<£)2  (Er£.)2 


(5.45) 


(5.46) 


Thus  we  see  a further  application  of  the  lower  bounds  to  E 


1 


£-  4 


1 + (AH)2/^-^2 


E < - 


(5.47) 


aaSSxq  XxqqSxgs  aq  pqnoqs  3 put?  xi^s  sx  (H^)  ‘suoxqoung  3ab«  aqBxnooB 

o 

XxqBuosBax  aog*asnBoaq  Suxqsaqaqux  AxaA  uiaas  qou  saop  (6tf'£)  uoxqBnbs 


• g(HV)  < g(<H>-T3) 

qsqq  papxAoqd  <H>  puB  ^3  go  aSsqaAB  aqq  usqq  xaggBras  sx  3 

(oC-Q 

<h>  > a > ^a  pub  £<h>  > a >Ia 

(6rO 

:uaqq  aABq  a ft 

(<H>-Ta)/g(H^)  + <H>  - Ta 
3 ^HV)  <H7  - 

s(<H>-T3)/g(HV)  + I 

X <h>  — a 

(<H>-3)/s(HT7)  + <H>  - <H> 

:uaqi 

*(8'tT*^)  *^3  UT  ~^3  -*°J  cI,a  ‘punoq  agduiag  aqq  asn  ‘aouBqsux 

qo3 

•XSxaua  aqq  SuxqaqoBqq  gBAxaqux  aqq  uxqqxw  qnq  ‘punoq  jawog 

b ao  aaddn  ub  XxxxBSsaoau  qou  ‘ASxaua  aqBuixxo:tddB  ub  SuxuxBqqo  snqq 
punoq  aswox  b aoEgdax  sa  ‘sX^xtbab  qou  sx  sxqq  sy  *3  uBqq  qaxxuuis  sx 
qoxqw  ^3  oq  punoq  xaddn  ub  paau  3M  Suxqsads  Xxqoxaqs  ‘(8^’^)  *b3  ux 
3°  sqBuixqsa  ub  aoBxdax  03  sx  axaq  quaradogBAap  go  auxg  agqxssod  y 


(W'S) 

?(<H>-T3)/  (H7)  + I 

<H7  7 ^ -SL <H?  > 3 

3~  <H> 

moqg  (9*3) 

uxBqqo  aw  ‘(6-S)  *t>a 
’ ba  ux  Bgnuuog  Sjagduiag  paAxqap  aa  uaqw  sb  Xbw  buibs  aqq  ux 

‘ <h>-  Ta  = 3 -Ta 

qas  aa  ‘^3  oq  spunoq  xaMog  SuxpxAOxd  go  puaqsux  gg 


Lo  I 
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than  the  Temple  bound  E^.  Therefore  we  may  expect  that  the  jj-h»s*  of 
Eq.  (5.^2)  should  enable  us  to  obtain  more  insight.  Let  be  any 

computed  lower  bound  to  the  bracketing  function: 


f - £' 

£ p . 1 

1 1 + q(£j) 


< E < £. 


(5-51) 


In  order  to  apply  Eq.  (5-5 l)  we  need  a lower  bound  to  q^, 
different  from  the  trivial  one,  equal  to  zero.  Here  we  may  appeal  to 
the  inner  projection  concept.  By  theorem  1.9,  as  T2  £_0. 


T2  > (-T)Q(-T)  > 0, 

where  Q is  a projection  operator.  Hence,  by  theorem  1.3  we  have 

HT2H  > h(-t)q(-t)h  > 0. 


Choose  Q - 0^  the  first  projector  in  the  spectral  resolution  of  (-T). 
Then  we  have: 


HT2H  > 


H01H 

(V£)2 


(5-52) 


We  may  derive  an  Eckart-type  theorem  for  0 


r 


PHP  - E 


ip'E(Ert)°KJ(VEi)(p-°i)?0  • 


k=2 


Then, 


PHP  - EjP 

e2  - if 


> P - 0, 


or 


°1> 


P(E  -H)P 
E'2  - Eq 


(5-53) 
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Replacing  Eq.  (5*53)  in  Eq.  (5.52)  we  obtain  for  q(£): 

<»|g(VH>PH|*>  . (5_; 

Equation  (5*5 ^0  should  be  studied  in  conjunction  with  Eq.  (5.5I). 

We  see  that  estimates  of  the  eigenvalues  of  PHP  appear  as  well 
as  the  expectation  value  of  H^,  which  seems  a rather  gloomy  perspective 
but  for  reference  functions  constructed  in  terms  of  eigenfunctions  of 
Hq  the  questions  of  operator  domains  may  be  coped  with.  This  seems  to 
indicate  that  the  energy  estimates  and  improved  bounds  through  compu- 
tation of  derivatives  of  the  bracketing  function  offer  possibilities 
of  development . and  application. 


CHAPTER  VI 


AN  EIGENVALUE  TRANSFORMATION 

In  the  preceding  chapters  we  have  studied  some  aspects  of  par- 
titioning technique  and  the  calculation  of  upper  and  lower  bounds  to 
eigenvalues.  An  all  pervading  feature  was  the  need  for  a splitting  of 
the  Hamiltonian  into  an  unperturbed  part,  H^,  and  a perturbation,  V. 

The  spectrum  of  the  unperturbed  Hamiltonian  or  at  least  its  lower  part, 
i.e.,  the  first  eigenvalues  and  eigenfunctions  had  to  be  known,  in 
order  to  construct  a tractable  reference  function.  The  spectrum  is 

not  actually  needed  to  compute  the  reduced  resolvent  T if  one  carries 

o 

the  transformation  to  the  LBwdin  space,  as  in  PTX  and  in  reference  [71] 
Of  course,  the  first  part  of  the  spectrum  is  needed  in  the  method  of 
truncations.  At  any  rate,  the  need  for  an  H^  with  a partially  discrete 
and  known  spectrum  is  not  necessarily  compatible  with  a positive  V 
prerequisite  for  inner  projections  and  all  the  transformations  employed 
previously. 

In  atoms  it  is  natural  to  split  the  Hamiltonian  into  an  unper- 
turbed part  with  the  desired  characteristics  and  a perturbation  V, 
which  is  positive  definite.  In  molecules,  if  one  considers  an  unper- 
turbed part  for  which  the  lower  part  of  the  spectrum  is  known,  the 
perturbation  is  not  positive.  If  one  chooses  the  perturbation  in  the 
same  way  as  in  atoms,  i.e.,  the  interelectronic  repulsion,  the  unper- 
turbed spectrum  has  no  known  spectrum.  This  dilemma  was  faced  by 
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r 72] 

Bazley  and  Fox  ‘ , who  suggested  a method  of  truncations  combined  with 

a splitting  of  the  unperturbed  Hamiltonian  into  parts  corresponding  to 
the  different  nuclei.  Johnson  and  Coulson^^  carried  out  a calcu- 
lation for  with  this  scheme,  and  the  results  were  not  very  good, 
the  reason  being  that  the  truncations  involved  threw  away  the  signifi- 
cant continuum  contribution. 

The  question  is  then,  is  it  possible  to  transform  the  original 
SchrBdinger  equation  in  such  a way  that  one  has  a complete  discrete 
spectrum?  Is  this  in  any  way  related  to  the  possibility  of  splitting 
the  Hamiltonian  into  two  parts,  one  of  which  is  positive  leaving  a 
tractable  HQ?  The  answer  to  the  first  problem  is  partially  affirmative. 
For  one  electron  problems,  it  is  possible  to  do  an  eigenvalue  trans- 
formation that  yields  a complete,  discrete  spectrum.  This  was  known 
to  SchrBdinger  1 ^ and  to  ..others.  ^5]  Lundqvist^^  discussed  the 
problem  in  terms  of  a Green's  function  approach  and  mentioned  that  for 
many  electron  problems  the  solutions  should  be  interesting, because, 
for  the  unperturbed  problem,  each  "state"  is  described  by  a product 
wave  function  in  which  each  "orbital"  or  one-electron  function  will 
depend  explicitly  on  the  quantum  numbers  of  all  the  other  electrons. 
Further  he  concluded  that  a purely  discrete  spectrum  is  not  obtained 
for  the  many  electron  problems. 

We  shall  discuss  an  eigenvalue  transformation  that  yields  an 
Hermitian  operator,  which  is  positive  definite.  It  includes  as  par- 
ticular cases,  the  transformations  of  Hylleraas  and  Lundqvist.  It  is 
closely  related  to  the  variation  iteration  schemes  discussed  in  Morse 
and  Feshbach^^  and  the  Sturmian  functions  of  Rotenberg^^  which  are 
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identical  with  the  functions  discussed  by  Hylleraas.  The  fact  that 
we  obtain  a positive  Hermitian  operator  implies  that  we  can  use  the 
techniques  of  the  preceding  chapters,  i.e.,  partitioning,  inner  projec- 
tions to  our  advantage.  The  results  will  be  applicable  to  positive 
or  negative  perturbations  and  in  general  to  any  perturbation  that  is 
decomposable  into  a positive  and  a negative  part,  for  bounded  H. 

-L: Transformation  of  the  SchrBdinger  Equation 

We  can  write  the  SchrBdinger  equation  for  an  atomic  system  in 
the  following  way: 

<We-E,X-2Vn  (6.1) 

where  T is  the  electronic  kinetic  energy  operator;  V , the  interelec- 

N e 

tronic  repulsion  operator,  and  VN  = £ r~  . The  ^ are  square  inte- 
grable  and  for  a given  nuclear  charge  Z there  is  a set  of  eigenvalues 
En  with  corresponding  eigenfunctions  We  shall  study  later  the  case 

of  molecules. 

Consider  Eq.  (6.1),  and  instead  of  Z being  given,  assume  that 
an  energy  parameter  £ is  at  our  disposal,  and  that  the  parameter  X is 
to  be  determined  for  a given  i: 

(T+Ve-Z)y  =x(6)VnT  . (6.2) 

For  £ equal  to  any  of  the  E in  Eq.  (6.1),  i.e.,  £ = E . 

n n 

% = yn»  and  Eq-  (6.2)  must  certainly  have  at  least  one  solution 

^(En)  ~ Z.  Equation  (6.2)  is  a typical  Sturm-Liouville  problem,  and 

it  should  be  noticed  that  Z is  the  common  eigenvalue  of  a family  of 

operators  L(E  ) with  corresponding  eigenfunctions  Y . 

11  n 


In  this  sense, 
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the  ^ are  not  degenerate,  because  they  are  eigenfunctions  to  different 
operators . 

Even  though  Eq.  (6.2)  is  an  equation  with  which  we  may  work, 
the  following  transformation  provides  some  insight  into  its  properties. 
Let 


% = V*  • (6.3) 

Then  Eq.  (6.2)  may  be  written 

L(£)0e  = v^(t+V2)vn"\  = 4(£)0£  • (6.4) 

Equations  (6.2)  and  (6.4)  have  the  same  eigenvalues  k(£). 

The  eigenfunctions  are  related  by  Eq.  (6.3).  Notice  that  L(£)  is  an 
Hermitian  operator,  therefore  if  we  drop  the  £ dependence  from  the 
notation,  Eq.  (6.4)  may  be  written: 


L0  = 1 0 
n n n 


(6.5) 


and  from  the  properties  of  Hermitian  operators,  L = L , it  follows 
that 


<0  |l|0  ) = \ (0  |0  >6 

m n n n 1 n nm 


(0  |0  ) = (0  |0  )5 

• ' m 1 n n 1 n'  nm 


(6.6) 


From  the  last  equation  and  Eq.  (6.3)  it  follows  that: 

KK'1  - (’•'JW  ■ <'*'JvN!V8m,  (fi-7) 


which  expresses  the  well-known  fact  that  the  eigenfunctions  of  the 
S turm-Liouville  problem  in  Eq.  (6.2)  are  orthogonal  with  respect  to  the 


weight  factor  = 


.>  0. 


For  £ £ 0,  as  T > 0 and  V 0 and  on  account  of  theorem  1.3 

e J 


we  have 


U0 


(6.8) 


which  is  a fundamental  property  in  what  follows.  Therefore  L is 
bounded  from  below,  and  if  it  has  a partially  discrete  spectrum,  the 
eigenvalues  are  positive. 

Before  discussing  the  spectrum  of  L,  we  will  relate  the  eigen- 
values of  H and  those  of  L,  through  an  application  of  the  Hellmann^^ 
and  Feynman^  ^ theorem  and  indicate  how  this  leads  to  upper  and  lower 

bounds  to  the  E . 

n 


Consider  an  eigenvalue  problem  of  the  form 


A = \BW 


(6.9) 


Then 


<!F|A-XB|y)  = 0 . (6.10) 

If  the  operators  depend  on  a parameter  £ we  may  differentiate 


Eq.  (6.10): 


(H|a-\b|¥)  + <y|A-x.B|H)  + <y|H|y> 


3Yv 


|9A  | 


ax  <y|B|y)  - \<y|§§|if> 


3£ 


= 0 


If 


3B 

d£. 


= 0,  we  obtain 


|A  = /wlIA 


= (vf|gT|T)/('F|B|i') 


(6.11) 


which  is  the  form  of  the  He llmann -Feynman  theorem  that  applies  in  our 
case.  Then. 

= . My)  < n 

di  (y|vN|y)  0 • 


(6.12) 
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This  is  not  at  all  surprising  because  the  function  £.  = £(x) 

has  negative  slope,  i.e.,  to  increase  the  nuclear  charge  decreases  the 

energy  of  an  atom.  The  function  X = X(<f)  is  related  to  the  function 

£ = £(x),  and  in  fact  in  certain  problemSjOne  can  be  obtained  from 

the  other.  In  general,  it  is  unwieldy  to  obtain  £.  = £.(x);  it  is 

the  problem  of  perturbation  theory;  but  in  certain  cases  it  is  possible 

to  obtain  the  dependence  of  the  interaction  parameter  X on  the  energy  £, 

and  this  is  the  basis  for  the  variation  iteration  methods  of  Chapter  9 

of  Morse  and  Feschbach.  Instead  of  discussing  their  approach  we  will 

relate  upper  and  lower  bounds  to  X(£.)  with  upper  and  lower  bounds  to 

the  E . For  a given  £ there  is  a series  of  solutions  X . The  lowest 
n n 

one  is  associated  with  the  ground  state  of  the  system  and  so  on.  If 
there  is  any  doubt  about  this,  it  should  be  realized  that  the  curves 
X(£)  and  £(  A,)  are  the  same,  and  for  a given  X,  there  are  many  values 
of  £ possible,  one  for  each  curve.  The  lowest  one  corresponding  to 
the  ground  state  curve.  In  fact,  £-,(-0  ■'-s  t'rie  isoelectronic  curve 
for  the  ground  state  as  a function  of  the  nuclear  charge  X.  See 
figure  7- 

Suppose  now  that  we  are  able,  by  any  method,  to  provide  upper 
and  lower  bounds  to  X(£),  i.e.,  that  for  any  £ we  may  obtain  X>  (£•) 
and  X^'  (£.),  such  that 

(e)  £ x(£)  < X>  (£)  . (6.13) 

We  should  stress  that  we  are  considering  a specific  state,  not 
necessarily  the  ground  state.  The  important  consequence  of  Eq.  (6.12) 
is  the  following: 


li  6 


Theorem  6.1: 

If  Z is  the  nuclear  charge  of  the  system,  and 

\<.(£)  £*•(£)  £ X^e) 

is  fulfilled  for  some  state,  then  any  £ that  satisfies  X>(£.)  £ Z is 
an  upper  bound  to  the  energy  E of  the  state  under  consideration  and 
any  £.  that  satisfies  X^S.)  2t  Z is  a lower  bound  to  it. 

The  proof  follows  from  the  fact  that  for  X>(£)£.  Z s X(E), 

E — £ because  of  the  negative  slope  of  X.  Similarly  for  the  lower 
bound:  if  X^E)  5'Z  =1X(E)  it  follows  that  £ ^ E.  See  figure  8. 

Therefore  upper  bounds  to  X(£)  are  related  to  upper  bounds  to  E 
and  the  same  applies  to  the  lower  bounds.  We  should  realize  now  that 
this  treatment  is  applicable  to  any  negative  part  of  the  Hamiltonian. 

One  has  to  repeat  the  derivation  that  starts  in  Eq . (6.2),  where 
instead  of  X being  the  nuclear  charge  and  V = E ^r,  X is  an 
artificial  strength  parameter  and  -V  is  a negative  part  of  the  Hamil- 
tonian. The  pertinent  modification  of  theorem  6.1  is  that  Z is 
replaced  by  1,  or  by  the  appropriate  constant  according  to  the  circum- 
stances . 

So  far  there  is  nothing  done  towards  solving  Eq.  (6.1),  because 
Eqs . (6.2)  or  (6.5)  are  just  as  difficult.  But  as  L is  positive, 
actually,  as  any  part  of  L is  positive  by  construction,  it  is  possible 
to  carry  inner  projections  of  L,  or  of  parts  of  L,  which  are  a necessary 
tool  in  obtaining  lower  bounds  to  X(£).  In  addition  it  is  possible  to 
split  L into  an  unperturbed  part,  with  known  spectrum  and  a positive 
perturbation.  In  the  case  of  one-electron  systems,  the  spectrum  is 
discrete.  For  many  electron  systems,  the  spectrum  is  partially  discrete 
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and  partially  continuous.  In  the  next  section,  we  study  the  spectrum 
with  some  detail. 

6.2.  Nature  of  the  Spectrum 

The  spectrum  associated  with  Eq.  (6.2)  and  Eq.  (6.4)  can  be 
studied  by  considering  first  the  unperturbed  problem  associated  to  them, 
namely,  the  one  obtained  by  neglecting  the  interelectronic  repulsion. 
Therefore  we  set: 

L°  " V^(T-<5)V"2  , £<0  . (6.14) 

Both  T and  V are  sum  of  one  electron  operators: 

T - £ T(  i) ; VB  • £ V(i)  . (6.15) 

i i 

Equation  (6.2)  is  then 

^ [T(i)-\°v(i)]Y0  = £y°,  L<  o . (6.16) 

i 

The  problem  is  to  find  solutions  4°  and  Y°  for  a given  £.  As 
antisymmetry  was  a condition  on  the  solutions  of  Eq.  (6.1)  we  shall 
impose  the  same  requirement  on  7°.  We  purposely  will  ignore  all  other 
questions  of  symmetry  of  the  solutions  of  Eq.  (6.l6).  As  the  operator 
involved  is  a sum  of  one  electron  operators,  we  may  seek  solutions  in 
the  form  of  antisymmetized  products  of  spin  orbitals,  i.e.,  construct 
determinants,  taxe  linear  combinations  of  them  in  order  to  construct 
eigenfunctions  of  L ,S  , etc.  But  for  our  purposes  it  is  enough  to 

assume  that  all  that  can  be  done  following  the  treatment  given  in  the 
book  by  Slater 
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Let  A be  the  antisymmetrizer  operator.  We  seek  solutions  of 

the  form 

y0  = aTTu  (j)  . (6.17) 

j J 

Replacing  Eq.  (6.I7)  in  Eq.  (6.l6)  we  obtain 

A^TT  c j ) Ct(  i) -X.°v( i)]U;L(  i)  = 6y°  , (6.18) 

i j^i 

where  we  have  used  the  fact  that  the  antisymmetrizer  commutes  with 
the  symmetric  operators  T and  VN«  Let  us  set  £ = £ £ where  the 
Vs  are  entirely  arbitrary.  Then  Eq.  (6.18)  becomes 

a^TT  uj(j)CT(i)-X.°V(i)-  £.3u.(i)  =0  . (6.19) 

i 

This  allows  us  to  assert  that  Eq.  (6.I7)  is  a solution  of 
Eq.  (6.16)  provided  that  the  following  conditions  are  fulfilled: 

(0  t = Yi  6.  ^ 0 

• jL 

x 

(ii)  [t(  1) ~x?v(  1) - £.]u.(l)  = 0 (6.20) 

(iii)  = X,°  for  all  i 

Consider  condition  (ii).  If  V(l)  = r^\  it  is  entirely  anal- 
ogous to  the  hydrogen  atom  eigenvalue  problem,  but  with  the  energy  as 
a fixed  arbitrary  parameter  and  the  nuclear  charge  as  the  eigenvalue 
to  be  determined. 

If  one  uses  the  method  of  Frobenius  to  solve  (ii)  one  finds 


the  same  terminating  condition  for  the  series  as  in  the  hydrogen 
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atom  case,  namely  that 


> • • • y 


oo 


(6.21) 


The  eigenfunctions  are  hydrogen-like  functions,  but  with  no 
dependence  on  the  quantum  number  n in  the  argument: 


fact  that  they  form  a complete,  discrete,  set  was  used  by  him  and 


mation  of  the  SchrBdinger  equation  for  the  S-states  of  the  hydrogen 
atom  which  yields  essentially  the  same  solutions  as  Eqs.  (6.21)  and 


and  contribution  of  the  continuum,  see  page  27^-  of  reference  [36] . 

In  the^  present  context  the  significant  point  is  that  for  a one-electron 

problem  and  one  center  the  transformation  being  discussed  yields  a 

complete  discrete  spectrum.  Unfortunately  this  does  not  follow  for 

many  electron  problems.  Consider  condition  (i)  of  Eq.  (6. 20) . The 

sum  of  "one-electron  energies"  can  be  negative  even  if  one,  or  at  the 

most,  n-1  of  them  are  positive.  Therefore,  there  are  solutions  of 

Eq.  (6.l6)  which  correspond  to  positive  and  these  form  a continuum. 

At  any  rate  let  us  consider  the  discrete  part  of  the  spectrum  and 

examine  the  solutions  in  more  detail.  We  have  to  impose  condition  (iii) 

Let  the  <£..  * s satisfy 
1 J 


(6.22) 


These  functions  were  already  known  by  SchrBdinger . The 


by  others,  ^-L  There  is  a discussion  by  Hylleraas^2^ 


of  a transfor- 


(6.22),  and  where  he  proves  the  completeness  of  the  discrete  spectrum 
of  ( ii)  for  <£^  < 0.  Clearly,  they  are  orthogonal  with  weight  factor 
V(l)  = r^  . For  a more  detailed  discussion  of  different  basis  sets, 
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d.  = 6P.2,  £<  o.  (6.23) 

Then  Eq.  (6.21)  becomes 

\°  = n .p  . V" (6.24) 

lix  . 

The  p^’s  are  integers  at  our  disposal.  Let's  denote  a set  of  p/s 

by  p = [p/^.  Then,  given  a set  p,  we  impose  a restriction  on  the 

associated  n.’s: 
l 

nipi  = nM(p),  n=l,  2,  . . . , 00  s (6.25) 

where  M(p)  is  the  lowest  common  multiple  of  the  set  p.  In  such  a case 
Eq.  (6.24)  b ecomes 


L°  = L°  =nM(p)V^r,  for  all  i. 
l n,p 


(6.26) 


We  see  then  that  the  set  p determines  a "configuration,"  and  for  each 

configuration  there  is  a set  of  eigenvalues  X°  , and  a set  of  eigen- 

n,  p 

functions,  antisymmetized  products  of  spin  orbitals,  where  the  radial 
functions  are  given  by  Eq.  (6.22)  provided  the  condition  in  Eq.  (6.25) 
is  imposed  on  the  indices  ck.  The  eigenfunctions  are  given  by 

5,0  * C - AYRn.i.».(p ^ ri)Yi.».(ei'*i)fp(i)]  <6-27) 

c 111!  11  l J 


where 

1 pi 

We  should  notice  that  this  shows  explicitly  how  each  one-elec 
tron  function  which  appears  in  the  solutions  of  Eq.  (6.2)  depends  on 
the  quantum  numbers  of  all  the  other  spin  orbitals  as  suggested  by 
Lundqvist.  It  should  be  emphasized  that  not  all  choices  of  the  set  p 
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are  independent.  For  a given  there  is  a restriction  on  any  set  p 


• . 


c ■ eE  pi2  - e'  E pi£ 


(6.28) 


or 


< - eE  E 


(6.29) 


Therefore,  for  any  set  p',  the  eigenvalues  X°  « are 


related  to  the  eigenvalues  X.  by: 

n,  p 


° , = nM( p ' = nM(p)V^2£  wjj: 


n,p 


(6.30) 


Any  choice  p'  such  that  p!^  = api  should  be  discarded  because,  by 
Eq.  (6.30)  does  not  give  anything  new.  If  pi  = 0$^;  M(p')  = QM(p)  and 


A°  , = X°  ■ 
n,p  n,p 


are  : 


For  the  two-electron  problem  the  choices  of  "configurations  p 

(1,1),  (1,2),  (1,3),  (2,3),  (2,5),  (2,7),  •••;  (3,4),  (3,5), 


(3,7),  . ..,  etc.  Each  configuration  will  lead  to  a set  of  solutions 

by  letting  n assume  the  values  1,2,3, •••>“  • 

The  discussion  of  the  full  operator  L should  have  the  preceding 
as  a starting  point,  but  for  our  purposes  it  is  sufficient  to  have 
found  an  unperturbed  operator  with  a partially  discrete  spectrum,  with 
known  solutions  and  therefore  tractable  by  the  methods  of  the  preceding 
chapter.  We  have  not  obtained  a similar  solution  for  the  molecular 
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case,  but  we  shall  take  advantage  of  the  completeness  of  the  one- 
electron  solutions  in  conjunction  with  the  method  of  truncations  in 
a later  section. 

6.3.  Atomic  Systems 


The  fact  that  starting  from  the  SchrBdinger  equation  we  have 
obtained  a positive  definite  Hermitian  operator  has  several  useful 
consequences  that  we  will  discuss  in  this  section. 

First  we  shall  construct  an  inner  projection  of  the  operator  L: 

L = V^)T+Ve-£)V^ 


_!  _! 

= V 2MV  2 
N N 


(6.31) 


For  £,  < 0,  M > 0.  Therefore,  by  theorem  1.9>  we  can  construct 
an  inner  projection  M'  which  fulfills  the  following: 

JL  JL 

M'  = M2QM2  M , (6.32) 

where  Q is  an  arbitrary  projection  operator. 

We  can  define  now 


L ' s VM'V  ’ (6-33) 

which,  by  theorem  1.3,  implies 

L' ^ L . (6.34) 

From  theorem  1.2,  it  follows  that  the  eigenvalues  of  L'  are, 
in  order,  smaller  than  those  of  L,  i.e. 


L *0 ! = X\0\ ; L0. 

1 11  1 


x.0. ; 

1 1 


k!  £ 

1 


X. 

x 


in  order 


(6.35) 
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The  importance  of  this  is  that  we  can,  if  we  wish,  incorporate 
the  SchrBdinger  equation  to  the  host  of  problems  that  can  be  treated 
directly  by  the  method  of  intermediate  problems,  without  appealing  to 
truncations  or  splitting  the  Hamiltonian  into  a perturbed  and  an  unper- 
turbed part.  Thus  we  can  set,  Q = (J£|  > where  iE : is  an  arbi- 

trary linear  manifold,  and  it's  overlap,  or  metric  matrix.  Let  s 
transform  to  the  following  manifold:  M2  1 it)  = M[g),  the  Aronszajn  space. 
Then  we  may  write: 

L'  = v"%|g)^“1(g|MV^  (6.36) 

where  the  metric  matrix  is  given  by: 

^ = ( g| m|  g)  . (6.37) 

With  this  projection,  the  eigenvalue  problem  of  L'  becomes: 

V^M|g)A'l(g|MV^|0'>  * Xr|«’>  (6-38> 

or 

MV^Mlg)^-1  <g|MV^|<$')  = X'MV^|0') 

which  implies: 

[(glMV^Mlg)^"1  - X'  U]  (g|MV^]0))  = 0 . (6.39) 

Therefore,  we  have  non-trivial  solutions  for 
det({g|MV”1M|g)^"1  - X'l  1 = 0 

or  explicitly: 

(<s|(T+Ve-£)V^1(T+Ve-£)|g>  “ X'(s|T+Ve-£!s)}  =0. 

(6.40) 
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The  roots  of  the  secular  determinant  Eq.  (6.40)  are,  in  order,  lower 
bounds  to  the  eigenvalues  of  Eq.  (6.31),  provided  £ < 0. 

Furthermore,  by  our  theorem  6.1,  any  lower  bound  X'  > Z 
obtained  for  a specific  value  of  £ , implies  that  £ is  a lower  bound 
to  the  energy  of  that  state.  In  particular,  the  lowest  root  of 
Eq.  (6.40)  is  associated  with  a lower  bound  £ , to  the  energy  of  the 
ground  state  of  the  system  provided  that  X’(£)  X(Eq)  = Z. 

Another  alternative  is  to  divide  L into  two  parts  L and  L.  , 

o V 


whe 


re 


L ■ Lo  + V l„-  V<T-  >v?' 


.-1. 


L = V 2V  V 2 = V V = V V 
V N e N N e vevN 


-1 


(6.41) 


The  operators  L and  LTT  are  both  positive.  L has  a known 
o V r o 

spectrum,  which,  as  we  saw  in  the  previous  section,  is  partially 
discrete.  It  would  be  interesting  to  compare  the  nature  of  truncations 
of  Lq  with  those  of  Hq.  In  fact,  we  may  apply  any  of  the  methods  dis- 
cussed previously,  perturbation  methods  combined  with  inner  projections, 
partitioning,  etc.  For  the  sake  of  comparison  with  the  intermediate 
problems  we  will  derive  the  secular  equation  corresponding  to  a trun- 
cation of  Lq.  We  might  just  as  well  obtain  the  partitioning  solution, 
but  we  want  to  show  in  one  case  the  characteristic  method  of  arriving 
at  the  secular  equation.  We  will  use  a Bazley  projection  of  L^: 

Ly  > = |(h)^"1(lh|,  zfl  = < Ih  1 VNV~ 1 1 Ch>  . 

The  linear  manifold  Ih  is  at  our  disposal,  and  for  simplicity,  though 
by  no  means  necessary,  we  will  choose  it  as  comprising  p eigenfunctions 
to  L : (h  = §°=  . . . ,0°}  . 
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A natural  truncation  of  L is,  in  this  case, 

o 


p 


L > L* 
o o 


£ WO 


n=l 

Compare  with  Eqs . (0.3),  (6.26)  and  (6.27).  Notice  that  the  trunca- 
tion yields  an  operator  L , smaller  than  L without  including  the 

O o 0 

orthogonal  complement  to  the  subspace  §° . This  is  because  L > o, 

o 

and  it  avoids  integrals  of  the  type  in  Eq.  (5.32). 


We  may  arrange  the  eigenvalues  of  L included  in  the  trunca- 

o 

tion  in  a square  matrix: 


• (6.42) 

Then  we  have 

L > L'  = |#°>*Vl  + \$)Al{t\  . (6.43) 

The  eigenvalue  problem  for  L'  is: 

|$o)(fto+^■1)<$o|0,)  = x<0>  (6.44) 

which  means 


(#°l^°)(^0+^“1  ) T X'  l|  - 0 


h°+^'1  -4'  <^°|^°)"1|  = 0 . (6.45) 

In  terms  of  the  set  Y°  = V 2 0°  , 

n N n 

l&°  +^°|v"1|  y)_1  - x.'1(^°|vNlY°)l  = 0.  (6.46) 
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We  have  to  remember  that  the  condition  on  the  eigenvalues  \\  is 
^i  - *-n  order,  so  if  we  use  theorem  6.1,  any  value  of  £ such  that 

^(6.)  ^ 4^(E)  — Z is  a lower  bound  to  the  ground  state. 

It  is  clear  that  we  can  develop  a partitioning  approach  for 
obtaining  upper  .and  lower  bounds  to  the  eigenvalues  of  L,  starting 
from  the  zero  order  operator  L . We  do  know  its  lowest  eigenvalues 
and  eigenfunctions,  and  the  perturbation  is  positive.  This  is  of  some 
interest  and  warrants  some  future  work,  yet  its  use  does  not  seem  to 
simplify  any  of  the  problems  that,  in  the  a tomic  case,  can  be  solved 
in  the  conventional  way.  On  the  other  hand,  in  the  case  of  negative 
perturbations  it  may  prove  to  be  very  useful. 

6.4.  Molecular  Problems 

To  finish  this  brief  review  of  applications  of  the  transforma- 
tion indicated  in  section  6.1,  we  consider  now  molecular  problems. 

First  of  all,  we  should  point  out  that  the  main  difficulty  in 
treating  molecular  problems  is  the  absence  of  a splitting  of  the 
Hamiltonian  into  an  unperturbed  part  with  partially  known  spectrum  and 
a positive  perturbation.  From  what  was  said  in  the  previous  sections, 
we  should  be  able  to  treat  splittings  involving  negative  perturba- 
tions. 

Bazley  and  Fox1-  suggested  a treatment  for  molecules  that 
seem  to  bypass  the  problem,  but  it  pays  the  price  of  introducing  trun- 
cations and  this,  as  we  say  in  Chapter  V,  restricts  the  accuracy  of  the 
approach  when  the  operator  being  truncated  does  not  have  a complete 
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discrete  spectrum.  It  should  be  clear  by  now  what  we  look  for:  an 
eigenvalue  transformation  that  should  allow  the  use  of  truncations, 
with  hope  of  convergence  to  the  exact  eigenvalues.  We'll  see  how  this 
is  achieved  and  how  truncations  may  be  avoided  completely. 


There  are  several  qualifications  on  the  applicability  of  the  method  in 
terms  of  spectrum,  domains,  limit  points,  etc.,  of  the  Hamiltonians 
involved,  but  we  shall  refer  to  their  paper  for  details.  We  shall 


set  of  low  lying,  discrete  eigenvalues  below  the  first  accumulation 
point  of  H.  For  simplicity  we  shall  consider  a molecular  system  for  a 
given  configuration  of  the  nuclei  and  we  thus  are  treating  the  elec- 
tronic molecular  problem  within  the  Born-Oppenheimer  approximations, 
and  this,  if  solved  exactly  would  lead  to  a lower  bound  to  the  molecular 
energy1  ’ . Since  the  Laplacian  is  invariant  under  translations  of 

the  coordinate  origin  we  may  write  for  the  electronic  Hamiltonian: 


We  shall  first  sketch  the  treatment  due  to  Bazley  and  Fox. 


assume  that  the  spectra  of  H°  and  H,  as  we  will  define  later,  have  a 


0=1  i=l 


(6.47) 


n m 


m 


where 


n 


(6.48) 


a 
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The  d s label  the  n nuclei,  and  the  i’s  label  the  m electrons. 


The  operator  i-s  Laplacian  of  electron  i expressed  with  the 

nucleus  CH  as  original.  The  a^  are  arbitrary  positive  real  numbers 
that  satisfy  Eq.  (6.48). 

We  certainly  have  a splitting  into  an  unperturbed  Hamiltonian 
and  a positive  perturbation,  but  the  unperturbed  part  has  no  known 
spectrum  or  analytic  solutions.  In  fact  for  two  centers  and  one 
electron,  our  zero  order  Hamiltonian  is  just  the  R + Hamiltonian. 

We  may  write  it: 


Even  though  we  do  not  have  analytic  solutions  for  Eq.  (6.49),  we  can 
introduce  truncations  in  the  operators  Hfl  and  H^,  both  with  known, 
hydrogen-like  spectrum. 


(6.49) 


where  we  have  set  a 


1 


la  alb  2 


What  Bazley  and  Fox  have  done  is  suggested  by  Eq.  (6.49). 


The  natural  truncation  will  be  in  this  case: 


P 


P 


(6.50) 


where  E 


we  may  write: 


n m 


m 


(6.5I) 
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where 


E(a)0(a) 


where 


(6.52) 


( a ) = 


and  the  projectors 


are  given  by 


(6.53) 


°ja)  - ><*(“>  I 


(6.54) 


( c0 

where  Yiy  are  just  hydrogen-like'  functions  centered  at  the  nucleus  a. 

It  is  clear  that  this  method  of  truncations,  in  conjunction 
with  the  formulae  of  Chapter  V will  lead  to  lower  bounds  to  the  elec- 
tronic energy.  Further,  if  one  applies  inner  projections  to  V,  and 
chooses  suitable  linear  manifolds,  lower  bound  formulae  with  simple 
integrals  to  compute  are  obtained^^. 

Instead  of  proceeding  further  along  this  line,  and  we  refer  to 
[72] 

Bazley  and  Fox  , we  want  to  point  out  that  this  procedure  cannot 
converge  to  the  exact  results  because  the  truncations  leave  out  the 
whole  continuum  portion  of  H , and  increasing  p will  lead  to  a 
limited  improvement. 

Johnson  and  CoulsonL  have  applied  the  truncations  defined 
by  Eq.  (6.52)  to  the  molecule  Hamiltonian,  Eq . (6.4 9),  and  obtained 
very  poor  convergence.  We  list  their  results  for  g and  u states  as  a 
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function  of  p in  Eq.  (6.52) 


order  of  truncation  p: 

1 

3 

6 

Exact  Ene 

g states(l) 

-1.390 

-1.320 

-1.300 

-1.103 

(2) 

-0.451 

-0.438 

-O.36I 

(3) 

-0.300 

-0.271 

-0.236 

(4) 

-0.214 

-0.178 

u states(l) 

-1.110 

-1.004 

-O.969 

-0 . 668 

(2) 

-0.370 

-0.346 

-0.255 

(3) 

-0.224 

-0.193 

-0.137 

(4) 

-0.136 

-0.127 

The  exact  energies  are 

taken  from 

Bates,  et 

The  method  is  appealingly  simple,  i.e. , one  has  to  evaluate  a 
set  of  overlap  integrals  only,  because  the  operator  involved  is  just 

But  the  price  of  neglecting  the  continuum  in  the  trunca- 
tions. seems  to  be  too  high. 

What  we  should  notice  is  that  we  may  write  the  Schrbdinger 
equation  for  the  H^+  molecule  in  the  form 

((-  - 4 - E)  + <-  5Vlb  - 4 - E>]  * ■ <6-55) 


Now  instead  of  considering  E as  the  eigenvalue,  we  replace  it  by  an 
energy  variable  £,  at  our  disposal  and  replace  the  "nuclear  charge"  2 
by  a parameter  \ to  be  determined.  We  have  then 


v 2 la 


X,  C \ / 1 __  2 

-E)  + (-  2Vlb 


(6.56) 


But  each  term  in  parenthesis  is  exactly  of  the  same  nature  as  the 
operator  in  the  second  equation  of  Eq.  (6.20).  It  has  a complete, 
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discrete  spectrum.  We  may  apply  truncations  and  expect  to  converge 
faster  than  before  to  the  exact  eigenvalue,  determined  by  theorem  6.1. 

Calculations  of  this  sort  were  reported  by  Walmsley^^  and 
the  results  were  the  following: 

p=l  p = 3 p = 6 exact 

-1.267  -1.110  -1.119  -1.103 


Only  the  ground  state  was  treated  and  the  results  should  be 
compared  with  the  ones  obtained  by  Coulson  and  Johnson.  We  see  then, 
that  this  particular  case  of  the  eigenvalue  transformation  in  Eq.  (6.2) 
applied  to  molecules  shows  fast  convergence,  and  this  is  very  prom- 
ising, especially  since  the  truncation  method,  following  the  Bazley 
and  Fox  idea  is  not  the  only  way  to  proceed.  In  fact,  the  only 
requirement  is  that  theorem  6.1  should  be  implemented  in  some  way  or 
another.  So  we  may  write  the  eigenvalue  problem  for  the  molecule  in 
the  form: 


(6.57) 


We  carry  out  an  eigenvalue  transformation  of  this  equation,  and  replace 
E by  an  arbitrary  energy  £. , and  we  solve  for  an  eigenvalue  parameter 
that  multiplies  the  right  hand  side  of  Eq.  (6. 57)*  Theorem  6.1  is 
implemented  by  imposing 


MO  > ^(E)  = 1. 
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and  we  can  treat  the  modified  version  of  Eq.  (6.57)  in  several  ways. 
Further,  we  may  separate  any  part  of  the  nuclear  attraction  as  the 
negative  part  of  the  Hamiltonian.  An  interesting  form  is  given  by 
the  following,  where  for  simplicity  we  consider  the  E^  molecule: 

{("H-  ^ + ("7  a-  ^ + 4 1 ■ + 


or 


(<  v t> + ( vf> + f ‘ x(v + ?r)f 


(6.58) 


Each  of  these  splittings  has  to  be  explored  further.  At  the  moment 
we  may  say  that  they  may  be  treated  by  a wide  assortment  of  possible 
methods,  some  similar  to  the  one  derived  in  the  previous  section, 
which  was  obtained  by  an  inner  projection  of  the  positive  part  of  the 
Hamiltonian.  One  may  apply  perturbation  methods,  either  combined  with 
truncations,  with  results  similar  to  those  of  Chapter  V and  of  Eq.  (6.56). 
For  ' H^  , this  is  simple  and  successful.  If  the  difficult  integrals 
that  appear  when  one  uses  truncations  are  to  be  avoided  one  should 
carry  an  inner  projection  of  the  interelectronic  repulsion  and  obtain, 
for  conveniently  selected  manifolds,  simple  secular  equations . 

What  is  the  role  of  partitioning  in  this  case?  Twofold*  it 

provides  explicit  solutions  to  the  secular  problems  that  arise  from 

truncations  and  more  important,  it  can  bypass  the  intermediate  problems 

by  choosing  convenient  manifolds  in  the  inner  projections.  We  only 

need  an  unperturbed  Hamiltonian,  with  a partially  discrete  spectrum, 

and  a perturbation  which  must  be  just  a sum  of  a positive  and  a 
negative  part. 


The  positive  part  may  be  treated  by.  the  conventional  inner 
projection  technique,  and  the  negative  part  by  applying  the  method 
suggested  in  this  chapter.  This  is,  in  fact,  a starting  point  for 


future  work. 


CHAPTER  VII 


DISCUSSION 

In  the  preceding  chapters  we  have  made  a sketchy  survey  of  the 
present  status  of  the  partitioning  technique  and  lower  bound  theory. 
Actually  we  have  almost  neglected  some  of  the  most  important  aspects, 
like  multidimensional  partitioning  and  we  have  focused  the  attention 
on  some  specific  problems.  This  is  rather  preliminary  in  the  sense 
that  further  work,  both  formal  and  numerical,  should  be  done,  yet  it 
seems  worthy  to  review  the  results  obtained  so  far  because  they  hint 
at  new  problems  and  applications. 

Perhaps,  we  should  stress  that  lower  bound  calculations  have 
not  only  a formal  value,  but  they  may  settle  questions  with  direct 
physical  significance.  For  instance,  some  doubly  excited  states  of  Li 
have  been  spectroscopically  identified  with  the  help  of  simple  lower 

f Qy-j 

bound  criteria  j lower  bound  to  the  repulsion  energy  of  two  He  atoms 
in  the  short  distance  range  would  be  very  helpful,  if  accurate  enough, 
to  discuss  the  effects  of  non-adiabaticity . 

Properties  other  than  energy,  calculated  as  expectation  values 
over  variational  functions  or  by  low  order  perturbation  theory  span  a 
field  where  the  theory  of  upper  and  lower  bounds  should  be  applied  in 
order  to  have  a reasonable  confidence  in  the  significance  of  the  com- 
parison of  calculated  and  experimental  properties.  It  is  well-known 
that  improved  variational  functions  do  not,  by  any  means,  insure  better 


134 


135 


fOO  OqI 

expectation  values  ’ * It  is  hoped  that  further  study  of  operator 

inequalities,  of  the  type  indicated  in  Chapter  I will  be  useful.  In 
fact  some  of  the  conditions  being  used  are  too  strong,  i.e.,  instead  of 
positive  operators  we  need  only  to  consider  operators  in  their  domains 
of  positivity,  etc. 

In  Chapter  II,  we  showed  that  the  well-known  formula  for  lower 
bounds  to  the  ground  state  due  to  Temple  arises  from  a truncation  of 
the  spectral  resolution  of  the  reduced  resolvent  operator  T,  funda- 
mental tool  of  the  partitioning  approach.  This  indicates  under  what 
conditions  the  partitioning  approach  may  yield  a better  bound  and  a 
way  of  removing  one  of  the  drawbacks  of  Temple's  formula,  namely,  the 
requisite  of  knowing  the  first  excited  eigenvalue  of  H.  The  partitioning 
approach  shows  that  one  actually  needs  a lower  bound  to  the  ground 
state  of  PHP,  where  P is  the  projection  operator  of  the  orthogonal 
complement  to  i|r,  the  reference  function  under  consideration.  Lower 
bounds  to  the  eigenvalues  of  PHP  were  shown  to  be  useful  in  estimating 
the  Weinstein  function  Eq.  (1.26),  as  in  Eq.  (1.29)  and  Eq.  (2.10). 

The  Weinstein  function  is  the  key  quantity  in  the  full  resolvent  tech- 
nique and  recent  work  by  Sasakawa^^  using  Green's  functions  and 
dispersion  relations  may  be  implemented  through  the  explicit  calcula- 
tion of  the  Weinstein  function.  In  Chapter  V we  have  obtained,  by 
using  the  partitioning  technique  to  solve  the  intermediate  problems 
arising  from  truncations,  formulae  that  give  the  desired  lower  bounds 
to  the  eigenvalues  of  PHP.  A problem  related  to  obtaining  bounds  to 
the  eigenvalues  of  H and  PHP  is  the  accuracy  of  approximate  wave 
functions  for  the  ground  and  excited  states.  We  have  reviewed  briefly 
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several  formulae  existing  in  the  literature  and  trivally  verified  that 
the  bounds  are  not  necessarily  very  good  in  a simple  test  case.  The 
formulae  require  lower  bounds  to  the  excited  states,  and  hence  this 
warrants  calculations  of  lower  bound  for  excited  levels. 

In  Chapter  III  we  show  an  explicit  solution  of  the  inhomoge- 
neous equations  that  arise  from  a truncation  of  the  reduced  resolvent 
T associated  with  the  unperturbed  Hamiltonian  H . It  provides  a solu- 

u O 

tion  to  all  orders,  avoiding  the  pitfalls  of  Brillouin-Wigner  low  order 
perturbation  theory.  The  truncation  is  related  to  the  variational 
principle  and  by  applying  successive  partitioning  technique  we  obtain 
the  Feenberg  perturbation  series  directly.  This  derivation  stresses 
tha  fact  that  the  Feenberg  expansion  is  a bracketing  function  in  the 
sense  of  Eq.  (1.12),  i.e.,  it  is  an  implicit  relation  = f(S),  and 
and  £ bracket  the  eigenvalue  E.  As  it  is  well  known,  ^6]  successive 
partitioning  technique  assumes  a particular  simple  form  for  triple 
diagonal  matrices. 

There  are  certain  problems  where  the  Hamiltonian  matrix  has 
that  form:  Mathieu  problem,  anharmonic  oscillator,  linear  chains  with 
nearest  neighbour  interactions,  and  any  finite  symmetric  matrix  can  be 
brough  to  triple  diagonal  form  in  a finite  number  of  steps.  We  have 
shown  that  the  recurrence  relation  derived  by  Nordling^0^  can  be 
expressed  as  a continued  fraction,  and  this  is  an  alternate  derivation 
of  the  continued  fraction  representation  of  the  solutions  of  the 
Mathieu  problem.  We  believe  that  this  may  be  exploited  computationally 
in  the  solution  of  eigenvalue  problems  and  in  the  treatment  of  systems 
with  nearest  neighbour  interactions. 


137 


• In  Chapter  IV  we  have  reviewed  the  theory  of  reciprocal 

variational  problems  and  how  it  may  be  applied  to  quantum  mechanical 
problems  through  inhomogeneous  equations  like  those  appearing  in 
perturbation  theory.  This  allowed  a new  derivation  of  the  lower  bound 
formulae  for  the  even  order  energies  in  Rayleigh-SchrBdinger  theory, 
originally  due  to  Hirschfelder  and  Prager  by  analogy  with  electromag- 
netic theory.  The  upper  bound  formulae  for  the  even  order  energies, 
in  both  Ray leigh-SchrBdinger  and  Brillouin-Wigner  theory,  due  orig- 
inally to  Hylleraas,  find  a natural  place  in  the  present  formalism, 
and  the  extension  to  excited  states  follows  naturally  from  this  essen- 
tially geometrical  approach. 

Chapter  V contains  a discussion  of  another  truncation  of  T , 

o 

which  leads  to  a set  of  partitioning  formulae  that  yield  explicit 
solutions  to  the  secular  equations  of  the  intermediate  problems.  In 
that  sense  they  are  formally  equivalent,  yet  we  should  point  out  that 
the  partitioning  solutions  allow  us  to  bypass  the  secular  problem. 
Furthermore  they  provide,  also,  a set  of  lower  bounds  to  the  eigen- 
values of  PHP . This  comes  about  by  searching  the  asymptotes  to  the 
bracketing  function.  This  important  byproduct,  useful  within  the 
context  of  Chapter  II,  is  not  obtained  in  the  conventional  method  of 
solving  the  secular  equation  of  the  intermediate  problems. 

We  should  emphasize  here  that  the  truncations  of  Chapter  III 
and  V are  not  the  only  tools  available  to  compute  the  bracketing  func- 
tions of  the  partitioning  technique.  The  multiple  dimensional  parti- 
tioning technique  combined  with  inner  projections  of  positive  operators 
has  proven  to  be  extremely  successful,  as  shown  by  Choi  He  applied 
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transformations  of  the  inner  projection  space  suggested  by  LBwdin  in 
PTXI,  and  obtained  very  good  results  for  He  like  ions  including  H~, 
which  is  a pathological  case  in  the  intermediate  Hamiltonian  approach 
due  to  the  ordering  theorem  1.2. 

The  eigenvalue  transformation  introduced  in  Chapter  VI  includes, 

as  particular  cases,  transformations  known  to  Schrbdinger^^  and 

r op  i 

Hylleraas  . . For  one-electron  problems,  it  yields  a complete  discrete 

spectrum  when  applied  to  a central  field.  We  showed  that  in  the  many- 
electron  case  the  one-electron  functions  depend  on  the  quantum  numbers 
of  the  set  of  spin  orbitals  defining  the  state,  as  stated  by 
Lundqvist^^  and  that  the  transformation  does  not  give  a complete 
discrete  spectrum  in  this  case,  but  it  yields  a continuum,  too. 

As  an  Hermitian  positive  operator  is  generated,  we  may  apply 
the  standard  techniques.  This  leads  to  several  formulae  for  lower 
bounds.  We  have  shown  how  the  method  is  applicable,  in  principle,  to 
molecules.  They  may  be  treated  in  several  ways,  one  of  which  is  by 
truncations  of  the  type  introduced  by  Bazley  and  Fox^^  for  the  calcu- 
lation of  lower  bounds  to  molecules.  The  main  difference,  and  it  is 
not  a trivial  one,  is  that  the  truncations  are  made  on  operators  with 
a complete,  discrete  spectrum.  It  was  drawn  to  our  attention  that 
calculations  on  H2+,  made  by  Walmsley/ fall  into  this  scheme.  If 

we  compare  them  with  an  analogous  truncation  method  for  H2+,  carried 

[73] 

by  Johnson  and  Coulson.,  , we  see  the  superiority  of  the  former 
approach,  and  it  can  be  ascribed  to  the  neglect  of  the  continuum  in 
the  latter. 
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We  have  discussed  a host  of  problems  that  arise  in  the  appli- 
cation of  partitioning  technique  to  the  calculation  of  upper  and  lower 
bounds  to  atomic  and  molecular  systems.  It  is  felt  that  much  work  has 
to  be  performed  in  the  future,  in  the  areas  of  properties  other  than 
energy  and  in  the  molecular  systems.  Further  development  of  the  formal 
inequalities  and  spectral  resolutions  is  needed.  The  inhomogeneous 
equation  approach  is  very  forceful,  and  some  systems  may  be  treated  by 
the  continued  fraction  approach  inherent  to  triple  diagonal  matrices. 

The  case  of  negative  perturbations  and  the  approach  to  molecular 
systems  that  follow  from  the  last  chapter  are  certainly  of  much 
interest  and  will  be  further  studied. 

After  completing  this  dissertation,  the  work  by  Swalen  and 

[9 1] 

Pierce  was  brought  to  our  attention.  Their  treatment  of  triple 
diagonal  matrices  leads  to  a continued  fraction  calculation  of  the 
eigenvalues,  and  in  fact,  the  results  are  identidical  to  ours.  On 
the  other  hand,  the  successive  pivotal  condensations  they  employ  or 
the  elimination  method,  can  be  cast  within  the  framework  of  succes- 
sive partitioning  technique,  and  this  is  the  basis  for  the  study  of 
the  convergence  properties,  and  it  would  be  valuable  to  repeat  and 
extend  their  discussion  taking  advantage  of  this  fact. 


APPENDIX  I 


In  this  appendix,  we  illustrate  some  of  the  properties  of  the 
bracketing  function  = f(£),  in  particular  the  convergence  as  deter- 
mined by  the  derivative  f 1 ( £.) . For  the  present  purpose,  it  is  conven- 
ient to  choose  a case  where  the  exact  solutions  are  known.  Thus  we 
choose  a simple  secular  equation  for  the  He  atom,  and  use  antisymme- 
tized  products  of  hydrogen-like  functions  as  a basis.  We  proved  in 

Chapter  II  that  the  partitioning  solution  is  given  by  the  truncation 

! 

T0(p)  as  defined  in  Eq.  (3-36)  and  Eq.  (3-37).  The  corresponding 
system  of  linear  equations  in  Eq.  (3*39)  can  be  solved  by  iteration. 

These  calculations  are  described  in  reference  [92] . It  is 
convenient  but  by  no  means  necessary  to  express  the  bracketing  func- 
tion corresponding  to  Eq.  (3*37)  in  terms  of  the  energy  shift  only, 
i.e.,  for  a state  g: 


£l,q  = Eq  + ^qlv"1-^)!^)  = E°  + ' 


(A.l) 


where 

P 0° 

Tc>>  £ - £°  (A-2) 

k=o  k 

k^q 

Then  by  applying  successive  partitioning,  we  may  obtain  an 
implicit  equation  for  t (3).  This  is  done  by  triangularizing  the 
system  of  equations  in  Eq.  (3-44),  which  in  fact,  yields  the  Feenberg 
approximant  for  a limited  basis  of  eigenfunctions  to  Hq,  as  we  have 
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shown  in  Chapter  III.  Given  t^(p)  “ f ( t 1 qq( P ) ) » the  derivative  f' 
follows  easily  and  for  the  case  of  p = 3,  it  has  the  form  which  we 
gave  in  reference  [92]. 


f'(£)  = [V2  V2  + V2  v2  + V2./2  + V2  72 
Jjm  qm  JJm  qjj'qm  qm7q$ 


2 n-2 


(A. 3) 


+ 2VqmVmiV^q^7q^+7qm^^7q^  + 7qm  " 


where 


j^q 


and  the  subindices  q,  2,  and  m may  have  the  values  1,  2,  or  3 with 
the  restrictions  that  q f i,m,  . 

Using  as  a basis  the  configurations  (is)2,  (ls)(2s)  and  (2s)2, 
with  scale  factor  z — 2,  which  amounts  to  having  as  a perturbation 

v = T^>  we  plotted  (A. 3)  as  a function  of  £.  for  q - 1,  2,  and  3. 
From  (A. 3)  it  is  easy  to  derive  the  condition  for  f1  to  blow  up.  It 
is,  with  the  same  injunctions  on  the  indices  as  in  (A. 3) 


+ V ' ‘ 0 • (A-M 

We  will  see  that  the  values  of  <£■  for  which  the  bracketing  function 
£l,q  blows  up  are  very  close  to  .the  values  = f(E^),  k^q,  i.e.,  to 
the  other  roots  of  the  secular  equation. 

The  convergence  condition  for  the  iterative  method  of  solution 
is  given  by  Eq.  (I.I5),  i.e.,  the  procedure  converges  if  f'(£)  <;  1, 
and  this  is  what  determines  the  selective  nature  of  the  partitioning 
solution.  Let  us  examine  (A. 4).  If  we  neglect  the  off  diagonal 
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element  with  respect  to  (7^  +7^m)2  >(A.4)  may  be  replaced  by 

7 nl  — 0 
qH  qm 


(A. 5) 


which  means  that  either  factor  is  zero.  Hence, 


£ a E°  + V? . — E 1 
J JJ  J 


is  the  approximate  value  for  which  Pi  blows  ud. 

^iq  e 


As  Ej  + V°^  is  an  approximation,  given  by  first  order  pertur- 
bation theory,  to. the  eigenvalue  of  the  secular  equation  we  may 

expect  that  the  eigenvalues  E^»  k^q  are  close  to  the  points  at  which 
q l3  lows  up,  and  indeed  this  is  what  happens.  From  reference  [92]  > 
we  obtain  the  following  results. 

Eigenvalues  of  secular  equation  for  He  atom,  using  the  config- 

p p 

urations  (is)  , (s)(2s)  and  (2s)  as  a basis  (unsealed): 

= -2.830996 

e|  = -1.9569^8 
e'2  = -0.697606 

Values  of  £ that  make  blow  up: 

i>q 


For  q = 0: 
For  q = 1: 
For  q = 2: 


-2.036792, 

-2.750939 

-2.830438 


-O.698778 

-O.69828 

-1.9559134 


For  comparison  we  give  the  first  order  perturbation  energy 
approximation:  E ' E°  + V°  : 

q q qq 
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-2.750000 

E^  -2.036351 

-O.6992I9 

This  simple  study  of  the  convergence  of  the  iterative  proce- 
dure gives  us  the  conditions  under  which  we  may  expect  ^ 1 to  converge 

co  E - j-(E  ) and  not  to  ciny  other  root  of  the  secular  equation  E1  . 

H H lC 

More  details  are  given  in  reference  [92]. 
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Fig.  1.  Bracketing  Theorem 


1*5 


£ 

& 

-4- 

v 

'i  u) 


Fig.  2.  Bracketing  Function 
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Typical  Term  in  Br illouin-Wigner  Series 
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Fig.  5 . Reciprocal  Variational  Problems 
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Fig.  6.  Partitioning  Function  for  Truncated 
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Fig.  7.  Eigenvalue  Relation  X(£) 
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Fig.  8.  Upper  and  Lower  Bounds  to  \(d) 
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